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I. INTRODUCTION 



Data links used in many modem military communication systems re- 
quire computer-to-computer communications which transmit bits from one 
terminal to another with a low error rate. The communication channel may 
be a lowpass channel, if the transmission medium is wire or optical fiber, or a 
bandpass channel, if the transmission medium is radio frequency (RF) link or 
acoustic channel. The signal used to carry the data bits must be a properly 
modulated signal with sufficient energy. It must also be positioned in the fre- 
quency spectrum to propagate effectively through the intended medium. 
Multi-Frequency Modulation (MFM) is a signal modulation format that is 
readily adaptable to a variety of data link scenarios. MFM does not require 
special purpose MODEMS to translate between the digital and analog do- 
mains and can emulate most existing signal modulation formats as well as 
generate entirely new formats. MFM’s descriptive language is that of Digital 
Signal Processing (DSP) [Ref. 1]. Simulation and analysis of a Naval Post- 
graduate School (NPS)-developed MFM signal received through a bandpass 
channel is the subject of this thesis. 

A. BACKGROUND 

Classical modulation methods for bandpass channels use amplitude and/ 
or phase to carry signal information on a carrier wave in the channel. When 
the information source is a finite alphabet, as with data or quantized analog 
sources such as digitized speech or video, then only a finite number of signal 
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states are needed to represent or code the source. Phase shift keying (PSK) is 
an efficient method for coding these signal states or message [Ref 2]. 

MFM is a modulation technique that is ideally suited for computer-to- 
computer data transmission and reception because its basic structure is one 
of time and frequency slots. In MFM, the signals are directly encoded, modu- 
lated, decoded, and demodulated using DSP techniques within the host com- 
puter. In MFM, signal “packets” located in a given frequency spectrum and 
time are created as shown in Figure 1. Figure 1 shows a single signal packet. 
Each packet consists of L bauds of K tones. The length of a baud is AT seconds 
during which K discrete tones are transmitted. The phase of each tone is the 
coded information. These LK subsignals form an orthogonal signal set. Each 
subsignal may be independently modulated with phase information. 

Dr. P. H. Moose at NPS, in conjunction with the Naval Ocean Systems 
Center (NOSC) in San Diego, has developed Multi-Frequency Quadrature 
Phase Shift Keying (MFQPSK) signals to be used in high data-rate acoustic 
burst communications from moving platforms using medium frequency, rela- 
tively wideband (25 percent bandwidths) links [Ref. 1, 3]. A simulation of the 
MFQPSK signal from a moving platform as seen through a bandpass channel 
was needed to analyze and test various Doppler, synchronization, and coding 
techniques and/or algorithms before final implementation of the modulation 
scheme is actually realized. Thus, the subject of this thesis is the simulation 
and analysis of a MFQPSK signal from a moving platform as seen through a 
bandpass channel. 
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Figure 1. MFM Signal Packet (after Ref. 1; p. 3.) 



B. OVERVIEW 

The NPS-designed MFQPSK signal will be transmitted from a moving 
platform. In Chapter II a simple model of the moving transmitter platform’s 
dynamics is developed; its performance is compared to actual track data. This 
model is implemented in the simulation. 

A mathematical description of the transmitted signal and the design 
parameters of the MFQPSK signal are introduced in Chapter III. The channel 
parameters which delay and compress/expand this transmitted signal and 
which characterize the received signal, the additive Gaussian noise which is 
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added to the received signal after sampling, and the final simulated sampled 
received signal are derived and presented in Chapter III. 

Once the simiilation was complete, testing and analysis was performed. 
First, the simulation was tested against MFQPSK signal theory to insure 
that the simulated received MFQPSK signal is consistent with theory. The 
simulated received signal should be consistent with theory so that the simu- 
lation may be used as an experimental tool for testing various Doppler, syn- 
chronization, and coding algorithms. Chapter IV consists of this analysis and 
simulation results. 

Second, the amount of residual Doppler mismatch which the signal can 
tolerate due to the moving transmitter was analyzed. Chapter V provides the 
analysis of the output signal-to-noise ratio degradation due to the Doppler 
mismatch. Then, in Chapter VI, an estimation of the Doppler compression 
factor within a baud is derived using the Discrete Fourier Transform of a re- 
ceived signal’s baud. The simulated outcome of this estimated Doppler versus 
the Doppler due to the moving transmitter within a baud is also illustrated in 
Chapter VI. 

In the simulation, the actual Doppler compression/expansion factor due 
to the dynamics of the moving transmitter is computed; therefore, it is 
known. The start time of the signal reception is also known; thus, the signal 
is always perfectly synchronized in the simulation. 
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II. MOVING TRANSMITTER PLATFORM DYNAMICS AND 
BEAM PATTERN 

A. BACKGROUND 

Since it is assumed that the MFQPSK signal will be transmitted from a 
moving platform, provision must be made at the receiver to compensate for 
Doppler shifts. The Doppler shift factor is computed based on the geometry in 
the simulation using the radial velocity of the moving transmitter relative to 
the stationary receiver. The complete set of equations for the Doppler shift 
factor is presented in the next chapter. The radial velocity is the derivative 
with respect to time of the slant range to the receiver. Realistic slant range 
values are necessary to compute realistic Doppler shift factors. Slant range is 
computed using the time-varying x-position, y-position, and z-position of the 
transmitter relative to the receiver. 

B. ACTUAL TRACK DATA VS. SIMULATED TRACK DATA 

To realistically model the dynamics of the moving transmitter, actual 
track data was provided by NOSC. The parameters that were given in the 
actual track data were x-position, y-position, z-position, and slant range rela- 
tive to the receiver. A simple random walk model was used to model the track 
data. This model assumes the transmitter is moving along a straight line 
with random fluctuations in the x, y, and z velocity components. Figures 2 
through 7 illustrate how well the random walk model corresponds to the 
actual track data on six different runs. The z-position seems to consistently 
have the largest error. The largest error in the z-position is 25 feet, (see 
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Figure 2a. Run 1: Actual and Simulated Track Data 
X— Position and Y— Position vs. Time 
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Figure 2b. Run 1: Actual and Simulated Track Data 
Z-Position and Slant Range vs. Time 
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Figure 3a. Run 2: Actual and Simulated Track Data 
X— Position and Y— Position vs. Time 
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Figure 3b. Run 2: Actual and Simulated Track Data 
Z-Position and Slant Range vs. Time 
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Figure 4a. Run 3: Actual and Simulated Track Data 
X-Position and Y-Position vs. Time 
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Figure 4b. Run 3: Actual and Simulated Track Data 
Z— Position and Slant Range vs. Time 
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Figure 5a. Run 4: Actual and Simulated Track Data 
X— Position and Y— Position vs. Time 
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Figure 5b. Run 4; Actuad and Simulated Track Data 
Z-Position and Slant Range vs. Time 



13 





Figure 6a. Run 5: Actual and Simulated Track Data 
X-Position and Y-Position vs. Time 
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Figure 6b. Run 5: Actual and Simulated Track Data 
Zr-Position and Slant Range vs. Time 
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Figure 7a. Run 6: Actual and Simulated Track Data 
X— Position and Y— Position vs. Time 
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Figure 6b). Note that the z-position graphs have a smaller vertical scale than 
the X- and y-position graphs. Ideally, since the transmitter is attempting to 
maintain a constant depth, the location should not fluctuate much in the 
z direction as was the case in Run 5 shown in Figure 6b. However, overall 
there is little error for these six runs in the model output for the slant range; 
therefore, the random walk model was considered to be an adequate model of 
the transmitter’s dynamics for the purpose of simulating Doppler shifts, time 
delays, and propagation losses of the MFM acoustic communication signal. 

C. TRANSMITTER PLATFORM DYNAMICS MODEL 

The following three equations are the random walk model used to simu- 
late the transmitter’s dynamics. The names of the variables used in the simu- 
lation are the same as those used in the equations below. When applicable, 
the parameters from Figure 1, which are printed in boldface, are equated to 
their simulation variables, printed in capital letters. The d)mamics model is 
described by; 

X(LL) = X(LL-l) -h (VX(LL) * DELT(LL)) (1) 

Y(LL) = Y(LL-l) -t- (VY(LL) * DELT(LL)) (2) 

Z(LL) = Z(LL-l) -h (VZ(LL) * DELT(LL)) . (3) 

The slant range to the receiver is computed using the following equation; 

R(LL) = ( X(LL)2 -h Y(LL)2 -h Z(LL)2 )0-5 (4) 

where, 

X(LL): The position of the transmitter relative to the receiver in the x 

direction during the baud 

Y(LL): The position of the transmitter relative to the receiver in the y 

direction during the baud 
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Z(LL): The position of the transmitter relative to the receiver in the z 

direction during the baud 

R(LL): Slant range of the transmitter relative to the receiver during 

the LL^h baud 

LL = 1: The baud number of the transmitted signal 

DELT(LL) = AT : The baud length in seconds 

VX(LL): The velocity in the x direction during time DELT(LL). VX(LL) is 
Gaussian distributed with mean = VXAVG and variance = 
VXVAR. VXAVG and VXVAR are input by the user. 

VY(LL): The velocity in the y direction during time DELT(LL). VY(LL) is 
Gaussian distributed with mean = VYAVG and variance = 
VYVAR. VYAVG and VYVAR are input by the user. 

VZ(LL): The velocity in the z direction during time DELT(LL). VZ(LL) is 
Gaussian distributed with mean = VZAVG and variance = 
VZVAR. VXAVG and VXVAR are input by the user. 

Note: LL and 1 are used interchangably throughout. 

The rectangular coordinate system of the transmitter’s motion is relative 
to the receiver. The origin (0, 0, 0) is perpendicular to the receiver in the 
plane of the transmitter (the x-y plane) as shown in Figure 8. 

D. TRANSMITTER BEAM PATTERN 

If the receiver is not within the transmitter’s transmission beam, then 
the message will not be received. The simulation does not use the actual 
beam pattern of the transmitter; therefore, it does not know if the receiver is 
within the beam pattern. The simulation assumes the receiver is within the 
transmitter’s beam the entire time of transmission (i.e., from the beginning of 
the first signal packet to the end of the last signal packet). 
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Figure 8. Coordinate System of the Transmitter and Receiver 



The angle between the transmitter’s slant range and z-position relative 
to the receiver is computed and output graphically by the simulation. This 
angle, THETAD(LL) in the simulation, is 

THETAD(LL) = arccos( abs(Z(LL)) / R(LL) ) * (180.0 / PI) (5) 

where PI = 3.141592654. Figure 9 illustrates the geometry of the line of sight 
angle, THETAD(LL) = 0, and THETAO = 0Q, the half beam width of the 

transmission beam. 

The simulation user can graphically observe the value of THETAD(LL) 
for each baud and determine for a given THETAO if the receiver would 
actually be able to receive the transmitted signal. The simulation does not 
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stop if THETAD(LL) is greater than the input half beam width of the trans- 
mitter, THETAO, but a warning message is output to the screen. The initial 
position (XO, YO, ZO) and THETAO are inputs to the simulation. 




Figure 9. Geometry of the Transmission Beam 
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III. DERIVATION OF THE SIMULATED RECEIVED 
MFQPSK SIGNAL 

A. DESCRIPTION OF THE TRANSMITTED MFQPSK SIGNAL 

Referring to Figure 1, the following definitions are used in MFQPSK 
[Ref 1]; 

T : Packet length in seconds 

AT : Baud length in seconds 

kx : Baud length in number of samples (not in Figure 1) 

L ; Number of bauds per signal packet 
At ; Time between samples in seconds 
fx = 1 / At : Sampling frequency in Hz 

Af = 1 / AT : Minimum frequency spacing between MFM tones in Hz 
K; Number of MFM tones in a baud 
1 : Baud number 

k : Harmonic number of the MFM tone 

0|]^ : Symbol or phase on the k^k tone of the pk baud 

Since At = AT / kx, the sampling frequency is fx = kx * Af. Consequently 
there are a maximum of kx/2 tones spaced Af Hz apart in frequency covering 
the range from 0 Hz to fx/2 Hz, Here fx/2 is the Nyquist frequency [Ref. 2]. 
The K tones carry the phase information during each baud. The relationship 
between fi, the i*^k frequency, and kj, the i^^ harmonic, is kj = fi / Af. Some of 
the tones may not be used (i.e., their amplitudes are set equal to zero) during 
any or all bauds of the packet. To generate the transmitted bandpass signals 
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between frequencies fi and ( 2 , only tones between harmonics ki = fi / Af, the 
minimum harmonic number or tone, and k 2 = f2 / Af, the maximum tone, will 
be given non-zero amplitudes. The harmonic numbers less than ki and 
greater than k 2 are given zero amplitudes. The number of tones in a baud is 
K = k 2 - ki -I- 1. These K contiguous tones are transmitted with non-zero am- 
plitudes. The signal bandwidth is W = K * Af. Thus, the time bandwidth 
product of the entire signal packet is TW = L * AT * Af * K = LK, which is the 
total number of symbols that can be sent in one signal packet. 

A mathematical description of the transmitted MFQPSK signal is neces- 
sary to understand some of the parameters used to simulate the received 
MFQPSK signal. The 1^^ baud of the transmitted signal is described by: 

^2 

Xj (u) = ^ (u) ; 0 < u < AT (6) 

k=k^ 

where xu^(u) = cos( 2 ji k Af u + ) [Ref. 1]. Here, u is time referenced to 

the beginning of the baud. Actual real time is t = t^ + (1*AT) + u where t(j is 
the time at the beginning of the 0^^ baud (i.e., the beginning time of the first 
signal packet). 

The discrete time signal corresponding to the baud is generated by 
sampling (6) at the sampling intervals At = l/f^. Thus, the discrete time signal 



where xn^Cn) = cos( (2k k n) / kx + ). Here, n is discrete time refer- 

enced to the beginning of the baud. 

Note that a baud interval of time AT seconds contains exactly k cycles of 
tone k. Therefore, adjacent tones differ by one in the number of cycles they 
generate during a baud. 

The phase, Oj]^, may be given one of the four values k/4, 3k/4, -3k/4, or 
-k/4, which are in quadrants 1 through 4 respectively. These four values make 
up the symbol set used to code the information or message on the signal. The 
simulation variable for is PHI(LL,K). 

The design parameters used in the MFQPSK signal designed and devel- 
oped by Dr. P. H. Moose at NPS are listed in Table I. These parameters are 
for a signal packet in a 16 to 20 KHz bandpass channel. These parameters 
are also used in the simulation to uniquely characterize a chosen baud type 
(i.e., baud types 1 through 5). The simulation variables in capital letters are 
equated to the parameters they represent. Recall that LL = 1, the baud 
number. 

B. TIME DELAY AND COMPRESSION/EXPANSION OF THE 
TRANSMITTED SIGNAL 

When the MFQPSK signal is transmitted, there are various factors 
which can affect the signal such as a moving transmitter, the channel or the 
medium, and the receiver. In the simulation, these parameters are computed 
and applied to the transmitted signal, thus producing a model of the received 
signal baud. By delaying and compressing/expanding the signal in time, the 
frequencies of the transmittedsignal are shifted in the received signal. Some 
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TABLE I 



DESIGN PARAMETERS FOR SIGNAL PACKET 
IN A 16-20KHZ BANDPASS CHANNEL 



BAUD TYPE: 


1 


2 


3 


4 


5 


AT = DELT(LL) 


1/240 


1/120 


1/60 


1/30 


1/15 


(sec) 


Af=DELF(LL) 


240 


120 


60 


30 


15 


(Hz) 


ki = KMIN(LL) 


68 


135 


269 


537 


1073 


fl (Hz) 


16320 


16200 


16140 


16110 


16095 


k£ = KMAX(LL) 


83 


166 


332 


664 


1328 


f£ (Hz) 


19920 


19920 


19920 


19920 


19920 


kx = KX(LL) 


256 


512 


1024 


2048 


4096 


fx (Hz) 


61440 


61440 


61440 


61440 


61440 



of the parameters or variables, which delay and compress or expand the 
transmitted signal in time, are illustrated in Figure 10 for the 1^^ baud. 

Recall xj(u) is the baud of the transmitted signal, which is described 
by (6). In Figure 10, X](u) and yi(u) are represented with a rectangle for illus- 
tration purposes only. Recall xj(u) is the superposition of xij^(u), for ki<k^2» 
over a total of K tones. The frequency of xjk(u) is 

coj^ = cOx * (k / kx) (8) 
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\ 




Figure 10. The Baud of the Transmitted Signal Delayed and 
Dilated in Time 



where cOx = 2n and fx is the sampling frequency of the transmitted signal. 
Table I lists the values of these parameters and the applicable variables 
which are used in the simulation. 

When xjij(u) propagates through the transmitter system, it may be 

delayed by the transmitter system electronics. In the simulation, this delay is 
denoted by = TAUK(K), and is a constant for all k. (In the simulation, 

TAUK(K) may be set to a constant or to any known function of k.) Acoustic 
transmission of xjCu) begins at u = and ends at u = AT + Tj^. 

When the signal, X](u - T(^), leaves the moving transmitter, it is delayed 
and compressed or expanded in time due to the movement of the transmitter 
and the distance it travels through the channel to the receiver. (The time 
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delay introduced by the electronics of the receiver has not been included in 
the simulation.) This time delay and compression/expansion factor will be re- 
ferred to as £(u). Derivation of £(u) is presented in Appendix A. The time 
delay and compression/expansion factor due to the medium is described by 

£(u) = (T] + ttju) / (1 + aj) (9) 

where 

X] = The time for the baud of the transmitted signal to travel through 

the channel and arrive at the output of the receiver. 

ai = The Doppler compression factor due to the moving transmitter 
platform. 

The parameters, X\ and ai, depend on the slant range to the receiver, 
R(LL), and the speed of sound of the acoustic channel, C(LL), for the 1*^^ baud. 
Recall R(LL) is computed using equation (4). The speed C(LL) is Gaussian 
distributed with mean = CO and variance = CVAR. The parameters CO and 
CVAR are input by the simulation user. The parameters, Xj and aj, are 

related to R(LL) and C(LL) in the following manner: 

XI = TAUL(LL) = R(LL) / C(LL) (10) 

and 

ai = ALPHA(LL) = R(LL) / C(LL) (11) 

where 

[(X(LL) * VXAVG) -t- (Y(LL) * VYAVG) + (Z(LL) * VZAVG)] 

R(LL) (12) 

Recall X(LL), VXAVG, Y(LL), VYAVG, Z(LL), and VZAVG were described in 
the previous chapter with equations (1) through (3). TAUL(LL) and 
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ALPHA(LL) are the simulation variables for Tj and aj, respectively. When the 
transmitted signal, X|(u - Tj^), arrives at the receiver output, it is compressed/ 
expanded in time by £ (u). 



The received signal is 




(13) 



where 



yik(u) = Fik A]i^ cos(2tc k Af (u - - £(u)) + Oi^), 



(14) 



and is the attenuation of tone k due to the propagation through the 
channel. The signal, yi(u), will begin at time uj and end at time U 2 . The 
signal, xi(0), transmitted at time u = 0 is the same signal, yi(ui), that arrives 
at the receiver at time u^; and likewise, the signal, xi(AT), transmitted at time 
u = AT is the same signal, yi(u 2 ), which arrives at the receiver at time U 2 (see 
Figure 10). Since it is the same signal, the times that this signal arrives at 
the output of the receiver may be equated to the corresponding times that 
this signal was sent to yield 



It is obvious that yj(u) is also a superposition of signals, yik(u). However, 
yjk(u) is at frequency 0)^', where 



ui = (1 + aj) Tk + 'T] 



(15) 



and 



U2 = (1+ «i) AT + (1 + tt]) ^k + X] . 



(16) 



COk' = 0>k / (1 + «]) 



(17) 
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( 18 ) 



due to the moving transmitter. At the receiver, y](u) will be sampled at 

0)y = C0x/(1 + Om) 

where is the Doppler factor associated with the m^^ Doppler channel. If 
tt] = (Xni, then the frequencies of the received signal are exactly in the center of 
Doppler channel m. However, in general, aj will fall between Doppler 
channels and there will be some residual Doppler mismatch a] - ajj,. If the 
channels are spaced at intervals Aa, then the maximum mismatch will be ± 
Aa/2. In the follovdng section, a derivation for the spacing of the Doppler 
channels Aa is presented. 

C. DOPPLER ESTIMATION AND DOPPLER CHANNELS 

Several processing approaches were considered for estimating the chan- 
nel model parameters. However, the ultimate goal is to send a signal through 
a bandpass channel using MFQPSK modulation on a succession of bauds 
which constitute a signal "packet". In the presence of white noise, the opti- 
miun receivers for these signals are filters matched to each of the tone/phase 
combinations. Recall the four phases are {ji/4, 3 k/ 4, -3ro'4, -5 t/4) for a MFQPSK 
signal. Therefore, a filter matched to the transmitted signal with = ti/4 

will have a positive output when Ji/4 is sent, a negative output when -3ti/4 is 
sent, and an output of zero when 3ji/4 or -ti/4 is sent, provided that the 
receiver is synchronized with the signal. Similarly, the output of a filter 
matched to the transmitted signal with phase equal to 3ti/4 is positive when 
3tc/4 is sent, negative when -ti/4 is sent, and zero otherwise. Also filter pairs 
matched to the i*^^ frequency of the transmitted signal, fj = kjAf, produce zero 
outputs for all phases at the frequency fj = kjAf when kj kj. 
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To summarize, a matched filter system for MFQPSK signals with K 
tones and baud length AT = 1/Af consists of K pairs of filters, one filter of each 
pair matched to phase niA and the other filter of the pair matched to phase 
3 tc/ 4. Each filter pair output demodulates the phase information encoded on a 
particular tone. 

Since the maximum Doppler shift will be present on the highest fre- 
quency or tone in the baud, the response of the filter pair matched to this tone 
will be analyzed to show how the Doppler channel spacing, Aa, is derived and 
how Oj„ is computed in the simulation. 

Let hi(u) and hq(u) be the filters matched to the signals: 

xik 2 ,i(u) = A]k2COs(2K lc 2 Af u + ti/4) ; 0 < u < AT (19a) 

and 

xjk 2 ,q(u) = A|k2COs(2;i lc 2 Af u + 3 tc/4) ; 0 < u < AT. (19b) 

where the indices i and q denote in-phase and quadrature phase, respectively. 
Thus, hj(u) and hq(u) are described as follows: 

hj(u) = 2 cos(2;: k 2 Af (AT - u) + n/4) ; 0 < u < AT (20a) 

and 

hq(u) = 2 cos(2rt lc 2 Af (AT - u) -i- 3 jc/4) ; 0 < u < AT. (20b) 

Now the received signal from equation (14), given 7t/4 is the phase of the tone 
k 2 , is 

yik2^u) = Fik2Aik2COs(2j:k2Af(u-Tk2-£(u)) + 7c/4);ui<u<U2 . (21) 

The received signal, given 3;i/4 is the phase of the tone k 2 , is similar to (21) 
with <3>ik = 3x/4 instead of ti/4 as above. 
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As illustrated in Figure 11, the output of the matched filter pair, in gen- 
eral, is the convolution of h(t) and yik(t) given by 



and 



X 

max 

0 



X 

max 

Jh^(t)yn(u-t)dt. 

0 



(22a) 



(22b) 



where Xj^ax = ^2 " ^i- Evaluating (22a) and (22b), the following expressions 
are obtained for the output of the matched filter pairs due to tone k 2 : 



X 

max 



^Ik, 



i = ^Ik 



lk„ 



J* COS [2;:k2Af(u - £(u - x) - Xj^ - AT)] dx 



(23a) 



and 



max 



3k, 



(u) = F„ A,, f sin [2;ik Afiu - £(u - x) - x, - AT)] dx 



(23b) 
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Figure 11. Matched Filter Pair for Tone k 2 



Inserting equation (9) for £(u) in the above two equations yields 

(24a) 



(24b) 

Recall that the received signal arrives when u^ < u < U 2 with 

ux = [(1 + «]) + Tj] and U 2 = [(1+ aj) AT + (1+ aj) + Tj] . 

Also recall from the previous section that aj is the Doppler compression/ 
expansion factor due to the moving transmitter. Note that for aj < 0, Tn^ax = 
U 2 - Ux ^ AT and for aj > 0, Tj^ax = ^2 - ux > AT. Using X^ax = ^ (1+oc]) for 

the upper limit in the integration in the convolution integrals (24a) and (24b), 
will not affect the integrals significantly because aj is of the order 10'3; 
therefore, diriding by (1+ aj) is very close to diriding by 1. 
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Sampling Z]i^ 2 ,i(u) and at u = U 2 and evaluating their integrals 

in (24a) and (24b) from 0 to yields 



u=u< 






AT (l+ttj) 



sin [27tk2(Xj] 

27tk2aj ^^2’^ 



(25a) 



and 



^k„,q I u=u - ^lk„^lk„ 



(1-COS [27ik20j]) 
27ck2(Xj 



^k2,q 



(25b) 



The total power, P]]^, in the tone of the 1^1* baud is one-half the ampli- 

tude, F]]^Ai]^, of the received signal squared. That is 



Flk^ik = (2 Pik)®-^- 



(26) 



Also assume l-t-aj = 1; therefore, AT (l-^al) is approximately equal to AT. With 
these substitutions, equations (25a) and (25b) become 



a c (sin [2jik«o,]) 
Z,, .(a,) = (2R, ^ ‘ 

^^2’' *^2 2;tk2aj 



and 



0.5 






(l-cos [ 2 jik 20 j]) 

27ck2aj 



(27a) 



(27b) 



A graph of Zi]^ 2 ^i(oc]) and Zu^ 2 ,q^^l^ from (27a) and (27b), respectively, is shown 
in Figure 12. 
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Figure 12. Matched Filter Outputs Versus Doppler Shift 



Note that at U] = 1/(4 k 2 ), I I = I Zii^ 2 ,q I ; this is illustrated in Figure 

12. This means that the response of the in-phase and quadrature channels 
are equal even though an in-phase tc/ 4 symbol was sent; therefore, a signifi- 
cant decoding error would be made if aj = 1/(4 k 2 ). It is obvious that, at aj = 0, 
Zjk 2 ,i( 0 ) is maximum and Zjk 2 ,q( 0 ) is zero; in which case, no decoding errors 
would be made. The ratio of to '^]k 2 ,q below, evaluated for various aj's 

less than 1/(4 k 2 ), produces the maximum width of the Doppler channels, Aa. 
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The ratio is described by: 




(28) 



For tt) = 1/(16 k 2 ), the ratio above is > 5.03 which is 14 dB. This is called inter- 
symbol separation. 

The value of 14 dB was considered the minimum amount of inter-symbol 
separation to be tolerated. Recall that the maximum amount of Doppler 
mismatch is ai-am=±Aa/2. Therefore, Aa, the Doppler channel spacing, 
is equal to 1/(8 k 2 ) and the magnitude of the maximum amount of 
residual Doppler mismatch is Aa/2 = 1/ (16 k 2 ), where ttm is the Doppler 

factor associated with the m^k Doppler channel. In the simulation, 
Aa = DELALF(LL) = 1/(8*KMAX(LL)). 

Solving tt] - Om = ± Aa/2 for a^, yields 



where aj„, the Doppler factor due to the m^l^ Doppler channel, is equal to 
m Aa. Substituting am=m Aa into (29) yields the following expressions for the 
lower bound, mj, and upper bound, m 2 , of m, the Doppler channel number as 



am = aj ± Ao/2 , 



(29) 



mi = (aj / Aa) - 0.5 



(30a) 



and 



m 2 = (aj / Aa) + 0.5 . 



(30b) 
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Since m is an integer between and m 2 , it is computed as 

m = int[(mi + m2)/2]. (31) 

In the simulation, m = M(LL) and is computed using (31); and 

= ALPHAM(LL) = M(LL) * DELALF(LL). 

D. THE ATTENUATION FACTOR DUE TO AN ACOUSTIC 
CHANNEL 

When the transmitted signal, X||j(u - Tj^), travels through the medium or 
channel, its amplitude will be attenuated by a factor Fjj^ due to the channel. 
In the simulation, the user has the input option to have this attenuation 
applied to the received signal or to let all the received signal tones have equal 
amplitudes. The quantity AA(LL,K) is the amplitude of the baud and 
tone of the simulation's received signal. If the user chooses to have equal 
amplitudes, then AA(LL,K) is set equal to 1 for all LL and all K. If the user 
desires the attenuation to be put on the received signal, then AA(LL,K) = Fjj^, 

which is computed below. 

In general, the transmission loss depends on frequency, depth, pressure, 
and temperature. However, when specific propagation conditions are of no 
interest and only a rough approximation of the transmission loss is adequate; 
the spherical-spreading law plus an added loss due to absorption is equal to 
the transmission loss [Ref. 4]. Therefore, the transmission loss in dB is 

TL = (20 logio R) -I- p R (32) 
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where 



R is the depth of the receiver from the surface in feet and 
P is the attenuation coefficient due to absorption in dB per foot. 

Thus, the following expression for Fij^ is used in the simulation: 

Fik = 10 -(TIV20). (33) 



The attenuation coefficient, P, is described by [Ref. 4] 



P = 



0.1 + (0.000275) + 0.003 






4.100 






(3000) 



J 



(34) 



where 

f = (k Af) / 1000 in kilohertz. 

This expression applies for a temperature of 39®F (4^0, which is an average 
seawater temperature, and a depth of about 3000 feet. The constant 0.003 is 
added to take care of the attenuation at very low frequencies [Ref. 4]. In the 
simulation, ABSORP = p. 

The absorption of seawater decreases by about 2 percent for every in- 
crease of 1000 feet in depth [Ref. 4]. Note in the simulation, the transmitter is 
at a depth of 1000 feet below the ocean's surface. This depth is set with the 
simulation variable TXDEP. Thus, the receiver is R = TXDEP + Z(LL) below 
the ocean's surface. TXDEP is not an input of the simulation; therefore, to 
change TXDEP from 1000, one must change TXDEP in the simulation code. 
The simulation computes P using equation (34), then increases P by 2 percent 
for every 1000 feet that the receiver's depth exceeds 3000 feet from the 
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ocean's surface or decreases p by 2 percent for every 1000 feet that the 
receiver's depth is less than 3000 feet from the ocean's surface. 



E. THE SAMPLED RECEIVED SIGNAL WITHOUT NOISE 

The received signal, yj(u), is sampled at an interval 



Au = [AT (1 + (Xm)] /kx . 



(35) 



(Recall cOy = 27T / [(1 + aj^) AT], (18), is the sampling frequency.) The 

parameters AT and kx are defined in section A. The Doppler factor, otjn=m Aa, 



trated in Figure 13. The synchronization error of the baud is Aui=ui-ui. 
The simulation variable names for the above parameters are Ui=Ul, 
ui=UHATl, and Aui=DELUl. 

The sampled received signal for the l^b baud is (13) with u = nAu + Ui. 
When uj = uj, the sampled received signal for the 1*^^ baud is (13), 



yjl^(nAu + Uj) = F]]^Aij^cos[27ik Af (n Au - (£(n Au + Uj) - £(uj))) + (37) 

Substituting (9) for £(u) and computing £(n Au + Uj) - £(u^) yields 



associated with the m*^^ Doppler channel was derived in the previous section. 
The first sample will be taken at u^, which is the best estimate of u^, the be- 
ginning or synchronization point of the 1*^1^ baud at the receiver. This is illus- 




k=k 



; 0 < n < (k^ - 1) 



(36) 



where 



£(n Au + uj) - £(uj) = (aj n Au) / (1 + aj). 



(38) 
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^ ^ y,(n) 





SAMPLED 
RECEIVER OUTPUT 

^ n ^ 



Figure 13. The Sampled Receiver Output 



Now substituting Au = [AT (1 + a^)] / into (38) produces 

£(n Au + ui) - £(ui) = [a] n AT (1 + am)] / kx (1 + a]). (39) 

Collecting the results from (39) and nAu = [n AT (1 + ttm)] / k^ together and 
substituting them into (37) yields the following sampled received signal for 
the baud and kth tone: 

yik(n) = FikA]]^ cos[(27rk/kx) ((l+ttm) / d+ai)) n + <J>id . (40) 

If a factor V]j^ is added to n to account for random timing jitters of the 1^^ baud 
and kth tone due to synchronization error, then (40) becomes 

yik(n) = FikAik cos[(2jik/kx) ((1+0^) / (!+«])) (n - vjk) + dD 

for 0 < n < (kx - 1) , 
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where 



^lk = 



Au, k 
1 X 

AT(l+a ) 

TTl 



Finally sununing (41) over all the tones in the baud gives 




; 0 < n < (k^-1) • 



( 42 ) 



(43) 



Summarizing, (41) is the final computational form of the received signal 
which is denoted by YY in the simulation. The received signal, YY, is then 
summed over all the k's for each baud from KMIN(LL) to KMAX(LL) to yield 
YRX in the simulation, which is yj(n) above. In the simulation, adl the packets 

are contiguous. The continuous sampled received signal is YYRX(II) where II 
ranges from 0, which is the first sample of the first baud (LL=1) in the first 
packet, to NPTS, which is the last sample of the last baud (LL = BDTOTL = 
the total number of bauds) in the last packet. YYRX(II) is equal to YRX plus 
NOISE with all the packets contiguous and the bauds within each packet 
contiguous, where NOISE is white Gaussian noise added to the received 
signal. NOISE, which depends on the input signal-to-noise ratio (SNR), is 
described in the next section. 



F. SIGNAL PLUS ADDITIVE NOISE 

Inevitably noise will be added to the signal, either from the environment, 
the electronics, or both. This noise is assumed to be additive white Gaussian 
noise. Since it is assumed that the received signal, yj(u), has been ideally 
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bandlimited to one-half the sampling frequency, there is no power in 
frequencies greater than or equal to iJ2. Recall that the signal is sampled at 
the transmitter at f^. If w(u) is white noise and has the white power spectral 
density equal to N(/2, then 

E [ w(u) ] = E [ w(n) ] = 0 (44) 

and 

o2 = VAR[ w(u) ] = VAR[ w(n) ] = Nq (45) 

where w(n) is the white noise sequence [Ref. 1]. 

Let the receiver input SNR be defined as the signal power in bandwidth 
W divided by the noise power in bandwidth W. The average tone signal power 
in the l^I* baud is defined as 

PAVG, = (1/K) 

k=k, 

^ (46) 

where Pii^ is the total power in the tone of the baud. (Recall that K is 

the total number of MFM tones in the baud.) Therefore, the wideband 
input SNR for the P^ baud [Ref 1] is 

SNRWB] = (K PAVGi) / (W Nq) = PAVG] / (Af Nq) = PAVG] k^ / (2 c\^), (47) 

because Af = f^ / and N^ = (2 a]2) / fjj. The narrowband input SNR for the 
pl’ baud is 

SNRNB] = Pii, / (Af No) = P]k kx / (2 ai2). (48) 

Note that if all the tones have equal amplitudes (i.e., AA(LL,K) = 1), then 
PAVG] is equal to Pjj^; in which case, SNRWBi = SNRNBj. 
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Solving (47) for q \^ yields 



^2 

of = [(1/K) X P]k^ \ I (2 SNRWBj) . 

(49) 

Recall Pik = (AA]k^) / 2, where AAj^ is the amplitude of the tone and 

baud of the received signal. In the simulation, NOISE is added to the received 
signal, y](n) = YRX, which is white Gaussian noise with zero mean and vari- 
ance = 0]2 computed in (49). A single input wideband SNR is input by the 
simulation user for all bauds, 1. The following simulation variables are used 
for the above parameters: 

Oi2 = NOSVAR(LL), 

K = KPTS(LL), 

AAik = AA(LL,K), 
kx = KX(LL), and 
SNRWBj = SNRIN. 
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IV. TESTING THE SIMULATION AGAINST MFQPSK THEORY 



A. BACKGROUND 

The purpose of developing this simulation was to create an experimental 
tool to test and analyze various Doppler, synchronization, and coding tech- 
niques and/or algorithms to support the implementation and testing of an 
MFQPSK acoustic link. Before the simulation could be used in this capacity, 
it had to be tested against MFQPSK signal theory. There was no actual test 
data to use for comparison; therefore, tests were performed to verify that the 
simulation output, the received signal, agrees with the theory of MFQPSK 
signals. It is anticipated that the simulation will be tested further when 
actual test data is available. 

It has been shown that, in theory, the output SNR (SNRqu't) equals the 
input narrowband SNR (SNRNBjjsf) for MFQPSK signals in additive white 
Gaussian noise memoryless channels that are demodulated coherently with a 
Discrete Fourier Transform (DFT) [Ref. 1]. The results of testing the simula- 
tion's ability to reproduce this result are described and presented below. 

B. THE TESTING METHODOLOGY 

The following approach was taken to verify that the simulation performs 
at an acceptable level (i.e., SNRqut = SNRNBjjsj). A set of 5 signal packets 

were generated with the simulation with a given input wideband SNR 
(SNRWBjj^). Each packet contained a single unique baud type; therefore, 

each simulation run consisted of the 5 different baud types. Referring to 
Table I, the 5 different baud types are 
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Baud Type 1: 256 sample points, 

Baud Type 2: 512 sample points, 

Baud Type 3; 1024 sample points. 

Baud Type 4: 2048 sample points, and 
Baud Type 5; 4096 sample points. 

Recall from Chapter III that if all the tones have the same power within a 
baud (i.e., the amplitudes of tones within a baud are all equal), then the input 
wideband SNR is equal to the input narrowband SNR for that baud. All the 
runs used for this analysis were produced with normalized amplitudes within 
each baud. The phases that were transmitted were generated randomly by 
the computer code. 

It is known that the minimum input SNR required to accurately decode 
MFQPSK is approximately 15 dB. This analysis considered input narrowband 
SNRs from 0 to 20 dB. The initial conditions for the transmitter's position 
and velocity were chosen from the NOSC track data. These initial conditions 
are 

XO = -5.0 feet, VXAVG = 1.0 ft/sec, VXVAR = 0.0025 (fVsec)2, 

YO = 5.0 feet, VYAVG = 5.4 ft/sec, VYVAR = 0.0025 (ft/sec)2, 

ZO = -5400.0 feet, VZAVG = 0.0 ft/sec, and VZVAR = 0.0025 (ft/sec)2. 

The initial conditions for the speed of sound was C0=4900 ft/sec and CVAR=0. 
Note the transmitter was placed almost directly above the receiver to reduce 
the amount of Doppler shift in the received signal. A description and 
examples of the simulation inputs and outputs are presented in Appendix B. 
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An average output SNR was estimated for each baud using the real and 
imaginary parts of the received signal's DFT. A kx*point DFT was performed 

on each baud of the received signal to decode the phase of the received signal. 
(The real and imaginary parts of the received signal's DFT contain the phase 
information sent on each of the tones between harmonic numbers and k£ 
in each baud.) If = ti/ 4 was transmitted (which is in the first quadrant), 
then the received signal DFT's real and imaginary parts should be located in 
the first quadrant (i.e., the real and imaginary parts should be positive). 
Likewise, for a transmitted phase angle of 3jr/4, -37i/4, or -tc/ 4 the received 
signal's DFT phase angle (i.e., real and imaginary parts) should lie in the 
second, third, or fourth quadrant, respectively. The conditional distributions 
within each baud were produced by categorizing the DFT's real and imagin- 
ary parts on a tone-by-tone basis given the phase that was transmitted on 
each tone between harmonic numbers kj and k 2 . There are two distributions 

of the received signal DFT in each of the four quadrants, a distribution of the 
real part and one of the imaginary part. The real and imaginary parts are 
statistically independent [Ref 1]. 

For example, in baud type 5 there are 4096 samples of the 256 tones en- 
coded with phase information between harmonic numbers ki=1073 and 
k2=1328. There should be approximately 64 (i.e., one fourth of 256) tones 
transmitted and received in each quadrant because the four transmitted 
phases were generated randomly from a imiform distribution. The distribu- 
tion of each DFT component should be close to a Gaussian distribution. 

Once the DFT of K tones in each baud were conditionally partitioned by 
knowing the transmitted phase (i.e., quadrant) into the four quadrants, then 
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the sample mean and variance were computed on each of 4 real part distribu- 
tions and each of the 4 imaginary part distributions. The sample means are, 
Rj and Ij, where the subscript i denotes the quadrant number. The sample 

means for each quadrant in a baud are 



N. 

1 

R. = (1/N.) ^ R.. for i = 1, 2, 3, or 4 

1 1 n 

j=l 



(50) 



and 



where 



N. 

1 

I. = (1/N.) ^ I., for i = 1, 2, 3, or 4 
1 1 ^ 11 
j=l 



(51) 



Ni = the number of times the phase corresponding to quadrant i was 
sent in that baud 

Rjj = the real part of the received signal's DFT corresponding to a 
harmonic number that was transmitted with phase in quadrant i 

Ijj = the imaginary part of the received signal's DFT correspond- 
ing to a harmonic number that was transmitted with phase in 
quadrant i. 



The sample variances of a given quadrant in a baud are 

N. 

1 

j=l 



and 




(52) 



(53) 
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The output SNRs for the quadrant within a baud were estimated by the 
sample mean squared divided by the sample variance as follows [Ref.l] 



SNROUTp. = (RJ^/4i 

R,i 1 R,i (54) 



and 



SNROUTj^ = (ip^/sj. . 



(55) 



The mean squared value of the real and imaginary parts of the DFT is a mea- 
sure of the signal power received in that quadrant. The variance of the real 
and imaginary parts is a measure of the noise received in that quadrant. 
These values form an estimate of the output SNR. Appendix C contains the 
statistics for the real and imaginary distributions of each of the four quad- 
rants and the overall statistics for each baud. 

To derive an estimate of the output SNR for each baud, the real and 
imaginary distributions of each quadrant were averaged over all four quad- 
rants. There were 2K points in the overall mean and in the overall variance 
due to the K real parts and K imaginary parts. SNROUT is the overall mean 
squared divided by the overall variance. Figure 14 is a plot of SNROUT ver- 
sus the input narrowband SNR, where SNROUT is described by 



snSut = 







+ N. 
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(56) 
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Zn-k 

i=l 

is the total number of MFM tones in each baud (see Chapter III, Section A). 

Note that the absolute value of the sample means must be used in the 
estimate (56), because, if they are not used, the overall mean will always be 
close to zero. This is due to the fact that the eight sample means are 
approximately equal in magnitude, but four have opposite signs due to their 
quadrant location. 
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C. RESULTS 

If the simulation results agree with the MFQPSK theory, then the input 
narrowband SNR should equal the output SNR. As can be seen in Figure 14, 
baud types 2, 3, and 4 agree within one dB with the MFQPSK theory. The 
longest baud, baud type 5, and the shortest baud, baud type 1, produce esti- 
mated output SNRs that are within two dB of the input SNR. As the input 
SNR increases, it can be observed from Figure 14 that the output SNRs for 
baud types 1 and 5 converge towards the value of the input SNR. 

In conclusion, the output SNR is approximately equal to the input nar- 
rowband SNR; therefore, the simulation reproduces what the MFQPSK the- 
ory predicts. 
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V. OUTPUT SNR DEGRADATION DUE TO DOPPLER 



A. BACKGROUND 

When the MFQPSK signal is transmitted from a moving platform, 
Doppler will shift the received signal's frequencies. The received signal (at 
frequency coj^', (17),) is sampled at cOy, (18), baud by baud. In general, aj, the 
time Doppler compression/expansion factor, does not equal ttm, the Doppler 
factor set for the m^l' Doppler channel; hence, there will be some Doppler mis- 
match. This mismatch will cause degradation in the output SNR. 

The output SNR degradation versus the Doppler mismatch was simu- 
lated and compared to theory. If inter-s)mibol interference (ISI), which is 
energy leaking from one tone to the next tone, is negligible, then the output 
SNR should degrade due to the Doppler mismatch as sinc(ne/4) from (27a) 
where e = oti/Aa. Derivation of the theoretical degradation and the comparison 

to the simulation results follows. 

B. ANALYSIS APPROACH 

The simulation was run several times. Each run generated one signal 
packet with an input narrowband SNR = 15 dB. (The values of 15 dB is con- 
sidered to be the minimum input SNR required to decode the MFQPSK 
signal.) A packet consisted of five bauds of baud type 3 (i.e., k^ = 1024 points 

and K = 64 tones). Baud type 3 was chosen based on the results of the 
simulation versus MFQPSK theory in the previous chapter. Recall that baud 
type 3 was in close agreement with the theory of MFQPSK signals (i.e., 
SNRqut = SNRNBjj^t). The Doppler channel, m, was set to zero for all the 
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packets. If m = 0, then a] - = aj - m Aa = aj. Thus, by increasing aj, the 

amount of Doppler mismatch increases. 

Each packet was generated with a different fixed value of aj = e Aa, 
where e was fixed at 0, 0.25, 0.5, 0.75, 1.0, ..., 2.5, 2.75 and Aa = l/(8k2). There 
were 12 packets generated, one for each value of e. The phases in the pass- 
band of each baud were selected randomly in the simulation. The initial con- 
ditions on the transmitter position and velocity and the speed of sound were 
the same as used in the earlier simulation analysis of the previous chapter. 

For each baud, the conditional statistics based on knowing the transmit- 
ted phase (i.e., quadrant) and an estimate of the output SNR (56) were com- 
puted with the method described in the previous chapter. A packet's 

A 

estimated output SNR, SNROUT(e), is simply the average of the estimated 

A 

output SNR’s from each of the five bauds within that packet. SNROUT(e) is 

SNROUT(e) = (1/5) ^ SNROUT(e). . 

(57) 

Theoretically, if ISI is neglected, the output SNR, SNRqut should de- 
grade due to Doppler, aj, from (26a) as (SNRjn sinc2(2jt k2 a])). Substituting 
a = e Aa = e / (8 k2) yields 

SNRqut = SNRjn + 20 logio [sin(7t€/4)/(7te/4)] (58) 

where SNRqut is in dB. For this analysis, the input SNR, SNRj^, was equal 

A 

to 15 dB. SNROUT(e) and SNRqut were plotted versus e to compare the 

simulation results with the theoretical degradation for an input SNR of 15 dB 
(see Figure 15). 
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Figure 15. SNRqut vs. ALPHA for SNRjn = 15 dB 



C. RESULTS AND CONCLUSIONS 

There is a significant difference between the simulation results and 
theory as shown in Figure 15. It was speculated that this difference is due to 
ISI. As a\ increases, the tones ’smear” into adjacent tones (i.e., ISI). To prove 

whether or not the difference was due to ISI, the analysis was repeated Avith 
an input SNR of 5 db. By lowering the input SNR to 5 dB, the signal level is 
not much higher than the noise level within the band; therefore, symbol inter- 
ference from adjacent tones is at about the same level as the noise level of 
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that tone. The noise level of each tone, at = 5 dB, did not increase as 

much relative to increases in aj as when SNRjn = 15 dB; thus, the output 
SNR at a given tone should not degrade as rapidly as aj increases. Figure 16 
is a plot of the theoretical degradation and the simulation results for both 
input SNR's (5 and 15 dB). By examining Figure 16 one can see that, indeed, 
the output SNR at 5 dB degrades at a much slower rate than the output SNR 
at 15 dB. The statistics of each packet for both 5 and 15 dB input SNR on a 
baud by baud basis are presented in Appendix D. 

The analysis of the simulation results for 15 and 5 dB seemed to sub- 
stantiate that ISI can not be neglected. An estimate and correction of Doppler 
shift may be used to reduce the ISI and lessen the output SNR degradation at 
higher SNRjn's. An estimate of the Doppler shift can be fed back into the re- 
ceiver sampling system to compensate for the Doppler shift. 
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Figure 16. SNRqut vs. ALPHA for SNRi.nj = 15 dB and 5 dB 
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VI. A DOPPLER ESTIMATE FOR FINE CORRECTIONS 



A. BACKGROUND 

The Doppler shifts due to the moving transmitter should be estimated so 
that they can be removed from the received signal in the receiver system to 
reduce the errors in decoding the MFQPSK signal. One technique for estimat- 
ing the Doppler compression/expansion factor was examined. An adaptation 
of a Doppler estimation method recently proposed for ordinary QPSK [Ref 5] 
was developed for use with MFQPSK and is the subject of this chapter. 

B. IMPLEMENTATION OF THE DOPPLER ESTIMATION 

This Doppler estimation method approximates the Doppler time 
compression/expansion factor within a single baud. To implement this esti- 
mator, several simulation runs were produced. A simulation rim consisted of 
a single packet. In each packet, five bauds of baud type 3 (i.e., k^ = 1024 

points and K = 64 tones between harmonics 269 and 332) were generated 
with randomly selected phases. The initial conditions of the transmitter's 
position and velocity are the same as used in the simulation versus MFQPSK 
theory in Chapter IV. Each packet was generated with m, the Doppler 
channel number, equal to zero (i.e., = 0, which means the residual Doppler 

mismatch, a i - otj^, is equal to aj) and a fixed a = aj = e Aa where e is between 
0 and 1.5 and Aa = l/( 8 k 2 ). For this analysis (using baud type 3), k 2 is 332; 
therefore, Aa = 1/2656. Recall a = tt] is the Doppler compression/expansion 
factor. A set of runs or packets were produced for four different input 
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narrowband SNR’s (10, 15, 20, and 40 dB). The amplitudes for all the packets 
were normalized (i.e., AA(LL,K) = 1 for all LL and K), 

Comparison of the phase for a given tone in the first half of a baud to the 
phase of the same tone in the second half of the baud yields an estimate of the 
Doppler compression/expansion factor within that baud. If there is no Doppler 
shift and no noise on the signal, the corresponding phases from the first and 
last half of the baud should be equal. The following method was used to 
compute aj which is an estimate of the Doppler compression/expansion factor 

within the baud. To obtain the phases of the first and second half of the 
received signal’s baud, a 512 (i.e., kx/2) point DFT was taken on the first 512 

points of the baud and another 512 point DFT was taken on the last 512 
points of the baud. Note that, when the bauds were generated with the simu- 
lation, the odd tones within the passband (harmonics between 269 and 332) 
were zeroed out by setting their amplitudes to zero (all the harmonics outside 
the passband are always zero). If the odd harmonics in the passband are not 
zero, then the output of the DFT's of the first and last halves of the baud will 
have interference between harmonic numbers because 1/2 of an odd harmonic 
does not produce an integer harmonic, causing an extra half cycle of signal 
between the harmonics, k', of the 512 point DFT's. With an ideal received 
signal (i.e., no noise and no Doppler), if the odd harmonics in the passband 
are nonzero, the corresponding phases of the harmonics of the first and last 
half of the baud will not be equal due to the interference of the signals 
between each harmonic in the 512 point DFT’s. 

A Doppler estimate, aj, was then computed for each baud by averaging 

the differences between the received signal’s phases from the first half and 



56 



the received signal's phases from the second half of the baud. Thus, the 
Doppler compression/expansion factor estimate for a baud is 



a = [l/(K/2)] X [PHASE of (S2*(k’) Sl(k'))] / (2rck’) 
k' 

where, 

S2*(k') = Conjugate of the DFT of the second half of the baud 
Sl(k') = The DFT of the first half of the baud 

k’ = A harmonic number of the first and second half DFT's 
(Note k’ = k/2, where k are the harmonics of the original 
kx point DFT of the baud.) 

K = The number of MFM tones in the original baud. 



(59) 



An average Doppler estimate, a, was calculated for each packet by averaging 
the five Doppler estimates, 6t], within each packet and was plotted versus a 

for each of the four input SNRs (see Figure 17). Appendix E contains the 
estimates and statistics for each packet at each input SNR. 



C. RESULTS AND CONCLUSIONS 

From the analysis for a = tt] between 0 and 1.5 Aa, a does not estimate a 

accurately for input SNRs less than 15 dB. Recall that an input SNR = 15 db 
is the minimum for successfully decoding this MFQPSK signal. At an input 
SNR of 15 dB, a estimates a well for a between 0 and 0.8 Aa; and for input 
SNRs of 20 dB or greater, a estimates a for a between 0 and Aa. These 
results are shown clearly in Figure 17. 

The estimate a is not acceptable for large Doppler shifts (i.e., Doppler 
shifts > Aa); therefore, this Doppler estimation/correction method may only be 
used for fine corrections of Doppler within a received signal baud and Doppler 
channel, where the maximum Doppler mismatch is ± Aa/2. Note from 
Figure 17 that the estimation is perfect for 0.5 Aa. 
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Figure 17. ALPHA-HAT vs. ALPHA for SNRin = 10, 15, 20, and 40 dB 
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VII. CONCLUSIONS AND RECOMMENDATIONS 



A numerical code to simulate a MFQPSK signal received from a moving 
transmitter through a bandpass channel has been developed. It resides 
presently on the NPS IBM mainframe computer. The simulation is written in 
FORTRAN 77 and the code is included in Appendix F. 

The model of the transmitter platform dynamics in the simulation is hm- 
ited to straight line motion with random fluctuations. Actual transmitter 
platforms have motion in the pitch, yaw, and roll planes too and thus a six- 
degree-of-freedom model may be more realistic for future studies of platform 
motion effects on MFQPSK signalling. 

The simulation should be tested further with actual test data to ensure 
that it produces a realistic received signal. Validation with actual test data 
may show that other channel characteristics or parameters which have not 
been considered in the simulation are important and must be included. 

An important discovery of the simulation was that the output SNR 
degradation due to Doppler shift was primarily due to ISI. ISI cannot be 
neglected; therefore a theoretical analysis of ISI should be developed to real- 
istically analyze MFQPSK in the presence of Doppler. 

The Doppler estimation/correction method analyzed and presented in 
chapter VI only appears useful to remove small Doppler shifts (i.e., for fine 
adjustments only). Analysis of methods to remove larger Doppler shifts (i.e., 
coarse adjustments) across several bauds or a packet should be performed for 
MFQPSK signals. 
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Synchronization error estimation and correction has not yet been 
addressed with this simulation, but should be in the future. 

This NFS developed MFQPSK signal has been implemented in hardware 
by T. Gantenbein and Dr. P. H. Moose using differential phase coding [Ref. 6]. 
It is conjectured that with differential coding the Doppler shifts, ISI, and syn- 
chronization errors will not affect the received signal as significantly as they 
affect the individually phase coded signal which is the coding technique 
presently implemented in the simulation. Therefore, the differential phase 
coding should be coded in the simulation. Doppler, ISI, and synchronization 
estimation/correction methods shoiild be studied that are compatible with dif- 
ferential phase coding. 
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APPENDIX A 



DERIVATION OF £(u) 



Given the model of the channel's time delays below in Figure 18, assume 
x(u) is a superposition of N impulses; 



N 

x(u) = ^ A. 5 (u-Uj) . 
i=l 



(60) 



x(u) 



FIXED DELAY: 




TIME VARYING 
DELAY: 






X| -h ttjU 



Figure 18. Time Delay Model of Channel 

An impulse sent at time u = 0 arrives at y at time There- 

fore, the output signal from the channel y(Tj^ + T; + ccjli^) is the impulse x(0); 
that is, y(Tjj -i- Tj -t- cxiT]^) = x(0). An impulse sent at time AT arrives at the 
output of the channel at time 

AT -I- Tk -t- X] -h ai(AT + Tk). (61) 

Therefore, the output signal from the channel 

y(AT + tk + ti -I- a](AT -i- tk)) (62) 
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which is the impulse x(AT); is 



y(AT + Tk + XI + aj(AT + Xk)) = x(AT) 



Now letting 



u = AT + Xk + X] + ai(AT + Xk) 



which can be rewritten as 



u = AT(1 + aj) + Xk(l + aj) + xi , 



then solving for AT yields 



Thus, 



AT = ( (u - Tj) / (1 + aj) ) - Tk = u - Tk - ((xi + aju) / (1 + a])). 



y(u) = x( AT ) = x( u - Tk - ((■^1 + oqu) / (1 + «])) ) 



or 



y(u) = x( u - Xk - £(u)) 



where 



£(u) = ((Xj + ttju) / (1 + tti)). 



(63) 

(64) 

(65) 

( 66 ) 

(67) 

( 68 ) 
(69) 
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APPENDIX B. AN EXAMPLE OF A SIMULATION RUN 



A. INPUT FROM THE SCREEN 



PLEASE ENTER THE INITIAL POSITION X0,Y0,Z0 (FT) 
OF THE RECEIVER RELATIVE TO THE TRANSMITTER . . . 
XO = -10.0000 YO = 20.0000 ZO = -6000.000 



ENTER THE TRANSMITTER"S VELOCITY AVERAGE AND VARIANCE 
IN THE X, Y, Z-DIRECTIONS (FT/SEC) AND (FT/SEO^XZ ... 
(I.E. VXAVG, VXVAR, VYAVG, VYVAR, VZAVG, VZVAR) 

VXAVG = 2.00000000 VXVAR = 0 . 500 0 00 0 07 E-0 1 

VYAVG = A. 00000000 VYVAR = 0 . 499999896 E-02 

VZAVG = 0 . OOOOOOOOOE+00 VZVAR = 0 . 249999994E-02 



ENTER THE AVERAGE SPEED OF SOUND IN FT/SEC 
AND THE VARIANCE IN (FT/SEC)5(X2 ... 

CO = 5000.00000 CVAR = 0 . OOOOOOOOOE+00 



ENTER THE TRANSMITTER"S DOWN LINK TRANSMITTED 
HALF BEAM WIDTH ANGLE (DEG) ... 

THETAO = 4.00000000 



DO YOU WANT NORMALIZED AMPLITUDE FOR THE RECEIVED SIGNAL? 

PLEASE ENTER IsYES OR 0:NO ... 
lAMP = 1 

ENTER THE # OF PACKETS IN THE TRANSMITTED SIGNAL ... 

NPAKS = 2 
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THIS PROGRAM ENCODES A QPSK MULT I FREQUENCY SIGNAL. 
THE PHASES ARE SHONN BELON FOR ONE FREQUENCY . . . 



2 



1 



3 



A 



SELECT ONE OF THE FOLLONING METHODS FOR ENCODING 
THE PHASES FOR EACH FREQUENCY FOR EVERY BAUD 
MITHIN EACH PACKET ... 

ENTER . . . 

1: THE PROGRAM RANDOMLY SELECTS ALL THE PHASES 
FOR ALL 2 PACKETS 
2: YOU INDIVIDUALLY SELECT THE PHASES 
METHOD = 1 



ENTER THE BAUD TYPE ^ CORRESPOtJDI NG TO THE FOLLOWING 
BAUD LENGTH FOR PACKET NUMBER: 1 ... 



1 : 


BAUD 


LENGTH 


(DELT) = 1/2^0 


SECONDS 


2 : 


BAUD 


LENGTH 


(DELT) = 1/120 


SECOf^DS 


3 : 


BAUD 


LENGTH 


(DELT) = 1/60 


SECONDS 


^ : 


BAUD 


LENGTH 


(DELT) = 1/30 


SECONDS 


5 : 

BDTYPE(l) = 


BAUD 


LENGTH 

1 


(DELT) = 1/15 


SECONDS 



ENTER THE f^UMBER OF BAUDS IN PACKET NUMBER: 1 ... 

NBAUDS(l) = 1 



ENTER THE BAUD TYPE # CORRESPONDING TO THE FOLLOWING 
BAUD LENGTH FOR PACKET NUMBER: 2 ... 



1 

2 

3 

A 

5 

BDTYPEC2) = 



BAUD LENGTH (DELT) 
BAUD LENGTH (DELT) 
BAUD LEfJGTH (DELT) 
BAUD LENGTH (DELT) 
BAUD LENGTH (DELT) 
2 



1/2^0 SECONDS 
1/120 SECONDS 
1/60 SECONDS 
1/50 SECONDS 
1/15 SECONDS 



ENTER THE NUMBER OF BAUDS IN PACKET NUMBER: 2 ... 

NBAUDS(2) = 1 
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WOULD YOU LIKE THE DISCETE FOURIER TRANSFORM 

OF THE OUTPUT SIGNAL ? (ENTER I^YES OR 0:NO) ... 
IDFT = I 



WOULD YOU LIKE THE DFT OUTPUT WINDOWED ? 

ENTER 1:YES OR 0:NO ... 

IWNDOW = 1 



PLEASE ENTER THE DESIRED INPUT WIDE BAND 
SIGNAL-TO-HOISE RATIO IN DB .... 

SNRDB = 15.0000 
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B. LIST OUTPUT TO DISK FILE 30 
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LL 


K 


PHI(LL,K) 


IPHI 


1 


68 


-0.785 


9 


1 


69 


0.785 


1 


1 


•70 


0.785 


1 


1 


71 


2.556 


2 


1 


72 


-2.356 


3 


1 


75 


-2.356 


3 


1 


7^ 


-2.556 


5 


1 


75 


-0.785 


9 


1 


76 


2.356 


2 


1 


77 


-0.785 


9 


1 


78 


2.356 


2 


1 


79 


0.785 


1 


1 


80 


-0.785 


9 


1 


81 


0.785 


1 


1 


82 


0.785 


1 


1 


83 


-2.356 


3 


2 


135 


-0.785 


9 


2 


136 


-2.356 


3 


2 


137 


-2.356 


3 


2 


138 


2.356 


2 


2 


139 


-2.356 


3 


2 


190 


2.356 


2 


2 


191 


-0.785 


9 


2 


192 


-0.785 


9 


2 


193 


0.735 


1 


2 


199 


-0.785 


9 


2 


195 


2.356 


2 


2 


196 


-0.785 


9 


2 


197 


2 . 556 


2 


2 


198 


0 .785 


1 


2 


199 


0.785 


1 


2 


150 


0.785 


1 


2 


151 


-2.356 


3 


2 


152 


-0.785 


9 


2 


153 


-0.785 


9 


2 


159 


2 . 356 


2 


2 


155 


-2.356 


3 


2 


156 


2.356 


2 


2 


157 


2,356 


2 


2 


158 


-2.356 


3 


2 


159 


2.356 


2 


2 


160 


2.356 


2 


2 


161 


, 0.785 


1 


2 


162 


2.356 


2 


2 


163 


2.356 


2 


2 


169 


-2.356 


3 


2 


165 


-2.356 


3 


2 


166 


2.356 


2 
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u 

0 

1 

2 

3 

5 

6 

7 

8 

9 

10 

11 

12 

13 

19 

15 

16 

17 

18 

19 

20 

21 

22 

23 

29 

25 

26 

27 

28 

29 

30 

31 

32 

33 

39 

35 

36 

37 

38 

39 

90 

91 

92 

93 

99 

95 

96 

97 

98 

99 

50 

51 

52 

53 

59 

55 

56 

57 

58 

59 

60 

61 

62 

65 

69 



INPUT DATA FOR 
REAL PART 
0.200989E+01 
-.291531E+01 
122278E+01 
0 . 225909E + 01 
0 .176306E+01 
1 18927E+01 
-.213088E-01 
-.255857E+00 
“.615879E+00 
0.199238E+01 
-.130629E+01 
-.392880E+00 
0 . 207660E+01 
-.329395E+00 
“.369812E+01 
0.398756E+01 
0 . 12591 5E+01 
-.985057E+01 
0.209897E+01 
0 .725176E+01 
-.267107E+01 
“.327209E+01 
0.5026S7E+01 
-.287833E+01 
-.696552E+01 
0.376308E+01 
0 .269107E+01 
-.925671E+01 
192092E+01 
0.917586E+01 

- .685801E+00 
-.910082E-01 
0.2571 53E^-01 
0 .686129E+00 
“ . 128562E+01 
0 . 271385E+01 
0 . 193209E+01 
0 . 122657E+01 
0 . 156959E+01 
-.215376E+01 
“.159369E+00 
0 .906599E+01 
0 . 1 15762E+01 
-.275395E+01 
0.252170E+01 
“.110620E+01 
-.378736E+00 
0 . 253587E+01 
-.999909E+00 

336829E+01 
0 . 299053E+01 
-.839522E+00 
159595E+01 
“.102779E+01 
0 . 128035E+01 
565980E+00 
-.550359E401 
-.231S39E400 
0 . 125277E+01 
-.697582E+00 
0 .705833E+00 
0. 133619E401 
0.229186E400 

- . 556922E+00 
0 .999179E+00 



BAUD H 1 
IMAG PART 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+OO 
0 . OOOOOOE+OO 
0 . OOOOOOE+-00 
0 . OOOOOOE+OO 
0 . OOOOOOE+OO 
0 .000000E+-00 
0 .OOOOOOE+OO 
0. OOOOOOE + OO* 
0 .OOOOOOE+OO 
0 . OOOOOOE + OO 
0 .OOOOOOE+OO 
0 . OOOOOOE+ 00 
0. OOOOOOE+OO 
0 .OOOOOOE+OO 
0 . OOOOOOE+OO 
0 . OOOOOOE+OO 
0. OOOOOOE+OO 
0 .OOOOOOE+OO 
0 .OOOOOOE + OO 
0 . OOOOOOE+OO 
0 .OOOOOOE+OO 
0 .OOOOOOE+OO 
0 .OOOOOOE + OO 
0 .OOOOOOE+OO 
0 .OOOOOOE+OO 
0. OOOOOOE+OO 
0. OOOOOOE+OO 
0 .OOOOOOE+OO 
0 . OOOOOOE+OO 
0. OOOOOOE+OO 
0. OOOOOOE+OO 
0 .OOOOOOE+OO 
0 . OOOOOOE+OO 
0 .OOOOOOE+OO 
0 . OOOOOOE+ 00 
0 .OOOOOOE+OO 
0 .OOOOOOE+OO 
0 .OOOOOOE+OO 
0 . OOOOOOE+OO 
0 .OOOOOOE+OO 
0 . OOOOOOE + OO 
0 . OOOOOOE + OO 
0 . OOOOOOE+OO 
0 . OOOOOOE+OO 
0 .OOOOOOE+OO 
0 .OOOOOOE+OO 
0 . OOOOOOE+ 00 
0 . OOOOOOE+OO 
0 . OOOOOOE+OO 
0 .OOOOOOE+OO 
0 . OOOOOOE+OO 
0 . OOOOOOE+OO 
0 . OOOOOOE+OO 
0 .OOOOOOE+OO 
0 .OOOOOOE+OO 
0 .OOOOOOE+OO 
0 . OOOOOOE+OO 
0 . 000000E+ 00 
0 . OOOOOOE+OO 
0 . OOOOOOE+OO 
0 .OOOOOOE+OO 
0 .OOOOOOE+OO 
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65 


0.^13^70E+01 


0 . OOOOOOE+00 


66 


175258E+01 


0 . OOOOOOE+00 


67 


-.939^83E+00 


0 . OOOOOOE+00 


68 


0 .279167E+00 


0 . OOOOOOE+00 


69 


-.21137^E+00 


0 . OOOOOOE + 00 


70 


0 .7^9620E+00 


0 . OOOOOOE+00 


71 


0 .227079E + 01 


0 . 000000E^^00 


72 


0 . 16996^E+01 


0 . OOOOOOE + 00 


73 


360312E+01 


0 . 000000E4 00 


7^ 


0 . 181599E+01 


0 . OOOOOOE + 00 


75 


0 . 192265E+01 


0 . OOOOOOE+00 


76 


“.293610E+01 


0 . OOO00OE-^OO 


77 


- . 571738E+00 


0 . OOOOOOE+00 


78 


0.171358E+01 


0 . OOOOOOE+00 


79 


-.227510E+01 


0 . OOOOOOE+00 


80 


0 . 15^18^E+01 


0 . OOOOOOE+00 


81 


0.2^0^58E+01 


0 . OOOOOOE+00 


82 


192191E+01 


0 . OOOOOOE+00 


83 


-.17113^E+01 


0 . OOOOOOE^-00 • 


8^ 


0.3^9656E+01 


0 . OOOOOOE^ 00 


85 


>.21^113E+01 


0 . OOOOOOE + 00 


86 


0 .219070E+01 


0 . OOOOOOE+00 


87 


0.365189E401 


0 . OOOOOOE+00 


88 


0.713565E+00 


0 . OOOOOOE+00 


89 


-.2^0271E+01 


O.OOOOOOE+00 


90 


0 .625637E+00 


0 . OOOOOOE+00 


91 


0 . 177^16E + 01 


0 . OOOOOOE+00 


92 


“.563617E*01 


0 . OOOOOOEfOO 


93 


0 . 517219E^01 


0 . OOOOOOE+ 00 


9^ 


0.379215E+01 


0 . OOOOOOE+OO 


95 


-.5^078^E+01 


0 . OOOOOOE400 


96 


0 . ^280^7E+01 


0. OOOOOOE+00 


97 


0.539627E+01 


0 . OOOOOOE+00 


98 


-.^33^80E+01 


0 . OOOOOOE+00 


99 


“.32033^^E+00 


0 . OOOOOOE+00 


100 


0 .50^‘^50E■^01 


O.OOOOOOE-fOO 


101 


-.399863E+01 


0 . OOOOOOE+OO 


102 


180^50E+01 


0 . OOOOOOE+OO 


103 


0 . 731280E+01 


0. OOOOOOE+OO 


10^ 


-.358992E+01 


0. OOOOOOE+OO 


1 05 


“.6866^1E+01 


0 .OOOOOOE+OO 


106 


0.6782^i5E+01 


0 .OOOOOOE + OO 


107 


0.353709E'»01 


0 .OOOOOOE + OO 


108 


“.738260E»01 


0 .OOOOOOE + OO 


109 


0 . 855238E^ 00 


0 .OOOOOOE+OO 


no 


0 . 565053E+ 01 


0 . OOOOOOE+OO 


111 


^66336E^01 


0 .OOOOOOE+OO 


112 


“.2026^6E+01 


0 . OOOOOOE+OO 


113 


0.^il6927E>01 


0 .OOOOOOE+OO 


11^ 


0.963901E^00 


0. OOOOOOE+OO 


115 


129393E+01 


0 .OOOOOOE+OO 


116 


- . 182671E+00 


0 . OOOOOOE+OO 


117 


0 .31^196E+01 


0 . OOOOOOE+OO 


118 


-.6^1^27E+00 


0. OOOOOOE + OO 


119 


-.229205E+01 


0 .OOOOOOE+OO 


120 


0 .89^^7 37 E-02 


0 .OOOOOOE+OO 


121 


“.208729E+01 


0 .OOOOOOE + OO 


122 


0 . 127^l76E^01 


0 .OOOOOOE + OO 


123 


0 . 189270E* 01 


0. OOOOOOE+OO 




0.827031E^OO 


0 . OOOOOOE + OO 


125 


180^2^E+01 


0 .OOOOOOE+OO 


126 


0 . 181878E-01 


0 . OOOOOOE+OO 


127 


0 .7 38333EfOO 


0 . OOOOOOE + OO 


128 


0 . 1 30366E» 01 


0 . OOOOOOE + OO 


129 


0 . 626686E+00 


0 . OOOOOOE+OO 


130 


0 . 20001 lE+Ol 


0. OOOOOOE+OO 


131 


“.176676E+01 


0 . OOOOOOE+OO 


132 


- . 18^777 E-^01 


0. OOOOOOE+OO 


133 


0.305751E»01 


0 . OOOOOOE+OO 


13^ 


-.^4l575^Ef 00 


0 .OOOOOOE+OO 


135 


-.35667^E+01 


0 . OOOOOOE+OO 


1 36 


0.^29929E+01 


0 . OOOOOOE+OO 
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137 

138 

139 

1^0 

191 

192 

195 

199 

195 

196 

197 

198 

199 

150 

151 

152 

153 

159 

155 

156 

157 

158 

159 

160 

161 

162 

165 

169 

165 

166 

167 

168 

169 

170 

171 

172 

173 

179 

175 

176 

177 

178 

179 

180 

181 

182 

183 

189 

185 

186 

187 

188 

189 

190 

191 

192 

193 

199 

195 

196 

197 

198 

199 

200 

201 

202 

203 

209 

205 

206 

207 

208 



0 .33197 0E + 01 
-.989755E+01 

127909E+01 
0 .602902E+01 
0 .897187E-01 
-.553555E+01 
0 . 508767E+01 
0 . 591502E+01 
-.296052E+01 
0 .850205E+00 
0.319615E+01 
- .907688E+01 
-.112058E+00 
0 .956206E+01 
0.205910E+01 
-.292958E+01 
0.360312E+01 
0.238991E+01 
0.885219E+00 
-,507907E+01 
0 .230592E+01 
0 . 173215E+01 
-,932735E+00 
0 .965825E+00 
0 . 193239E+01 
0 .579S99E-01 
- . 168030E+01 
0 . 350895E+01 

- .280785E+01 

376679E+01 
0.258339E+01 
0 .592977E+00 
“.919950E+01 
0 . 111968E+01 
0 .27 1280E+01 
-.502696E+01 
-.321359E+01 
0 .732952E+01 
“.217958E+01 
-.501155E+01 
0.501020E+01 
0 . 179601E+01 
“ . 189239E+01 
0 . 396009E^ 01 
-.109039E+01 
-.591859E+00 
0 .917718E + 00 
-.272576E+01 

153335E+01 
0 . 59091 OE+01 

- . 250186E+00 
-.513901E+01 
0 . 968525E+01 
0 . 31 0399E+01 

- .999965E+01 
“.127512E+01 
0 .250906EJ 01 

- .666329E+00 

1 39897E+01 
0 . 570889E+01 

- .896835E-01 
-,293930EI-01 
-.259588E»01 
-.512660E»00 
0.871717E+00 
“.115026E401 

- . 399199E401 
0 .286601 E + 01 
0 .951813E+01 
-.516909E^01 
-.216182E+01 
0 .682289E+01 



0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE + 00 
0 . OOOOOOE+00 
0 . OOOOOOE + 00 
0 . OOOOOOE^OO 
0 . OOOOOOE + 00 
0 . OOOOOOE + 00 
0 . OOOOOOE+00 
0 . OOOOOOE + 00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . O0OOOOE4OO 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
O.OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+OO 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 ,000000E + 00 
0 . OOOOOOE+00 
O.OOOOOOE+00 
0 . OOOOOOE+00 
O.OOOOOOE+00 
O.OOOOOOE^OO 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . 000000E-» 00 
0 , OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+OO 
0 . OOOOOOE+00 
0 . 000000E+ 00 
0 . OOOOOOEf 00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
O.OOOOOOEfOO 
0 . OOOOOOE+00 
0 . OOOOOOE HOO 
0 . OOOOOOE+00 
0 . OOOOOOE^-00 
0 . 000000E400 
0. 000000 E^OO 
0 . 000000E»00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . 000000E4 00 
0 . OOOOOOE^ 00 
0 . OOOOOOE + 00 
0 . OOOOOOE+00 
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209 0.316336E+00 

210 ~.910015E^01 

211 0.373189E+01 

212 0.^59230E+01 

213 -.709038E+01 

21^ -.790167E-01 

215 0 .640<^98E + 01 

216 0.885365E+00 

217 - . 10191^E+02 

218 0 . 376‘^23E+01 

219 -.109551E+00 

220 -.386055E+01 

221 -.3807^9E+01 

222 0.589190E+01 

223 0 .23<^^06E+00 

22^ -.^89929E+01 

225 “.257231E+01 

226 0.668137E+01 

227 0.256053E+00 

223 -.820556E^01 

229 0.317285E+01 

230 0.728581E+01 

231 -.353177E+01 

232 ~.^73028E+01 

233 0.651311E+01 

23^ -.125509E+01 

235 -.925538E+01 

236 0.317662E+01 

237 0.596795E+01 

238 -.^15065E+01 

239 0.927513E+00 

2^i0 0.638612E+01 

2^41 0.1 16869E + 01 

2^i2 -.316552E + 01 

2^3 0.217262E+01 

2^^ 0.15^»U2E + 01 

2^5 -.508769E+01 

2^6 0.2199^3E+01 

2^7 0.279915E+01 

2^48 -.269927E^01 

2^9 0.25^^80E+00 

250 0.117177E+01 

251 -.U5215E+01 

252 -.553^69E+00 

253 0 . 37 09^»7E + 01 

25^ 0 . 1^4689<4E+01 

255 0.176155E>01 



O . OOOOOOE ^00 
0 . OOOOOOE + 00 
0 . OOOOOOE^ 00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . 000000E+ 00 
0 . OOOOOOE + 00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE + 00 
0 . OOOOOOE + 00 
0 . OOOOOOE + 00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE + 00 
0 . OOOOOOE+00 
0 . OOOOOOE^ 00 
0 . 000000E+ 00 
0 . OOOOOOE+00 
0 . OOOOOOE + 00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . 000000E+ 00 
0 . OOOOOOE+00 
0 . 000000E+ 00 
0 . 000000E+ 00 
O.OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+ 00 
0 . OOOOOOE + 00 
0 . OOOOOOE't-00 
0 . OOOOOOE + 00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE + 00 
0 . OOOOOOE^OO 
0 . OOOOOOE+ 00 



DFT OUTPUT DATA FOR 
K REAL PART 

68 0 .906573A^iE+02 

69 0 . 12512337E + 03 

70 0 .85056396E+02 

71 -.95022675E+02 

72 -.80296AA8E+02 

73 11772A76E+03 

7A 10582A92E+03 

75 0 .9507A265E+02 

76 -.7 39953A6E-»02 

77 0 .95286A69E+02 

78 -.7 55723A2E-^02 

79 0.8A2A6826E+02 

80 0.1 1 1 2A681E-»03 

81 0 . 10516679E + 03 

82 0 . 1 1A861 5AE+03 

83 -.89A3A708E+02 



BAUD S 1 

I MAG PART 

-.900A1351E+02 
0 .87728928E + 02 
0 . 1 0223833E 4 03 
0.10510611E+03 
-.10A72539E403 
1 0389015E + 03 
-.87098938E+02 
-.n 210881 E4()3 
0.71338562E402 
-.91975923E402 
0 .88835632E+02 
0.97707855E+02 
-.82108078E+02 
0.83761368E+02 
0.10973633E403 
-.95062195E+02 



MAGNITUDE 

0 . 1277 5930E4 03 
0 . 152818AAE4 03 
0 . 1 32993A7E4 03 
0 . 1 A169193E + 03 
0 . 13196561E403 
0.157010AAL403 
0. 1 3705887E + 03 
0 . 1 A699A89E403 
0. 1027837AE403 
0 . 132A3378E»03 
0.11663167E+03 
0 . 12901299E+03 
0.13826636E+03 
0 . 13AAA7 07 E+03 
0 . 15885600E403 
0 . 13051968E403 



PHASE 

(DEG) 

AA81 1005E402 - 

0. 3503A729E + 02 0 

0 . 502A] A70E4 02 0 

0 . 1321 1565E403 0 

-.127A7855E+03 - 

1385721AE+03 - 

1A05A595E+05 - 

-.A9700256E402 - 

0. 1 360A727E403 0 

A5986572E+02 - 

0 . 1 3038779E + 03 0 

0.A9231018E+02 0 

-.36A29962E+02 - 

0.38535995E402 0 

0 . A369276AE402 0 

-.1 33252,96E«03 - 



NOTE: The DFT is windowed (lUT^DOW = 1) 



(RAD) 

78209955E+00 
.611A71A1E+00 
.87687898E+00 
.23058529E+01 
.222A9212E+01 
.2A185A00E+01 
.2A529552E+01 
.867A3307E400 
. 237AA726E + 01 
. 76 77 1 0A5E400 

22756953E+01 

.8592A359E+00 

.63582295E+00 

.6725799AE+00 

.76258266E+00 

.23257027E+01 
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DFT 


INPUT DATA FOR 


BAUD « 2 


U 


REAL PART 


IMAG PART 


0 


-.526725E+01 


0 . 000000E400 


1 


-.399182E+00 


O.OOOOOOE+00 


2 


0 .696870E+01 


0 . OOOOOOE+00 


3 


-.290971E+01 


0,000OO0E+00 




-.359591E+01 


0 . OOOOOOE+00 


5 


0 . 968083E+01 


0 . 000000E400 


6 


-.939102E+01 


0 . OOOOOOE+00 


7 


115969E401 


0 . OOOOOOE+00 


8 


0 .533558E + 01 


0 . OOOOOOE+00 


9 


-.290291E+01 


0 . OOOOOOE+00 


10 


-.667569E+01 


0 . OOOOOOE+00 


11 


0 .728567E+01 


0 . OOOOOOE+00 


12 


0 .979552E+00 


0 . OOOOOOE + 00 


13 


-.380783E+01 


0 . OOOOOOE+OO 


1^ 


0.311805E+01 


0 . OOOOOOE+OO 


15 


0.950093E+01 


0 . OOOOOOE+OO 


16 


-.996919E+01 


0 . OOOOOOE+OO 


17 


0.199608E+01 


0. OOOOOOE+OO 


18 


0.155680E+01 


0. OOOOOOE+OO 


19 


-.970798E+01 


0 .OOOOOOE+OO 


20 


0.990136E+01 


0 . OOOOOOE + OO 


21 


0 .601283E+01 


0 .OOOOOOE+OO 


22 


-.750963E+01 


0 .OOOOOOE+OO 


23 


-.399910E+01 


0 .OOOOOOE+OO 


2^4 


0 . 110179E+02 


0 . OOOOOOE+OO 


25 


0 .299212E+01 


0 .OOOOOOE+OO 


26 


-.511219E>01 


0 .OOOOOOE + OO 


27 


0.563213E+01 


0. OOOOOOE+OO 


25 


0.576661E401 


0 .OOOOOOE+OO 


29 


-.539395E+01 


0 .OOOOOOE+OO 


30 


330958E+01 


0 . OOOOOOE+OO 


31 


0 . 113159E+02 


0 . OOOOOOE+OO 


32 


0 . 199053E+ 01 


0 .OOOOOOE+OO 


33 


- . 199059E+01 


0 . OOOOOOE + OO 


3^ 


0 .998989E+01 


0 . OOOOOOE+OO 


35 


0.992170E+01 


0 .OOOOOOE+OO 


36 


-.811950E+01 


0 .OOOOOOE+OO 


37 


0 . 138805E+01 


0 . OOOOOOE+OO 


38 


0.912952E+01 


0 .OOOOOOE+OO 


39 


- . 3867 19E + 01 


0 .OOOOOOE+OO 


^0 


“.538585E+01 


0 .OOOOOOE + OO 


^1 


0.650927E+01 


0 .OOOOOOE+OO 


^42 


0 . 396169E+00 


0 .OOOOOOE+OO 


^3 


0 . 101979E+00 


0 . OOOOOOE+OO 


^49 


0 .922129E+01 


0 .OOOOOOE + OO 


^5 


-.261779E+01 


0 .OOOOOOE + OO 


^46 


-.652925E+01 


0 . OOOOOOE+OO 


^7 


0.692599E+01 


0 .OOOOOOE+OO 


^5 


0 .912038E+00 


0 .OOOOOOE+OO 


^9 


-.929529E^01 


0 . OOOOOOE+OO 


50 


0.202519E+01 


0 . OOOOOOE+OO 


51 


0 . 106559E402 


0 .OOOOOOE+OO 


52 


-.710919E+01 


0 .OOOOOOE+OO 


53 


- .957 1 07E+01 


0 .OOOOOOE+OO 


59 


0.109375E402 


0 . OOOOOOE+00 


55 


-.109101E401 


0 .OOOOOOE + OO 


56 


- .7 1 26 30Ef01 


0. OOOOOOE + OO 


57 


0 .622093E^01 


0 .OOOOOOE + OO 


55 


0 . 965999E+01 


0 .OOOOOOE+OO 
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0 . OOOOOOE+00 
0 . 000000E4 00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
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0.5735a^iE + 01 
^27 -.8^52‘^9E + 01 

^28 0.516992E>00 

<»29 0 . 1193^i3E+02 

<^30 -,<^<<2223E+01 

<431 -.389617E+01 

<432 0.306069E+01 

<433 -.7<4l859E+00 

<i3<4 -.601171E+01 

<435 0.91521 0E+00 

<436 “,6628<40E+00 

<437 -.279501E+01 

<438 0.1 02183E+01 

<i39 0.129399E+01 

<160 0.962253E+00 

661 -.105563E+00 

662 0.622626E+01 

663 -.606331E+00 

666 -.160753E+01 

665 0.689188E+00 

666 -.570710E+01 

667 -.202896E+01 

668 0.676189E+01 

669 -.360101E+01 

650 -.632951E+01 

651 0.869652E+01 

652 0.619665E+01 

653 -.520679E^01 

656 0.318501E+01 

655 0.661761E+01 

656 -,293308E+01 

657 -.569965E+00 

658 0.337600E+01 

659 -.188377E+01 

660 -.153699E+01 

661 0.728273E+00 

662 0.173535E+01 

663 0.862127E-01 

666 ~,321595E+01 

665 0.280656E+01 

666 0.515863E+01 

667 0.127567E+01 

668 -.271868E+01 

669 0.316677E+01 

670 -.686675E+00 

671 -.370720E+01 

672 0.685160E+00 

673 0.658657E+01 

676 -.623002E+01 

675 -.598999E+01 

676 0.536761E+01 

677 -.155600E+01 

678 -.832896E+01 

679 0.668226E+00 

680 0.308166E+01 

681 -.259356E+01 

682 0.101956E+01 

683 -.183806E+00 

686 0.916696E+00 

685 0.103868E+00 

686 ~.568762E+01 

687 -.193591E+01 

688 0.333655E+01 

689 -.607597E^01 

690 -.881921E+01 

691 0.956705E+01 

692 0.926565E+00 

693 -.905360E+01 

696 0.602092E+01 

695 0.527295E+01 

696 -.906265E+01 

697 -.233179E+01 



0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0. 000000 E+ 00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
O.OOOOOOE+00 . 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 .OOOOOOE+00 
O.OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . 000000E400 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
O.OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE + 00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 .OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE + 00 
0 . OOOOOOE + 00 
0 . OOOOOOE + 00 
0 . OOOOOOE+00 
0 . OOOOOOE + 00 
0 . OOOOOOE+00 
0 . OOOOOOE+OO 
0 . OOOOOOE+00 
0 . 000000E^ 00 
0 . OOOOOOE+00 
0 . OOOOOOE + 00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE-J 00 
0 . OOOOOOE+00 
0 . OOOOOOE + 00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
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0.880906E+01 
^99 -.595897E+01 

500 -.y^S-'iSlE+Ol 

501 0 . 78^i850E^01 

502 0.331569EI01 

503 -.686388Et01 

509 0.161638E+01 

505 0.279332E+01 

506 “.107526E+02 

507 0.952893E+01 

508 0.396753E+01 

509 -.579705E+01 

510 0.397997E+01 

511 0.292562E+01 



0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . 000000E^ 00 
0 . 000000E+ 00 
0 . 0 0 0 0 0 0 E H) 0 
0 . 000000E-» 00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
0 . OOOOOOE+00 
O.OOOOOOE+00 
0 . OOOOOOE+00 
O.OOOOOOE+00 
0 . OOOOOOE+00 
O.OOOOOOE+00 



DPT OUTPUT DATA FOR 
K REAL PART 

135 0.21869619E+03 

136 -.20999591E+03 

137 -.1<4575168E+03 

138 -.21255990E+03 

139 -.15577037E+03 

190 13250822E^03 

191 0 . 173927 02E + 03 

192 0 . 13396306E + 03 

193 0. 16901073E+03 

199 0.19780035E+03 

195 -.20230358E+03 

196 0.20827090E+03 

197 - . 19098985E+03 

198 0 .1557 1500E + 03 

199 0.22990698E+03 

150 0.18717580E+03 

151 -.13020739E+03 

152 0.21125522E+03 

153 0.16326271E^03 

159 18001019E+03 

155 -.19258061E+03 

156 -.1902/963E+03 

157 -.21912985E+03 

158 -.20591913E+03 

159 19955086E+05 

160 “.21075299E+03 

161 0.17198572E+03 

162 17657633E+03 

163 18819537E + 03 

169 -.17739966E+03 

165 18872591 E+03 

166 -.1S190575E+03 



BAUD 8 2 

IMAG PART 

-.13937286E+03 

“.20220082E+03 

- . 15299678E+03 
0 . 19105999E+03 
-.13739233E+03 
0.19138666E+03 
-.22279791E+03 

13225531E+03 
0 . 19665099E + 03 
15813335E + 03 
0.272923S3E+03 
18 137956E-t03 
0 .22690936E+03 
0 . 19765985E+03 
0 . 18650917E+03 
0. 16107190E403 
-.21927161E+03 
-.22968059E+03 

- . 1 7 587959E+03 
0.15380961E+03 
- . 13529289E+03 
0.16097937E+03 
0 . 18872530E+03 
-.22383812E+03 
0.20583072E+03 
0.22362857E+03 
0 . 1 788 1 97 3E+03 
0 . 15956850E+03 
0 . 19992966E + 03 

17915202E+03 

- . 18359022E + 03 
0. 128 15350E'^03 



MAGNITUDE 

0.25667896ET03 
0.28790259E+03 
0.21127296E+03 
0.285S0273E+03 
0.20770925E+03 
0 . 19377969E+03 
0.28269697E+03 
0. 1 9183899E+03 
0.22376578E+03 
0.25327560E+03 
0.33932995E^03 
0.27617651E+03 
0.29620269E+03 
0 .21959901E+03 
0.29217578E+03 
0 .29693380E+03 
0.25501762E+03 
0 . 31206079E+03 
0.23997571E+03 
0.23676892E»03 
0 .23535905E+03 
0 .299237 32E + 03 
0 -28592799E^03 
0 . 30381006E+03 
0 .28668237E+03 
0.30728906E+03 
0.29810397E+03 
0.23799971E+03 
0.23718896E+03 
0.25208797E+03 
0.26325732E+03 
0.22251521E403 



PHASE 

(DEG) 

-.31567688E+02 
-.13538632E+03 
13362015E + 03 - 

0. 13809989E+03 0 

-.13858707E+03 - 

0 . 13319339Et03 0 

-.52022568E+02 - 

-.93583252E+02 - 

0.90998292E+02 0 

38650635E + 02 - 

0 . 12659782E+03 0 

-.91051259ET02 - 

0 . 13015030E^03 0 

0 . 93979050E+02 0 

0 . 39667282E + 02 0 

0 .907 13165E+ 02 0 

1207 0269E + 03 - 

-.97392807E+02 - 

-.97130529E+02 - 

0 . 13998868E+03 0 

19991 098E+03 - 

0 . 1 3976907 E-» 03 0 

0.13860855E+03 0 

-.15259298E+03 - 

0.13911296E^03 0 

0 . 13330223E+03 0 

0 .961 16013E+02 0 

0. 1 3789659E + 03 0 

0 . 1 9298850E+03 0 

1 397 1 030Ef03 - 

- . 1 3579803E+ 03 - 

0 . 19983507E + 03 0 



NOTE: The DFT is windowed (WINDOW = 1) 



(RAD) 

55096029E^00 
23629370E+01 
.23321 119E+01 
.29099298E+01 
.29188013E+01 
.23237906E+01 
.90796512E+00 
.76067126E+00 
.71968163E+00 
.67958087E+00 
.22095999E+01 
.71697972E+00 
.22715511E+01 
.75835259ETOO 
.69232982E+00 
.71057868E+00 
.21066589E+01 
.82716098E+00 
.82258289E+00 
.29395369E+01 
.25291739E+01 
.29399302E+01 
.29191723E+01 
.23133116E+01 
.23907099E+01 
.23265629E+01 
.80987633E+00 
.29067997E+01 
.29868903E+01 
.23511391E+01 
.23701239E+01 
.25278992E+01 
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C. GRAPHICAL OUTPUT USING DISSPLA 

1. DFT Windowed Between ki and k2 (IWNDOW=l) 

DPT OUTPUT OP THP RPCEIVED 51GNPL 
rOR BnUD NUI1BCR 1 




DPT OUTPUT OP THE RECEIVED SIGNAL 
rOR BAUD NUIIBCR 1 




PFEO II ) 



79 



PHFISE (DEG) (IFlGNI'nJOE 

10XO-i:S.O -50.0 -C.o 0.0 45.0 ^I^.o ITi.O NO.O 0.0 5D.0 lOD.U 191 0 200.0 S0.0 JCO.O ao.D 



DFT OUTPUT or THU RECEIVED SlGNhL 

rOR BOUD NUI1BER 2 
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135,0 137.5 ItO.Q 112.5 M5.0 H7.S ISD.O 152 S 155 0 IS*. 5 IBO.O 16:. 5 Iffi 0 16V5 

rFEQ n ) 



DFT OUTPUT OF THE RECEIVED S19NRL 

rOF: BFlUD NUIIBER 2 




TFEO I t ) 
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X(LL) (FT) 



RECEIVED SIGNAL 




X-P03ITI0N 




BRl.O (LU 
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Y -POSIT I ON 



s 

ri' 

r-4 

W 

Es 

ov" 



l.< 15 I « 

BRL^? (LU 



Z-P05ITI0N 



CO 




ERi.O ILU 
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C(LU (KFT/SECI R(LU (FT) 

S.a'H] 5.0CQ5 5. DO JO 5.D315 5.0020 5.DD25 5590.0 5923.5 SDOO.O tTUD.5 0CO1.O 6001.5 GDC12.0 



SLANT RANGE TO RECEIVER 



l.D X.l 1.: 13 11 ].s I .5 J 7 1 .« 1.9 : 0 

BAUD (LU 



SPEED or SOUND 




83 



RLFHfl(LL) 

IKClfllLL) (DEQ) -0.20 -0.1S -0.10 -O.CG O.Ci5 0.C6 O.IO 0.15 



COMPRESSION mCTOR DUE TO THE IIOVIMG TX 




hNGLE BETNEEIJ RlLLl hi ID ZD 




Pqi.O ILL I 
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phase: (DEG) nflQrdHJDE 

-193.0 -L:5.0 -90.0 -6.0 0.0 '^l.O [Th.U l:)D D 0 0 2S.0 50.0 75.0 lOO.O 125.0 150.0 ITS.D 



2. DFT Not Windowed (IWNDOW=0) 



DFT OUTPUT OF THE RECEIVED 51 CURL 

rOR BRUO NUMBER 1 




DFT OUTPUT OF THE RECEIVED 51 CURL 

rOR BRUD NUMBER 1 




0.0 rs.D sc.D 75 D 100.0 iTs o ISO 0 jrs c :CD.O ZS 0 so 0 r-5 D 

FREQ U J 
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DPT OUIT'LIT Lir TIIC 1 YQi ..JPNill 

rOR BRUD NUI1BCR 2 




DPT OUTPUT OF THF RFP.EJvPD 
rOR BAUD MUIIBBR 2 



S I Gl JFIL 
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APPENDIX C. STATISTICS OF SNRour VS. SNRin ANALYSIS 
A. INPUT SNR = 0 dB 

INPUT SNRNB = 1.000 = 0.000 DB 

BAUD TYPE 1: KX = 256 SAMPLE POINTS; K = 16 TOfiES 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 1ST QUADRANT XX 



* 


NUMBER 


OF ■ 


TONES 


= 


A 






MEAN 


OF 


DFT 


REAL 


PART 


= 


97.73993 






VARIANCE 


OF 


DFT 


REAL 


PART 


= 


8919.37 






QUADRANT SNROUT 


OF 


DFT 


REAL 


PART 


= 


1.1397 = 


0.5986 


DB 


MEAN 


OF 


DFT 


IMAG 


PART 


= 


129.16291 






VARIANCE 


OF 


DFT 


IMAG 


PART 


= 


6690.81 






QUADRANT SNROUT 


OF 


DFT 


IMAG 


PART 


r 


2.5122 = 


9.0005 


DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT XX 





NUMBER 


OF 


TONES 


= 


2 






MEAN 


OF 


DFT 


REAL 


PART 


= 


-85.21500 






VARIAtJCE 


OF 


DFT 


REAL 


PART 


r 


A29 .5A 






QUADRANT SNROUT 


OF 


DFT 


REAL 


PART 




16.9056 = 


12.2803 


DB 


MEAN 


OF 


DFT 


IMAG 


PART 




82.06A99 






VARIANCE 


OF 


DFT 


IMAG 


PART 


= 


2787 .80 






QUADRANT SNROUT 


OF 


DFT 


IMAG 


PART 


= 


2.A158 = 


3.8305 


DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 3RD QUADRANT XX 





NUMBER 


OF ‘ 


TONES 


= 


6 






MEAN 


OF 


DFT 


REAL 


PART 


= 


-61.63330 






VARIANCE 


OF 


DFT 


REAL 


PART 


= 


12897 .87 






QUADRANT SNROUT 


OF 


DFT 


REAL 


PART 


— 


0.29A5 = 


-5.3089 


DB 


MEAfJ 


OF 


DFT 


IMAG 


PART 


- 


-157 .1A986 






VARIANCE 


OF 


DFT 


IMAG 


PART 


= 


1129.39 






QUADRANT SNROUT 


OF 


DFT 


IMAG 


PART 




21.8667 = 


13.3978 


DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE ATH QUADRANT XX 





NUMBER 


OF 


TONES 


= 


A 






MEAN 


OF 


DFT 


REAL 


PART 


= 


135.552A3 






VARIANCE 


OF 


DFT 


REAL 


PART 




7555. A1 






QUADRANT SNROUT 


OF 


DFT 


REAL 


PART 


= 


2.A320 = 


3.8596 


DB 


MEAN 


OF 


DFT 


IMAG 


PART 


= 


-13A. 99750 






VARIANCE 


OF 


DFT 


IMAG 


PART 


= 


11628 . 12 






QUADRANT SfJROUT 


OF 


DFT 


IMAG 


PART 




1 .5673 = 


1.9519 


DB 



XXX OVERALL (REAL + IMAG) STATISTICS FOR BAUD TYPE 1 XXX 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



32 

113.65836 
5680 .15 

2.27A3 = 3.568A DB 
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INPUT SNRf^B = 1.000 = 0.000 DB 

BAUD TYPE 2: KX = 512 SAMPLE POINTS; K = 32 TONES 



XK GIVEN THE TRANSMITTED PHASE IS IN THE 1ST QUADRANT 





NUMBER 


: OF 


TOfJES 


= 


5 






MEAN 


OF 


DFT 


REAL 


PART 


= 


137.35571 






VARIANCE 


OF 


DFT 


REAL 


PART 


= 


1A93A.A1 






QUADRANT Sf^ROUT 


OF 


DFT 


REAL 


PART 




1.2633 = 


1.0151 


DB 


MEAN 


OF 


DFT 


IMAG 


PART 


= 


263.66382 






VARIANCE 


OF 


DFT 


IMAG 


PART 


= 


5A65A.25 






QUADRANT SNROUT 


OF 


DFT 


IMAG 


PART 


= 


1.2720 = 


1 . 0AA8 


DB 



GIVEN THE TRAf^SMITTED PHASE IS IN THE 2ND QUADRANT 





NUMBER 


: OF 


TONES 


= 


8 






MEAN 


OF 


DFT 


REAL 


PART 




-126 .1A920 






VARIANCE 


OF 


DFT 


REAL 


PART 


= 


60110.95 






QUADRANT SNROUT 


OF 


DFT 


REAL 


PART 


= 


0.26A7 = 


-5.7718 


DB 


MEAN 


OF 


DFT 


IMAG 


PART 


= 


137 .62115 






VARIANCE 


OF 


DFT 


IMAG 


PART 


= 


19688.26 






QUADRANT SNROUT 


OF 


DFT 


IMAG 


PART 


= 


0.9620 = 


-0 . 168<i 


DB 



GIVEN THE TRAfiSMITTED PHASE IS IN THE 3RD QUADRANT 





fJUMBER 


OF 


TOfJES 




10 






MEAfi 


OF 


DFT 


REAL 


PART 


= 


-191 .8702A 






VARIANCE 


OF 


DFT 


REAL 


PART 


= 


196A1 .78 






QUADRAfJT SrJROUT 


OF 


DFT 


REAL 


PART 


— 


1.87A3 = 


2.7283 


DB 


MEAfJ 


OF 


DFT 


IMAG 


PART 




-288 .79A92 






VARIAtICE 


OF 


DFT 


IMAG 


PART 


= 


30312.55 






QUADRANT SfJROUT 


OF 


DFT 


IMAG 


PART 




2.751A = 


A .3956 


DB 



XX GIVEf^ THE TRANSMITTED PHASE IS IN THE ATH QUADRANT xx 



NUMBER OF TONES 
MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 

QUADRANT SfJROUT OF DFT REAL PART 



9 

15A .75328 
A0AAA.66 

0.5921 = -2.2758 DB 



MEAN OF DFT IMAG PART 

VARIANCE OF DFT IMAG PART 

QUADRANT SNROUT OF DFT IMAG PART 



-215.355A8 
21217 .86 

2.1858 = 3.3961 DB 



XXX OVERALL (REAL + IMAG) STATISTICS FOR BAUp TYPE 2 xxx 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAfJ 
BAUD VARIANCE 
BAUD SNROUT 



6A 

191 . A5129 
28251 . 36 

1 .297A = 1 .1308 DB 
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INPUT SNRNB = 1.000 = 0.000 DB 

BAUD TYPE 5: KX = 1024 SAMPLE POINTS; K = 64 TONES 



GIVEN THE TRANSMITTED PHASE IS IN THE 1ST QUADRANT xx 





NUMBER 


OF ■ 


rONES 


= 


15 






MEAN 


OF 


DFT 


REAL 


PART 


= 


295.59131 






VARIANCE 


OF 


DFT 


REAL 


PART 


= 


171068.06 






QUADRANT SNROUT 


OF 


DFT 


REAL 


PART 


= 


0.5108 = 


-2.9179 


DB 


MEAN 


OF 


DFT 


IMAG 


PART 


= 


272.87988 






VARIANCE 


OF 


DFT 


IMAG 


PART 


= 


75251.81 






QUADRANT SNROUT 


OF 


DFT 


IMAG 


PART 


= 


0.9398 = 


-0.0AA6 


DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT XX 





NUMBER 


: OF 


TONES 


- 


15 






MEAN 


OF 


DFT 


REAL 


PART 


= 


-502.71587 






VARIANCE 


OF 


DFT 


REAL 


PART 


= 


107979.00 






QUADRANT SNROUT 


OF 


DFT 


REAL 


PART 




2.5405 = 


5.6930 


DB 


MEAN 


OF 


DFT 


IMAG 


PART 




595.86084 






VARIANCE 


OF 


DFT 


IMAG 


PART 


= 


91919.19 






QUADRANT SNROUT 


OF 


DFT 


IMAG 


PART 


- 


1.6876 = 


2.2728 


DB 



XX GIVEN THE TRAr^SMITTED PHASE IS IN THE 3RD QUADRANT XX 



NUMBER OF TONES 
MEAN OF DFT REAL PART 

VARIArJCE OF DFT REAL PART 

QUADRANT SNROUT OF DFT REAL PART 



20 

-565.57617 

158015.50 

0.8458 = -0.7274 DB 



MEAN OF DFT IMAG PART 

VARIANCE OF DFT IMAG PART 

QUADRANT SNROUT OF DFT IMAG PART 



-287.26196 

110806.51 

0.7447 = -1.2801 DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 4TH QUADRANT XX 



NUMBER OF TONES 
MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 

QUADRANT SNROUT OF DFT REAL PART 



16 

528.55247 

97164.81 

1.1108= 0 . 4565 DB 



MEAN OF DFT IMAG PART 

VARIANCE OF DFT IMAG PART 

QUADRANT SrJROUT OF DFT IMAG PART 



-504.95586 

125754.50 

2.0279 = 5.0705 DB 



XXX OVERALL (REAL + IMAG) STATISTICS FOR BAUD TYPE 3 XXX 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

368.99458 

111852.19 

1.2175 = 0.8540 DB 
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INPUT SNRNB = 1.000 = 0.000 DB 

BAUD TYPE A: KX = 20A8 SAMPLE POINTS; K = 128 TONES 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 1ST QUADRANT XX 





NUMBER 


OF TONES 


= 


38 






MEAN 


OF 


DFT 


REAL 


PART 


= 


765.552^9 






VARIAtJCE 


OF 


DFT 


REAL 


PART 


= 


^480832. 37 






QUADRANT SNROUT 


OF 


DFT 


REAL 


PART 


= 


1.2189 = 


0.8596 


DB 


MEAN 


OF 


DFT 


IMAG 


PART 


= 


^87 .95532 






VARIANCE 


OF 


DFT 


IMAG 


PART 


= 


^36912.56 






QUADRANT SNROUT 


OF 


DFT 


IMAG 


PART 


= 


0.5A50 = 


-2.6363 


DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT XX 





NUMBER 


OF 


TONES 


= 


3^ 






MEAN 


OF 


DFT 


REAL 


PART 


= 


-919.^1^06 






VARIANCE 


OF 


DFT 


REAL 


PART 


= 


^73379.37 






QUADRANT SNROUT 


OF 


DFT 


REAL 


PART 


= 


1.7857 = 


2.5181 


DB 


MEAti 


OF 


DFT 


IMAG 


PART 




909 . 3^058 






VARIANCE 


OF 


DFT 


IMAG 


PART 


= 


^03385.00 






QUADRANT SNROUT 


OF 


DFT 


IMAG 


PART 


= 


2. 0*^99 = 


3.1173 


DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 3RD QUADRANT xx 





NUMBER 


OF 


TONES 


= 


35 






MEAN 


OF 


DFT 


REAL 


PART 


= 


-7^^.60986 






VARIAtJCE 


OF 


DFT 


REAL 


PART 


= 


5355^1 .62 






QUADRANT SNROUT 


OF 


DFT 


REAL 


PART 


= 


1 .0353 = 


0.1506 


DB 


MEAtJ 


OF 


DFT 


IMAG 


PART 


= 


-6^5.88306 






VARIANCE 


OF 


DFT 


IMAG 


PART 


= 


^93911 .50 






QUADRAtJT SNROUT 


OF 


DFT 


IMAG 


PART 


= 


0.8^^6 = 


-0 .733«i 


DB 



XX GIVEN THE TRANSMITTED PHASE IS lU THE *^TH QUADRANT XX 





NUMBER 


OF ■ 


TOr^ES 


r 


21 






MEAN 


OF 


DFT 


REAL 


PART 


= 


8^5.7^316 






VARIAtJCE 


OF 


DFT 


REAL 


PART 


= 


5067^7.9^ 






QUADRAf^T SrjROUT 


OF 


DFT 


REAL 


PART 




1.^115 = 


1 .^i968 


DB 


MEAN 


OF 


DFT 


IMAG 


PART 


= 


-607 .39^53 






VARIANCE 


OF 


DFT 


IMAG 


PART 


= 


^72505.12 






QUADRANT SNROUT 


OF 


DFT 


IMAG 


PART 


= 


0.7808 = 


-1 . 0746 


DB 



XXX OVERALL (REAL + IMAG) STATISTICS FOR BAUD TYPE A XXX 



TOTAL NUMBER OF POIt^TS, 2K 
BAUD MEAN 
BAUD VARIAtJCE 
BAUD StJROUT 



256 

738.25732 

^60690.81 

1.1831 = 0.7301 DB 
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INPUT SNRNB = 1.000 = 0.000 DB 

BAUD TYPE 5*. KX = A096 SAMPLE POINTS; K = 256 TONES 



GIVEN THE TRANSMITTED PHASE IS IN THE 1ST QUADRANT XX 



NUMBER OF TONES 
MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 

QUADRANT SNROUT OF DFT REAL PART 



6A 

1507.82715 

2880655.00 

0.7893 = -1.0279 DB 



MEAN OF DFT IMAG PART 

VARIANCE OF DFT IMAG PART 

QUADRANT SNROUT OF DFT IMAG PART 



167A .69238 
2577853 . 00 

1.0880 = 0.3661 DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT XX 



NUMBER OF TONES 
MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 

QUADRANT SNROUT OF DFT REAL PART 



61 

-1636 .89526 
2A06180 . 00 

1.1136 = 0.A671 DB 



MEAN OF DFT IMAG PART 

VARIANCE OF DFT IMAG PART 

QUADRANT SNROUT OF DFT IMAG PART 



13A5.8A106 
327<^A81 .00 

0.5552 = -2.5715 DB 



XX GIVEN THE TRAr^SMITTED PHASE IS IN THE 3RD QUADRANT XX 



NUMBER OF TONES = 59 

MEAN OF DFT REAL PART = -IA^7. 5^858 
VARIANCE OF DFT REAL PART = 171015A.00 
QUADRANT SNROUT OF DFT REAL PART = 1.2255 = 



0.8825 DB 



MEAti OF DFT IMAG PART = -1208.50A15 
VARIANCE OF DFT IMAG PART = 268626A.00 

QUADRANT SNROUT OF DFT IMAG PART = 0.5A37 = -2.6A65 DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE ATH QUADRANT XX 



NUMBER OF TOl^ES 
MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 

QUADRANT SNROUT OF DFT REAL PART 



72 

1315.15088 
1992730 .00 

0.8680 = -0.6150 DB 



MEAN OF DFT IMAG PART 

VARIANCE OF DFT IMAG PART 

QUADRANT SNROUT OF DFT IMAG PART 



-15A1 .86792 
293962A . 00 

0.8087 = -0.9220 DB 



XXX OVERALL (REAL + IMAG) STATISTICS FOR BAUD TYPE 5 XXX 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



512 

1A61 .01611 
252A291 . 00 

0.8A56 = -0.7285 DB 
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B. INPUT SNR = 5 dB 



INPUT SNRNB = 3.162 = 5.000 DB 

BAUD TYPE 1: KX = 256 SAMPLE POINTS; K = 16 TONES 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 1ST QUADRAf^T xx 



NUMBER OF TONES 
MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 

QUADRANT SNROUT OF DFT REAL PART 



A 

9 ^. 572^5 
2661 .9A 

3.3599 = 5.2633 DB 



MEAN OF DFT IMAG PART 

VARIANCE OF DFT IMAG PART 

QUADRANT SNROUT OF DFT IMAG PART 



112.22A91 

2098.79 

6.0008 = 7.7821 DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT XX 





[JUMBER 


OF 


TONES 


= 


2 






MEAN 


OF 


DFT 


REAL 


PART 


= 


-87.53000 






VARIAf^CE 


OF 


DFT 


REAL 


PART 


= 


135.80 






QUADRANT SNROUT 


OF 


DFT 


REAL 


PART 


= 


56.A19A = 


17 .51^3 


DB 


MEAfJ 


OF 


DFT 


IMAG 


PART 


= 


85.78000 






VARIANCE 


OF 


DFT 


IMAG 


PART 


= 


883.68 






QUADRANT SNROUT 


OF 


DFT 


IMAG 


PART 


= 


8.3263 = 


9.20^8 


DB 



XX GIVEN THE TRAfJSMITTED PHASE IS IN THE 3RD QUADRANT XX 
[DUMBER OF TONES = 6 



MEAN 


OF 


DFT 


REAL 


PART 


= 


-7A .27998 






VARIANCE 


OF 


DFT 


REAL 


PART 


= 


^080. A6 






QUADRAfJT SrmOUT 


OF 


DFT 


REAL 


PART 


= 


1 .3522 = 


1 .3103 


DB 


MEAN 


OF 


DFT 


IMAG 


PART 


= 


-127.96658 






VARIANCE 


OF 


DFT 


IMAG 


PART 


= 


357 .95 






QUADRANT SfJROUT 


OF 


DFT 


IMAG 


PART 


- 


95.8118 = 


16.6098 


DB 


XX GIVEN THE TRANSMITTED PHASE IS 




IN THE 9TH QUADRANT 


XX 




[DUMBER 


: OF 


TONES 




9 






MEAf^ 


OF 


DFT 


REAL 


PART 




115.89991 






VARIAf^CE 


OF 


DFT 


REAL 


PART 


= 


2387 .79 






QUADRANT SfJROUT 


OF 


DFT 


REAL 


PART 


= 


5.6209 = 


7.9977 


DB 


MEAf^ 


OF 


DFT 


IMAG 


PART 




-115.59292 






VARIAfjCE 


OF 


DFT 


IMAG 


PART 


= 


3679.88 






QUADRANT SNROUT 


OF 


DFT 


IMAG 


PART 




3.6278 = 


5.5965 


DB 



XXX OVERALL (REAL + IMAG) STATISTICS FOR BAUD TYPE 1 xxx 



TOTAL NUMBER OF POItJTS, 2K = 


32 


BAUD MEAf^ = 


103.52617 


BAUD VARIAf^CE = 


1796.58 


BAUD SrJROUT = 


5.9656 = 7.7565 DB 
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INPUT SNRNB = 3.162 = 5.000 DB 

BAUD TYPE 2: KX = 512 SAMPLE POINTS; K = 52 TONES 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 1ST QUADRANT xx 





NUMBER 


; OF 


TONES 


= 


5 






MEAN 


OF 


DFT 


REAL 


PART 


- 


156 .3A995 






VARIANCE 


OF 


DFT 


REAL 


PART 


= 


A717 .56 






QUADRANT SNROUT 


OF 


DFT 


REAL 


PART 


— 


5.1818 = 


7 . 1AA8 


DB 


MEAN 


OF 


DFT 


IMAG 


PART 


- 


227 . A059A 






VARIANCE 


OF 


DFT 


IMAG 


PART 


= 


17291 .69 






QUADRANT SNROUT 


OF 


DFT 


IMAG 


PAR‘T 


= 


2.9907 = 


A .7577 


DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT xx 



NUMBER OF TONES = 8 



MEAN 


OF 


DFT 


REAL 


PART = 


-150 . 082A0 






VARIANCE 


OF 


DFT 


REAL 


PART = 


19011.11 






QUADRANT SNROUT 


OF 


DFT 


REAL 


PART = 


1.18A8 = 


0.7365 


DB 


MEAN 


OF 


DFT 


IMAG 


PART = 


156 .5197A 






VARIANCE 


OF 


DFT 


IMAG 


PART = 


6250.11 






QUADRANT StmOUT 


OF 


DFT 


IMAG 


PART = 


3.9323 = 


5.9A6A 


DB 


XX GIVEN THE TRANSMITTED PHASE IS IN 


THE 3RD QUADRANT 






NUMBER 


: OF 


TOf^ES = 


10 






MEAN 


OF 


DFT 


REAL 


PART = 


-186 .9959A 






VARIANCE 


OF 


DFT 


REAL 


PART = 


6212.98 






QUADRANT SNROUT 


OF 


DFT 


REAL 


PART = 


5.6281 = 


7 .5036 


DB 


MEAN 


OF 


DFT 


IMAG 


PART = 


-2A1 .53^93 






VARIANCE 


OF 


DFT 


IMAG 


PART = 


9587 .26 






QUADRANT SNROUT 


OF 


DFT 


IMAG 


PART = 


6.0851 = 


7.8427 


DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE ATH QUADRANT XX 



NUMBER OF TONES 
MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 

QUADRANT SNROUT OF DFT REAL PART 



9 

166 .15880 
12788.79 

2.1588 = 3.3A22 DB 



MEAN OF DFT IMAG PART 

VARIANCE OF DFT IMAG PART 

QUADRANT SNROUT OF DFT IMAG PART 



-200.21880 
6707 .95 

5.9761 = 7 .76A2 DB 



XXX OVERALL (REAL + IMAG) STATISTICS FOR BAUD TYPE 2 XXX 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



6A 

186.78583 • 

893A .9A 

3.90A8 = 5.9160 DB 
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INPUT SfJRNB = 3.162 = 5.000 DB 

BAUD TYPE 3: KX = 102A SAMPLE POINTS; K = 6A TONES 



GIVEN THE TRANSMITTED PHASE IS IN THE 1ST QUADRANT 



NUMBER OF TONES 
MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 

QUADRANT SNROUT OF DFT REAL PART 



13 

32^.77A17 

5A069.19 

1.9508 = 2.9021 DB 



MEAN OF DFT IMAG PART 
VARIANCE OF DFT IMAG PART 
QUADRANT StIROUT OF DFT IMAG PART 



312.1567A 

23785.73 

A. 0967 = 6.12A3 DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT xx 





NUMBER 


OF ' 


TONES 


= 


15 






MEAN 


OF 


DFT 


REAL 


PART 


= 


-AAl . A2603 






VARIANCE 


OF 


DFT 


REAL 


PART 


= 


3A121 .99 






QUADRANT SNROUT 


OF 


DFT 


REAL 


PART 




5.7106 = 


7.5668 


DB 


MEAN 


OF 


DFT 


IMAG 


PART 


= 


380.12280 






VARIAfiCE 


OF 


DFT 


IMAG 


PART 


= 


29067 .91 






QUADRANT SNROUT 


OF 


DFT 


IMAG 


PART 


= 


A. 9709 = 


6 .96A3 


DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 3RD QUADRANT xx 





NUMBER 


OF 


TONES 


= 


20 






MEAN 


OF 


DFT 


REAL 


PART 


- 


“36A. 12720 






VARIAf^CE 


OF 


DFT 


REAL 


PART 


= 


A9952. A7 






QUADRANT SNROUT 


OF 


DFT 


REAL 


PART 




2.65A3 = 


.2395 


DB 


MEAN 


OF 


DFT 


IMAG 


PART 


- 


“320.32812 






VARIAfiCE 


OF 


DFT 


IMAG 


PART 


= 


35015.78 






QUADRANT SNROUT 


OF 


DFT 


IMAG 


PART 


= 


2.930A = 


<^.669Z 


DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE ATH QUADRANT xx 





NUMBER 


OF 


TOflES 


= 


16 






MEAN 


OF 


DFT 


REAL 


PART 


= 


3A3. 55005 






VARIAI^CE 


OF 


DFT 


REAL 


PART 


= 


30738.97 






QUADRANT SNROUT 


OF 


DFT 


REAL 


PART 


— 


3.8396 = 


5.8A29 


DB 


MEAN 


OF 


DFT 


IMAG 


PART 




-AA2. A60A5 






VARIAfiCE 


OF 


DFT 


IMAG 


PART 


= 


3973A .58 






QUADRANT SNROUT 


OF 


DFT 


IMAG 


PART 


= 


A. 9270 = 


6,9258 


DB 



XXX OVERALL (REAL + IMAG) STATISTICS FOR BAUD TYPE 3 XXX 



TOTAL [DUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIAfiCE 
BAUD SimOUT 



128 

366 . 16089 
35357 .58 

3.7919 = 5.7886 DB 
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INPUT SNRNB = 3.162 = 5.000 DB 

BAUD TYPE A: KX = 20A8 SAMPLE POINTS; K = 128 TONES 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 1ST QUADRANT xx 



NUMBER OF TONES 
MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 

QUADRANT SNROUT OF DFT REAL PART 



38 

75A . 15A30 
151360.81 

3.7576 = 5.7A91 DB 



MEAN OF DFT IMAG PART 

VARIANCE OF DFT IMAG PART 

QUADRANT SNROUT OF DFT IMAG PART 



598.70117 

138265.00 

2.592^4 = A. 1371 DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT XX 



NUMBER OF TONES 
MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 

QUADRANT SNROUT OF DFT REAL PART 



3A 

-8A2. 0A663 
1503A9.12 

A. 7160 = 6.7357 DB 



MEAN OF DFT IMAG PART 
VARIAfJCE OF DFT IMAG PART 
QUADRANT SNROUT OF DFT IMAG PART 



83A. 70776 
127103.50 

5.A817 = 7.3891 DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 3RD QUADRANT xx 



NUMBER OF TONES 
MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 

QUADRANT SNROUT OF DFT REAL PART 



35 

“7A1 .81006 
168815.31 

3.2597 = 5.1317 DB 



MEAN OF DFT IMAG PART 

VARIANCE OF DFT IMAG PART 

QUADRANT SNROUT OF DFT IMAG PART 



-688 . 02A17 
155789.25 

3.0586 = 9.8267 DB 



XX GIVEN THE TRAtJSMITTED PHASE IS IN THE 9TH QUADRANT XX 



NUMBER OF TONES 
MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 

QUADRANT SNROUT OF DFT REAL PART 



21 

799.97021 

159680.31 

9 . 0027 = 6 . 0235 DB 



MEAN OF DFT IMAG PART 

VARIANCE OF DFT IMAG PART 

QUADRANT SNROUT OF DFT IMAG PART 



-665.15771 
199759 . 06 

2.9599 = 9.7097 DB 



XXX OVERALL (REAL + IMAG) STATISTICS FOR BAUD TYPE 9 XXX 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



256 

739.13799 ’ 

195979.56 

3.7555 = 5.7965 DB 
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INPUT SNRNB = 3.162 = 5.000 DB 

BAUD TYPE 5: KX = A096 SAMPLE POINTS; K = 256 TONES 



GIVEN THE TRANSMITTED PHASE IS IN THE 1ST QUADRANT 



NUMBER OF TONES 
MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 

QUADRANT SNROUT OF DFT REAL PART 



6A 

1508.21289 

9156A9.00 

2.A8A5 = 5.9520 DB 



MEAN OF DFT IMAG PART 

VARIANCE OF DFT IMAG PART 

QUADRANT SNROUT OF DFT IMAG PART 



1605.33105 
8232A9 .25 

3.130A = A. 9560 DB 



XX GIVEN THE TRAFJSMITTED PHASE IS IN THE 2ND QUADRANT xx 



NUMBER OF TONES 
MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 

QUADRAtJT SNROUT OF DFT REAL PART 



61 

-1585.36865 

762098.57 

3.2980 = 5.1825 DB 



MEAN OF DFT IMAG PART 

VARIANCE OF DFT IMAG PART 

QUADRANT SNROUT OF DFT IMAG PART 



1A15. 77759 
1035021 .56 

1.9366 = 2.870A DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 3RD QUADRANT xx 





NUMBER 


OF 


TOriES 


= 


59 






MEAN 


OF 


DFT 


REAL 


PART 


= 


-1A73. 39185 






VARIANCE 


OF 


DFT 


REAL 


PART 


r 


558915.81 






QUADRANT SNROUT 


OF 


DFT 


REAL 


PART 


— 


A. 0232 = 


6 . 0512 


DB 


MEAN 


OF 


DFT 


IMAG 


PART 


= 


-13A2. A5337 






VARIANCE 


OF 


DFT 


IMAG 


PART 


= 


8A5222.31 






QUADRANT SNROUT 


OF 


DFT 


IMAG 


PART 


= 


2.1322 = 


3.2883 


DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE ATH QUADRANT xx 





NUMBER 


OF ■ 


TONES 


= 


72 






MEAN 


OF 


DFT 


REAL 


PART 


= 


lAOl .733A0 






VARIANCE 


OF 


DFT 


REAL 


PART 


= 


652763. 9A 






QUADRANT SNROUT 


OF 


DFT 


REAL 


PART 




3 .1052 = 


A .9209 


DB 


MEAN 


OF 


DFT 


IMAG 


PART 


- 


-1526.53638 






VARIANCE 


OF 


DFT 


IMAG 


PART 


- 


927957 .37 






QUADRANT SNROUT 


OF 


DFT 


IMAG 


PART 


= 


2.5112 = 


3.9989 


DB 



XXX OVERALL (REAL + IMAG) STATISTICS FOR BAUD TYPE 5 xxx 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAfJ 
BAUD VARIANCE 
BAUD SNROUT 



512 

1A35. 02075 » 

799351 .69 

2.751A = A. 3956 DB 
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C. INPUT SNR = 10 dB 



INPUT SNRNB = 10.000 = 10.000 DB 

BAUD TYPE 1: KX = 256 SAMPLE POIfJTS; K = 16 TONES 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 1ST QUADRANT xx 



NUMBER OF TONES 
MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 

QUADRANT SNROUT OF DFT REAL PART 



A 

92.78A91 
8A0 .8A 

10.2386 = 10.1029 DB 



MEAN OF DFT IMAG PARI 

VARIANCE OF DFT IMAG PART 

QUADRANT SNROUT OF DFT IMAG PART 



102.71296 

663.56 

15.8989 = 12.0137 DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT XX 



NUMBER OF TONES = 2 

MEAN OF DFT REAL PART = -88.83999 

VARIANCE OF DFT REAL PART = 92.97 

QUADRANT SNROUT OF DFT REAL PART = 183.6700 = 22.6909 DB 



MEAN OF DFT IMAG PART 

VARIANCE OF DFT IMAG PART 

QUADRANT SNROUT OF DFT IMAG PART 



87 .89000 
278.95 

27.6602 = 19.9185 DB 



XX GIVEN THE TRArJSMITTED PHASE IS IN THE 3RD QUADRAtJT xx 



NUMBER OF TONES 
MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 

QUADRANT SNROUT OF DFT REAL PART 



6 

-81 .36990 
1289.56 

5.1399 = 7.1099 DB 



MEAfi OF DFT IMAG PART 
VARIANCE OF DFT IMAG PART 
QUADRANT SNROUT OF DFT IMAG PART 



-111.58658 

113.19 

110.0076 = 20.9192 DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 9TH QUADRANT XX 



NUMBER OF TONES 
MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 

QUADRANT SNROUT OF DFT REAL PART 



9 

109.76790 

756.27 

19.5137 = 11.6178 DB 



MEAN OF DFT IMAG PART 
VARIAfJCE OF DFT IMAG PART 
QUADRANT SNROUT OF DFT IMAG PART 



-109.59991 
1169 . 08 

9.3981 = 9.7309 DB 



XXX OVERALL (REAL + IMAGl STATISTICS FOR BAUD TYPE 1 XXX 



TOTAL fJUMDER OF POINTS, 2K 
BAUD MEAtJ 
BAUD VARIANCE 
BAUD SNROUT 



32 

97 .82896 
568 . 06 

16 .8977 = 12.2659 DB 
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INPUT SNRNB = 10.000 = 10.000 DB 

BAUD TYPE 2: KX = 512 SAMPLE POINTS; K = 52 TONES 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 1ST QUADRANT XX 





NUMBER 


OF • 


rONES 


= 


5 






MEAN 


OF 


DFT 


REAL 


PART 


= 


167.05995 






VARIANCE 


OF 


DFT 


REAL 


PART 


= 


1A90.65 






QUADRANT SNROUT 


OF 


DFT 


REAL 


PART 


— 


18.7182 = 


12.7226 


DB 


MEAN 


OF 


DFT 


IMAG 


PART 


f 


207.00189 






VARIANCE 


OF 


DFT 


IMAG 


PART 


= 


5A71 .27 






QUADRANT SNROUT 


OF 


DFT 


IMAG 


part 


= 


7.8518 = 


8.9586 


DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT xx 





NUMBER 


OF ■ 


TONES 


= 


8 






MEAN 


OF 


DFT 


REAL 


PART 


= 


-165.52A95 






VARIANCE 


OF 


DFT 


REAL 


PART 


= 


6010. A2 






QUADRANT SNROUT 


OF 


DFT 


REAL 


PART 


= 


A.AA90 = 


6 . A826 


DB 


MEAN 


OF 


DFT 


IMAG 


PART 


= 


167 .12569 






VARIAriCE 


OF 


DFT 


IMAG 


PART 


= 


1969.51 






QUADRANT SNROUT 


OF 


DFT 


IMAG 


PART 




1A.1815 = 


11.5172 


DB 


XX GIVEN THE TRANSMITTED PHASE IS 


I 


N THE 5RD QUADRANT 


XX 




r^UMBER 


OF 


TONES 


= 


10 






MEAfJ 


OF 


DFT 


REAL 


PART 


= 


-ISA. 25990 






VARIAfiCE 


OF 


DFT 


REAL 


PART 


= 


1965.78 






QUADRANT SNROUT 


OF 


DFT 


REAL 


PART 


= 


17.271A = 


12.5755 


DB 


MEAN 


OF 


DFT 


IMAG 


PART 




-21A. 98996 






VARIANCE 


OF 


DFT 


IMAG 


PART 


= 


505A.56 






QUADRANT SNROUT 


OF 


DFT 


IMAG 


PART 


= 


15.251A = 


1— ' 

00 

ro 


DB 



XX GIVEN THE TRANSMITTED PHASE IS lU THE ATH QUADRANT XX 



[DUMBER OF TONES 
MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 

QUADRANT SNROUT OF DFT REAL PART 



9 

172.58655 
A0A6 .22 

7 .5615 = 8.6696 DB 



MEAr^ OF DFT IMAG PART 

VARIAr^CE OF DFT IMAG PART 

QUADRAfJT SNROUT OF DFT IMAG PART 



-191 .69989 
2121.56 

17.5216 = 12.5859 DB 



XXX OVERALL (REAL + IMAG) STATISTICS FOR BAUD TYPE 2 XXX 



TOTAL NUMBER OF POIfiTS, 2K 
BAUD MEAN 
BAUD VARIAf^CE 
BAUD SNROUT 



6A 

18A. 16551 > 

2826 .25 

12.0005 = 10.7920 DB 



98 



INPUT SNRMB = 10.000 = 10.000 DB 

BAUD TYPE 3: KX = 102A SAMPLE POINTS; K = 6A TONES 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 1ST QUADRANT xx 



MEAN 
VARIANCE 
QUADRANT SNROUT 



NUMBER OF TONES 
OF DFT REAL PART 
OF DFT REAL PART 
OF DFT REAL PART 



13 

3A1 .20361 
17078.35 

6.8168 = 8.3358 DB 



MEAN OF DFT IMAG PART 

VARIANCE OF DFT IMAG PART 

QUADRANT SNROUT OF DFT IMAG PART 



33A .22290 
7515.77 

1A.8627 = 11.7210 DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT xx 





NUMBER 


OF ■ 


TONES 


= 


15 






MEAN 


OF 


DFT 


REAL 


PART 


= 


-A06 .96606 






VARIANCE 


OF 


DFT 


REAL 


PART 


= 


10785.23 






QUADRANT SNROUT 


OF 


DFT 


REAL 


PART 


~ 


15.3563 = 


11 .8629 


DB 


MEAN 


OF 


DFT 


IMAG 


PART 


= 


372. A1919 






VARIANCE 


OF 


DFT 


IMAG 


PART 


= 


9200.02 






QUADRANT SNROUT 


OF 


DFT 


IMAG 


PART 


= 


15.0756 = 


11.7828 


DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 3RD QUADRANT XX 





NUMBER 


OF ■ 


TONES 


- 


20 






MEAN 


OF 


DFT 


REAL 


PART 


= 


-363.32935 






VARIAfJCE 


OF 


DFT 


REAL 


PART 


= 


1579A.66 






QUADRANT SNROUT 


OF 


DFT 


REAL 


PART 




8.3578 = 


9.2209 


DB 


MEAN 


OF 


DFT 


IMAG 


PART 


— 


-338.93921 






VARIANCE 


OF 


DFT 


IMAG 


PART 


r 


11063.66 






QUADRANT SNROUT 


OF 


DFT 


IMAG 


PART 


- 


10.3835 = 


10.1634 


DB 


XX GIVEN THE TRANSMITTED PHASE IS 


I 


N THE ATH QUADRANT 


XX 




NUMBER 


: OF 


TONES 




16 






MEAN 


OF 


DFT 


REAL 


PART 


= 


351 .993A1 






VARIANCE 


OF 


DFT 


REAL 


PART 


= 


9729.22 






QUADRANT SNROUT 


OF 


DFT 


REAL 


PART 




12.73A8 = 


11 .0A99 


DB 


MEAN 


OF 


DFT 


IMAG 


PART 


= 


-A07 .30591 






VARIANCE 


OF 


DFT 


IMAG 


PART 


= 


12552. A8 






QUADRANT SNROUT 


OF 


DFT 


IMAG 


PART 




13.216A = 


11 .2111 


DB 



XXX OVERALL (REAL + IMAG) STATISTICS FOR BAUD TYPE 3 XXX 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

36A .57397 
11176.81 
11 .8920 = 



10.7525 DB 



99 



INPUT SNRNB = 10.000 = 10.000 DB 

BAUD TYPE A: KX = 20A8 SAMPLE POINTS; K = 128 TONES 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 1ST QUADRANT xx 





NUMBER 


OF ■ 


TONES 


= 


38 






MEAN 


OF 


DFT 


REAL 


PART 


= 


7A7 .6960A 






VARIAfJCE 


OF 


DFT 


REAL 


PART 


= 


A7A69 .11 






QUADRANT SNROUT 


OF 


DFT 


REAL 


PART 




11.7771 = 


10.710^ 


DB 


MEAN 


OF 


DFT 


IMAG 


PART 


= 


660.9670A 






VARIANCE 


OF 


DFT 


IMAG 


PART 


= 


A3777 .00 






QUADRANT SNROUT 


OF 


DFT 


IMAG 


PART 


= 


9.9796 = 


9.9911 


DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT XX 



NUMBER OF TONES 
MEAN OF DFT REAL PART 
VARIAfJCE OF DFT REAL PART 
QUADRANT SMROUT OF DFT REAL PART 



3A 

-798.57202 

A7932.87 

13.30AA = 11.2399 DB 



MEAN OF DFT IMAG PART 

VARIANCE OF DFT IMAG PART 

QUADRANT SNROUT OF DFT IMAG PART 



792.75A39 

39958.12 

15.7280 = 11.9667 DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 3RD QUADRANT XX 



NUMBER OF TONES 
MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 

QUADRANT SfiROUT OF DFT REAL PART 



35 

-7A0 . 16992 
53078.29 

10.3216 = 10.1375 DB 



MEAN OF DFT IMAG PART 
VARIAfICE OF DFT IMAG PART 
QUADRANT SNROUT OF DFT IMAG PART 



-711 .70A35 
A9063.26 

10.3239 = 10.138A DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE ATH QUADRANT XX 





NUMBER 


OF 


TONES 


= 


21 






MEAN 


OF 


DFT 


REAL 


PART 


= 


773.37500 






VARIANCE 


OF 


DFT 


REAL 


PART 


= 


50136 .A3 






QUADRANT SNROUT 


OF 


DFT 


REAL 


PART 


- 


11.9296 = 


10.7663 


DB 


MEAfJ 


OF 


DFT 


IMAG 


PART 


= 


-697.57031 






VARIAfICE 


OF 


DFT 


IMAG 


PART 


= 


A7513.27 






QUADRAfJT SNROUT 


OF 


DFT 


IMAG 


PART 


= 


10.2A1A = 


10.1036 


DB 



XXX OVERALL (REAL + IMAG) STATISTICS FOR BAUD TYPE A xxx 



TOTAL NUMBER OF POIfITS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



256 

739.60791 * 

A5891 . A1 

11.9199 = 10.7627 DB 



100 



INPUT SNRNB = 10.000 = 10.000 DB 

BAUD TYPE 5: KX = A096 SAMPLE POINTS; K = 256 TONES 



XX given the TRANSMITTED PHASE IS IN THE 1ST QUADRANT XX 



NUMBER OF TONES 
MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 

QUADRANT SNROUT OF DFT REAL PART 



69 

1508.31250 

292952.31 

7.7791 = 8.9093 DB 



MEAN OF DFT IMAG PART 

VARIANCE OF DFT IMAG PART 

QUADRANT SNROUT OF DFT IMAG PART 



1566.18359 

265130.50 

9.2518 = 9.6623 DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT XX 



NUMBER OF TONES 
MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 

QUADRANT SNROUT OF DFT REAL PART 



61 

-1556.26519 

292012.25 

10.0076 = 10.0033 DB 



MEAN OF DFT IMAG PART 
VARIANCE OF DFT IflAG PART 
QUADRANT SNROUT OF DFT IMAG PART 



1955.11865 
327317 . 06 

6.9689 = 8.1083 DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 3RD QUADRANT XX 



NUMBER OF TONES 
MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 

QUADRANT SNROUT OF DFT REAL PART 



59 

-1987.87915 

169619.31 

13.0519 = 11.1567 DB 



MEAN OF DFT IMAG PART 

VARIANCE OF DFT IMAG PART 

QUADRANT SNROUT OF DFT IMAG PART 



-1917 .67261 
265298 .00 

7.5770 = 8.7950 DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 9TH QUADRANT XX 



NUMBER OF TONES 
MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 

QUADRANT SNROUT OF DFT REAL PART 



72 

1950.37061 
202109 .12 

10.9089 = 10.1738 DB 



MEAN OF DFT IMAG PART 

VARIANCE OF DFT IMAG PART 

QUADRANT SNROUT OF DFT IMAG PART 



-1517.82886 

292801.12 

7.8682 = 8.9587 DB 



XXX OVERALL (REAL + IMAG) STATISTICS FOR BAUD TYPE 5 XXX 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



512 

1995.31201 • 

253713.62 

8.8129 = 9.9512 DB 



101 



D. INPUT SNR = 15 dB 



INPUT SNRNB = 31.625 = 15.000 DB 

BAUD TYPE 1: KX = 256 SAMPLE POINTS; K = 16 TONES 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 1ST QUADRANT XX 

NUMBER OF TONES = A 
MEAtJ OF DFT REAL PART = 91.78998 

VARIANCE OF DFT REAL PART = 265.70 

QUADRAt^T StmOUT OF DFT REAL PART = 31.7 097 = 15.0119 DB 

MEAN OF DFT IMAG. PART = 97.37A9A 

VARIANCE OF DFT IMAG PART = 209.96 

QUADRANT SNROUT OF DFT IMAG PART = A5.1599 = 16.5A75 DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT XX 

NUMBER OF TONES = 2 

MEAt^ OF DFT REAL PART = “89.56999 

VARIAf^CE OF DFT REAL PART = 13.62 

QUADRANT SNROUT OF DFT REAL PART = 588.8657 = 27.7002 DB 

MEAN OF DFT IMAG PART = 89.00999 

VARIANCE OF DFT IMAG PART = 88.18 

QUADRANT SNROUT OF DFT IMAG PART = 89.8^86 = 19.5351 DB 



XX GIMEU THE TRAfiSMITTED PHASE IS IN THE 3RD QUADRANT XX 

r^UMBER OF TONES = 6 

MEAN OF DFT REAL PART = -85.36328 

VARIANCE OF DFT REAL PART = ^07.69 

QUADRAt^T StmOUT OF DFT REAL PART = 17.8737 = 12.5221 DB 

MEAN OF DFT IMAG PART = -102.36160 

VARIANCE OF DFT IMAG PART = 35.65 

QUADRANT SNROUT OF DFT IMAG PART = 29^.1055 = 2A.6850 DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE ^TH QUADRANT XX 
[DUMBER OF TONES = 

MEAfJ OF DFT REAL PART = 98.517^6 

VARIAtICE OF DFT REAL PART = 238.91 

QUADRAfiT StmOUT OF DFT REAL PART = ^0.62^1 = 16.0878 DB 

MEAN OF DFT IMAG PART = -98.39^90 

VARIANCE OF DFT IMAG PART = 366.81 

QUADRANT StmOUT OF DFT IMAG PART = 26.59^1 = 1^.2151 DB 



XXX OVERALL (REAL + IMAG) STATISTICS FOR BAUD TYPE 1 xxx 



TOTAL tiUMBER OF POItUS, 2K 
BAUD MEAfi 
BAUD VARIAtiCE 
BAUD SNROUT 



32 

9^.61931 

179 

^9.89^1 = 16.9805 DB 



102 



INPUT Sf^RNB = 31.623 = 15.000 DB 

BAUD TYPE 2: KX = 512 SAMPLE POINTS; K = 32 TONES 



KX GIVEN THE TRANSMITTED PHASE IS IN THE 1ST QUADRANT XX 





NUMBER 


OF 


TONES 


= 


5 






MEAN 


OF 


DFT 


REAL 


PART 


- 


173.05991 






VARIANCE 


OF 


DFT 


REAL 


PART 


= 


A69 .50 






QUADRANT SNROUT 


OF 


DFT 


REAL 


PART 


— 


63.7903 = 


18 . 0A75 


DB 


MEAN 


OF 


DFT 


IMAG 


PART 


— 


195.51991 






VARIAfiCE 


OF 


DFT 


IMAG 


PART 


= 


1735.51 






QUADRANT SNROUT 


OF 


DFT 


IMAG 


PART 


= 


22.0270 = 


13. A295 


DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT XX 



MEAN 
VARIANCE 
QUADRANT SNROUT 



NUMBER OF TONES 
OF DFT REAL PART 
OF DFT REAL PART 
OF DFT REAL PART 



8 

-171 . 09995 
1899.12 

15.A151 = 11.8795 DB 



MEAN OF DFT IMAG PART 

VARIANCE OF DFT IMAG PART 

QUADRANT SNROUT OF DFT IMAG PART 



173.09988 

623.56 

^8.052^ = 16.8171 DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 3RD QUADRANT xx 



NUMBER OF TOtiES 
MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 

QUADRANT SNROUT OF DFT REAL PART 



10 

-182.69992 
622. A3 

53.6277 = 17.2939 DB 



MEAN OF DFT IMAG PART 

VARIANCE OF DFT IMAG PART 

QUADRANT SNROUT OF DFT IMAG PART 



-200.01991 
961 . OA 

A1.6297 = 16.19A0 DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE ATH QUADRANT xx 





NUMBER 


OF 


TONES 


= 


9 






MEAN 


OF 


DFT 


REAL 


PART 


= 


176.17767 






VARIANCE 


OF 


DFT 


REAL 


PART 


= 


1279.21 






QUADRANT StiROUT 


OF 


DFT 


REAL 


PART 


— 


2^^.2639 = 


13.8^^96 


DB 


MEAN 


OF 


DFT 


IMAG 


PART 


— 


-186 .88879 






VARIANCE 


OF 


DFT 


IMAG 


PART 


= 


670 . 12 






QUADRANT SNROUT 


OF 


DFT 


IMAG 


PART 


= 


52.1213 = 


17 . 1702 


DB 



XXX OVERALL (REAL + IMAG) STATISTICS FOR BAUD TYPE 2 XXX 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



6A 

182.67627 ’ 

89A . OA 

37.3256 = 15.7201 DB 



103 



INPUT SNRNB = 31.623 = 15.000 DB 

BAUD TYPE 3: KX = 102A SAMPLE POINTS; K = 6A TONES 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 1ST QUADRANT xx 



NUMBER OF TONES 
MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 

QUADRANT SNROUT OF DFT REAL PART 



13 

350 . A3823 
5392.23 

22.77A8 = 13.5795 DB 



MEAN OF DFT IMAG PART 

VARIANCE OF DFT IMAG PART 

QUADRANT SNROUT OF DFT IMAG PART 



396.65332 

2375.91 

50.5885 = 17.0905 DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT xx 



NUMBER OF TONES 
MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 

QUADRANT SNROUT OF DFT REAL PART 



15 

-387 .58618 
3903.85 

99.1332 = 16.9976 DB 



MEAN OF DFT IMAG PART 

VARIANCE OF DFT IMAG PART 

QUADRANT SNROUT OF DFT IMAG PART 



368 . 07251 
2911.29 

96.5359 = 16.6779 DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 3RD QUADRANT xx 



NUMBER OF TONES 
MEAN OF DFT REAL PART 

VARIAf^CE OF DFT REAL PART 

QUADRANT SNROUT OF DFT REAL PART 



20 

-362.88928 

9993.27 

26.3725 = 19.2115 DB 



MEAN OF DFT IMAG PART 

VARIAtiCE OF DFT IMAG PART 

QUADRANT Sf^ROUT OF DFT IMAG PART 



-399.39929 

3993.32 

39.9957 = 15.9339 DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 9TH QUADRANT XX 





flUMBER 


OF 


TOf^ES 


= 


16 






MEAr] 


OF 


DFT 


REAL 


PART 


= 


356 .73682 






VARIAr^CE 


OF 


DFT 


REAL 


PART 


= 


3080.93 






QUADRANT SNROUT 


OF 


DFT 


REAL 


PART 


= 


91.3060 = 


16.1601 


DB 


MEAN 


OF 


DFT 


IMAG 


PART 


— 


-387 .59321 






VARIAf^CE 


OF 


DFT 


IMAG 


PART 


= 


3962.38 






QUADRANT SNROUT 


OF 


DFT 


IMAG 


PART 


= 


37.9039 = 


15.7868 


DB 



XXX OVERALL (REAL + IMAG) STATISTICS FOR BAUD TYPE 3 xxx 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

363.68092 

3531.69 

37.9510 = 15.7396 DB 



104 



INPUT SNRNB = 31.623 = 15.000 DB 

BAUD TYPE A: KX = 20A8 SAMPLE POINTS; K = 128 TONES 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 1ST QUADRANT xx 



NUMBER OF TONES 
MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 

QUADRANT SNROUT OF DFT REAL PART 



38 

. 1A038 
1A8A3. A6 

37 .3056 = 15.7177 DB 



MEAN OF DFT IMAG PART 

VARIANCE OF DFT IMAG PART 

QUADRANT SNROUT OF DFT IMAG PART 



696 .01^4^0 
13916.01 

3A.811A = 15.A172 DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT xx 



NUMBER OF TONES 
MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 

QUADRANT SNROUT OF DFT REAL PART 



3 ^ 

-77^.}2231 
15A03. A6 

38.90A6 = 15.9000 DB 



MEAN OF DFT IMAG PART 

VARIANCE OF DFT IMAG PART 

QUADRANT SNROUT OF DFT IMAG PART 



769 .157^7 
12525.41 

47.2322 = 16.7424 DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 3RD QUADRANT XX 



NUMBER OF TONES 
MEAN OF DFT REAL PART 

VARIAfiCE OF DFT REAL PART 

QUADRANT SNROUT OF DFT REAL PART 



35 

-739.27295 
16651 .30 

32.8217 = 15.1616 DB 



MEAN OF DFT IMAG PART 

VARIANCE OF DFT IMAG PART 

QUADRANT SNROUT OF DFT IMAG PART 



-725.06152 
15437 .08 

34.0553 = 15.3218 DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 4TH QUADRANT xx 



NUMBER OF TONES 
MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 

QUADRANT SNROUT OF DFT REAL PART 



21 

758.77905 
15687 .69 

36.7005 = 15.6467 DB 



MEAN OF DFT IMAG PART 

VARIANCE OF DFT IMAG PART 

QUADRANT SNROUT OF DFT IMAG PART 



-715.81812 

15135.91 

33.8530 = 15.2960 DB 



XXX OVERALL (REAL + IMAG) STATISTICS FOR BAUD TYPE 4 xxx 



TOTAL NUMBER OF POirjTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



256 

739.90381 
14483.24 
37 .7994 = 



15.7748 DB 



105 



INPUT SNRMB = 31.623 = 15.000 DB 

BAUD TYPE 5: KX = A096 SAMPLE POINTS; K = 256 TONES 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 1ST QUADRANT xx 



NUMBER OF TONES 
MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 

QUADRANT SNROUT OF DFT REAL PART 



6A 

1508.39551 

9A51A.62 

2A.0731 = 13.8153 DB 



MEAN OF DFT IMAG PART 

VARIANCE OF DFT IMAG PART 

QUADRANT SNROUT OF DFT IMAG PART 



15AA. 17187 
87073.12 

27.38A6 = 1A.3751 DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT xx 





NUMBER 


OF ■ 


TONES 




61 






MEAf] 


OF 


DFT 


REAL 


PART 


= 


-1539.88A03 






VARIANCE 


OF 


DFT 


REAL 


PART 


- 


77501 .AA 






QUADRANT SNROUT 


OF 


DFT 


REAL 


PART 


— 


30.5961 = 


14.8567 


DB 


MEAN 


OF 


DFT 


IMAG 


PART 


= 


1A77 .20898 






VARIAfJCE 


OF 


DFT 


IMAG 


PART 


- 


103902.56 






QUADRANT SNROUT 


OF 


DFT 


IMAG 


PART 


= 


21 .0018 = 


15.2226 


DB 



XX GIVEf^ THE TRANSMITTED PHASE IS IN THE 3RD QUADRANT XX 





NUMBER 


OF • 


rONES 


= 


59 






MEAf] 


OF 


DFT 


REAL 


PART 


= 


“1A96 .03369 






VARIAfJCE 


OF 


DFT 


REAL 


PART 


= 


53A83. A8 






QUADRANT SNROUT 


OF 


DFT 


REAL 


PART 


= 


A1.8A69 = 


16.2166 


DB 


MEAN 


OF 


DFT 


IMAG 


PART 


= 


-1A59. 93311 






VARIAtiCE 


OF 


DFT 


IMAG 


PART 


= 


83136 . 06 






QUADRANT SNROUT 


OF 


DFT 


IMAG 


PART 


= 


25.6375 = 


1A.0888 


DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE ATH QUADRANT xx 





NUMBER 


; OF 


TONES 


r 


72 






MEAtl 


OF 


DFT 


REAL 


PART 


= 


1A77 .83325 






VARIANCE 


OF 


DFT 


REAL 


PART 


r 


65A98 .62 






QUADRANT SNROUT 


OF 


DFT 


REAL 


PART 


X 


33.3AA1 = 


15.2302 


DB 


MEAf] 


OF 


DFT 


IMAG 


PART 


= 


-1512.91821 






VARIAfiCE 


OF 


DFT 


IMAG 


PART 


= 


92605.37 






QUADRANT SNROUT 


OF 


DFT 


IMAG 


PART 


= 


2A.7169 = 


13.9299 


DB 



XXX OVERALL (REAL + IMAG) STATISTICS FOR BAUD TYPE 5 XXX 



TOTAL NUMBER OF POIf^TS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



512 

1502.23193 ' 

81161.50 

27 .8051 = lA. AA12 DB 



106 



E. INPUT SNR = 20 dB 



INPUT SNRNB = 100.000 = 20.000 DB 

BAUD TYPE 1: KX = 256 SAMPLE POINTS; K = 16 TONES 



GIVEN THE TRANSMITTED PHASE IS IN THE 1ST QUADRANT XX 



NUMBER OF TONES 
MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 

QUADRANT SNROUT OF DFT REAL PART 



91 .237A3 
89.15 

98.9206 = 19.9529 DB 



MEAN OF DFT IMAG PART 

VARIANCE OF DFT IMAG PART 

QUADRANT SNROUT OF DFT IMAG PART 



99.37299 

66.91 

139.1027 = 21.2799 DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT XX 



NUMBER OF TOfJES = 2 

MEAN OF DFT REAL PART = -89.98000 

VARIANCE OF DFT REAL PART = 9.32 

QUADRANT SNROUT OF DFT REAL PART = xxxxxx = 32.7263 DB 



MEAN OF DFT IMAG PART 

VARIANCE OF DFT IMAG PART 

QUADRANT SNROUT OF DFT IMAG PART 



89.66999 
27 .90 

288.1609 = 29.5963 DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 3RD QUADRANT XX 



NUMBER OF TONES 
MEAN OF DFT REAL PART 

VARIAfJCE OF DFT REAL PART 

QUADRANT SNROUT OF DFT REAL PART 



6 

-87.61157 
129 . 06 

59.9730 = 17 .7932 DB 



MEAN OF DFT IMAG PART 

VARIANCE OF DFT IMAG PART 

QUADRANT SNROUT OF DFT IMAG PART 



-97.17159 

11.23 

890.5669 = 29.2957 DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 9TH QUADRANT XX 



NUMBER OF TONES 
MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 

QUADRANT SNROUT OF DFT REAL PART 



9 

95.00795 

75.29 

119.8836 = 20.7876 DB 



MEAN OF DFT IMAG PART 

VARIANCE OF DFT IMAG PART 

QUADRANT SNROUT OF DFT IMAG PART 



-99.95790 

116.38 

77.9779 = 18.8917 DB 



XXX OVERALL (REAL + IMAG) STATISTICS FOR BAUD TYPE 1 xxx 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



32 

92.82197 
56 .79 

151 .7199 = 21 .8109 DB 
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INPUT SfJRNB = 100.000 = 20.000 DB 

BAUD TYPE 2: KX = 512 SAMPLE POINTS; K = 32 TONES 



GIVEN THE TRANSMITTED PHASE IS IN THE 1ST QUADRANT 



NUMBER OF TONES 
MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 

QUADRANT SNROUT OF DFT REAL PART 



5 

176 . A5985 
1A7 .68 

210.8AA2 = 23.2396 DB 



MEAN OF DFT IMAG PART = 189.09990 

VARIANCE OF DFT IMAG PART = 5A9.10 

QUADRANT SNROUT OF DFT IMAG PART = 65.1219 = 18.1373 DB 



GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT XX 





NUMBER 


OF ■ 


TONES 


= 


8 






MEAN 


OF 


DFT 


REAL 


PART 


= 


-175.5A995 






VARIANCE 


OF 


DFT 


REAL 


PART 


= 


599.80 






QUADRANT SNROUT 


OF 


DFT 


REAL 


PART 


— 


51.2631 = 


17.0980 


DB 


MEAN 


OF 


DFT 


IMAG 


PART 


- 


176 . AA992 






VARIAtiCE 


OF 


DFT 


IMAG 


PART 


= 


197 .26 






QUADRANT SNROUT 


OF 


DFT 


IMAG 


PART 


= 


157.8330 = 


21 .9820 


DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 3RD QUADRANT xx 



NUMBER OF TOfiES 
MEAfJ OF DFT REAL PART 

VARIANCE OF DFT REAL PART 

QUADRAfn SNROUT OF DFT REAL PART 



10 

-181 .83990 
197 . 03 

167 .8235 = 22.2A85 DB 



MEAN OF DFT IMAG PART = -191.60988 
VARIAfiCE OF DFT IMAG PART = 30A.AA 

QUADRAfiT SNROUT OF DFT IMAG PART = 120.5977 = 20.813A DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE ATH QUADRANT xx 





NUMBER 


OF 


TOriES 


= 


9 






MEAN 


OF 


DFT 


REAL 


PART 


= 


178.21101 






VARIANCE 


OF 


DFT 


REAL 


PART 


= 


A05. 12 






QUADRANT SNROUT 


OF 


DFT 


REAL 


PART 


= 


78.3938 = 


18.9A28 


DB 


MEAf^ 


OF 


DFT 


IMAG 


PART 


r 


-18A .21101 






VARIANCE 


OF 


DFT 


IMAG 


PART 


= 


211 .55 






QUADRAfJT SNROUT 


OF 


DFT 


IMAG 


PART 


r 


160.A063 = 


22.0522 


DB 



XXX OVERALL (REAL + IMAG) STATISTICS FOR BAUD TYPE 2 xxx 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIAfiCE 
BAUD SrmOUT 



6A 

181.85138 ’ 

282.75 

116.9591 = 20.6803 DB 
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INPUT SMRNB = 100.000 = 20.000 DB 

BAUD TYPE 3: KX = 102A SAMPLE POINTS; K = 6A TONES 



XK GIVEN THE TRANSMITTED PHASE IS IN THE 1ST QUADRANT XX 



NUMBER OF TONES 
MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 

QUADRANT SNROUT OF DFT REAL PART 



13 

355.62256 

1701.75 

7A.3160 = 18.7108 DB 



MEAN OF DFT IMAG PART 
VARIANCE OF DFT IMAG PART 
QUADRANT SNROUT' OF DFT IMAG PART 



353 .6A551 
7A9 . lA 

166.9AA5 = 22.2257 DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT XX 





NUMBER 


OF ' 


TONES 


= 


15 






MEAN 


OF 


DFT 


REAL 


PART 


= 


"376 .6860A 






VARIANCE 


OF 


DFT 


REAL 


PART 




107A .68 






QUADRANT SNROUT 


OF 


DFT 


REAL 


PART 


— 


132.0519 = 


21.2068 


DB 


MEAN 


OF 


DFT 


IMAG 


PART 


r 


365.62622 






VARIANCE 


OF 


DFT 


IMAG 


PART 


= 


922.01 






QUADRANT SNROUT 


OF 


DFT 


IMAG 


PART 


= 


199.9395 = 


21.6139 


DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 3RD QUADRANT XX 





NUMBER 


OF ‘ 


TONES 


= 


20 






MEAN 


OF 


DFT 


REAL 


PART 


- 


-362.62378 






VARIANCE 


OF 


DFT 


REAL 


PART 


- 


1577.91 






QUADRANT SNROUT 


OF 


DFT 


REAL 


PART 


= 


83.3359 = 


19.2083 


DB 


MEAN 


OF 


DFT 


IMAG 


PART 


3 


-355.27930 






VARIANCE 


OF 


DFT 


IMAG 


PART 


3 


1101 .76 






QUADRANT SNROUT 


OF 


DFT 


IMAG 


PART 


3 


119.5697 = 


20.5905 


DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE ATH QUADRANT XX 



NUMBER OF TONES 
MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 

QUADRANT SNROUT OF DFT REAL PART 



16 

359 .399A1 
977 .57 

132.1315 = 21.2101 DB 



MEAN OF DFT IMAG PART 

VARIANCE OF DFT IMAG PART 

QUADRANT SNROUT OF DFT IMAG PART 



-376 . A3091 
12A8 . A5 

113.5008 = 20.5500 DB 



XXX OVERALL (REAL + IMAG) STATISTICS FOR BAUD TYPE 3 XXX 



TOTAL NUMBER OF POIfJTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

363.17578 • 

1115.50 

118.2399 = 20.7276 DB 
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INPUT SNRNB = 100.000 = 20.000 DB 

BAUD TYPE KX = 20A8 SAMPLE POINTS; K = 128 TONES 



GIVEN THE TRANSMITTED PHASE IS IN THE 1ST QUADRANT xx 





NUMBER 


OF 


TONES 


= 


38 






MEAN 


OF 


DFT 


REAL 


PART 


- 


792.11963 






VARIANCE 


OF 


DFT 


REAL 


PART 


= 


9625.69 






QUADRANT SNROUT 


OF 


DFT 


REAL 


PART 


= 


119.0627 = 


20.7578 


DB 


MEAN 


OF 


DFT 


IMAG 


PART 




715.70630 






VARIANCE 


OF 


DFT 


IMAG 


PART 


= 


9957.19 






QUADRANT SNROUT 


OF 


DFT 


IMAG 


part 


= 


119.9235 = 


20 .60^il 


DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT xx 

' NUMBER OF TONES = 3^ 

MEAN OF DFT REAL PART = -760.3691^ 

VARIAfJCE OF DFT REAL PART = 5032.17 

QUADRANT SNROUT OF DFT REAL PART = 11A.8930 = 20.6029 DB 

MEAN OF DFT IMAG PART = 755.90723 

VARIANCE OF DFT IMAG PART = 3920.06 

QUADRANT SNROUT OF DFT IMAG PART = 1A5.7620 = 21.636A DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 3RD QUADRANT XX 



NUMBER OF TONES 
MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 

QUADRANT SNROUT OF DFT REAL PART 



35 

>738.78^18 

5230.00 

10^4.3598 = 20.1853 DB 



MEAfJ OF DFT IMAG PART 
VARIANCE OF DFT IMAG PART 
QUADRANT SNROUT OF DFT IMAG PART 



-732.57861 
9859 . 96 

110.5907 = 20.9352 DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 9TH QUADRANT XX 





NUMBER 


OF 


TONES 


= 


21 






MEAN 


OF 


DFT 


REAL 


PART 


= 


750.57007 






VARIANCE 


OF 


DFT 


REAL 


PART 


= 


9882.30 






QUADRANT SrmOUT 


OF 


DFT 


REAL 


PART 




115.3872 = 


20.6216 


DB 


MEAtl 


OF 


DFT 


IMAG 


PART 


= 


-726 . 06567 






VARIANCE 


OF 


DFT 


IMAG 


PART 


= 


9869.59 






QUADRANT SNROUT 


OF 


DFT 


IMAG 


PART 


= 


108.3701 = 


20.3991 


DB 



XXX OVERALL (REAL + IMAG) STATISTICS FOR BAUD TYPE 9 XXX 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAf) 
BAUD VARIAfiCE 
BAUD SrmOUT 



256 

790.06909 > 

9585.59 

119.9913 = 20.7715 DB 
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INPUT SNRMD = 100.000 = 20.000 DD 

BAUD TYPE 5; KX = A096 SAMPLE POINTS; K = 256 TOMES 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 1ST QUADRANT xx 



MEAN 
VARIANCE 
QUADRANT SNROUT 



NUMBER OF TONES 
OF DFT REAL PART 
OF DFT REAL PART 
OF DFT REAL PART 



6 ^ 

1508.59062 
31A07 .17 

72.AA3A = 18.6000 DB 



MEAN OF DFT IMAG PART 

VARIANCE OF DFT IMAG PART 

QUADRANT SNROUT OF DFT IMAG PART 



1531.78125 

29779.73 

78.7905 = 18.96A7 DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT XX 





NUMBER 


OF • 


TONES 


= 


61 






MEAN 


OF 


DFT 


REAL 


PART 


= 


-1530.73755 






VARIANCE 


OF 


DFT 


REAL 


PART 


= 


25A65.21 






QUADRANT SNROUT 


OF 


DFT 


REAL 


PART 


— 


92.01A0 = 


19.6385 


DB 


MEAN 


OF 


DFT 


IMAG 


PART 


- 


1A89.81AA5 






VARIANCE 


OF 


DFT 


IMAG 


PART 


- 


3355A .58 






QUADRANT SNROUT 


OF 


DFT 


IMAG 


PART 


- 


66.1868 = 


18.2077 


DB 



XX GIVEN THE TRANSMITTED PHASE IS IN THE 3RD QUADRANT XX 





NUMBER 


OF ■ 


TONES 


= 


59 






MEAfJ 


OF 


DFT 


REAL 


PART 


= 


-1500.62695 






VARIAfJCE 


OF 


DFT 


REAL 


PART 


= 


17125.95 






QUADRANT Sf^OUT 


OF 


DFT 


REAL 


PART 




131.9895 = 


21 .1889 


DB 


MEAfJ 


OF 


DFT 


IMAG 


PART 


= 


-1983.72376 






VARIAfJCE 


OF 


DFT 


IMAG 


PART 


= 


26270.10 






QUADRANT SNROUT 


OF 


DFT 


IMAG 


PART 


- 


83.8006 = 


19.2325 


DB 



XX GIVEN THE TRAfJSMITTED PHASE IS IN THE ATH QUADRANT xx 



NUMBER OF TONES 
MEAfJ OF DFT REAL PART 

VARIAtICE OF DFT REAL PART 

QUADRANT SNROUT OF DFT REAL PART 



72 

1A93. 180A2 
2208A .69 

100.9562 = 20.0A13 DB 



MEAN OF DFT IMAG PART 

VARIANCE OF DFT IMAG PART 

QUADRANT SNROUT OF DFT IMAG PART 



-1510.180^42 
29661 .86 

76.8881 = 18.8586 DB 



XXX OVERALL (REAL + IMAG) STATISTICS FOR BAUD TYPE 5 xxx 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIAf]CE 
BAUD SNROUT 



512 

1506.13916 

26536.55 

85.323A = 19.5107 DB 
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APPENDIX D. STATISTICS OF SNRout VS. DOPPLER ANALYSIS 



A. INPUT SNR = 15 dB 



SNROUT VS. DOPPLER TEST -- EPSILON = 0.00 

INPUT SNRNB = 31.623 = 15.000 DB 

BAUD TYPE 3: KX = 102A SAMPLE POINTS; K = 6A TONES 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 1ST BAUD xx 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

358.19678 

361A.89 

35,A935 = 15.5015 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 2ND BAUD XX 



TOTAL NUMBER OF POULTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

359 . 38672 
A296 .25 

30.0631 = 1A.7803 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 3RD BAUD XX 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIAtiCE 
BAUD SNROUT 



128 

353.96387 

5632.05 

3A.‘^953 = 15.3777 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE ATH BAUD XX 



TOTAL NUMBER OF POIliTS, 2K 
BAUD MEAti 
BAUD VARIANCE 
BAUD SNROUT 



128 

360 . 0A858 
3868 . OA 

33.51AA = 15.2523 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 5TH BAUD xx 



TOTAL NUMBER OF POlUJS, 2K 
BAUD MEAJ] 
BAUD VARIANCE 
BAUD SNROUT 



128 

3^^ .^^60A 
520A . A3 

22.7965 = 13.5787 DB 



XXX TOTAL OVER ALL 5 BAUDS WITH EPSILON = 0.00 XXX 



TOTAL NUMBER OF BAUDS 
MEAfJ SNROUT 
STANDARD DEVIATIOfJ OF SNROUT 



5 

1A.8981 DB 
0.7865 DB 
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SNROUT VS. DOPPLER TEST -- EPSILON = 0.25 

INPUT SNRNB = 31.623 = 15.000 DB 

BAUD TYPE 3: KX = 102A SAMPLE POINTS; K = 6A TONES 



OVERALL (REAL + IMAG) STATISTICS FOR THE 1ST BAUD xx 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

356.60303 

A178.20 

30.A355 = 1A.8338 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 2ND BAUD xx 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIAfJCE 
BAUD SNROUT 



128 

357 .837^^0 
A579 . 19 

27.9630 = 1A.A658 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 3RD BAUD xx 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

352.0A83A 

3697.63 

33.5182 = 15.2528 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE ATH BAUD xx 



TOTAL NUMBER OF POItiTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

358 . 16A06 
A301 .95 

29.819A = 1A.7A50 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 5TH BAUD XX 



TOTAL NUMBER OF, POItiTS, 2K 
BAUD MEAN 
BAUD VARIAtICE 
BAUD StmOUT 



128 

3A2. 3A058 
590A.97 

19.8A72 = 12.9770 DB 



XXX TOTAL OVER ALL 5 BAUDS WITH EPSILON = 0.25 XXX 



TOTAL NUMBER OF BAUDS 
MEAN SriROUT 
STANDARD DEVIATION OF StiROUT 



5 

1^.^5^^9 DB 
0.8730 DB 



113 



SNROUT VS. DOPPLER TEST -- EPSILON = 0.50 

INPUT SNRNB = 31.623 = 15.000 DB 

BAUD TYPE 3: KX = 102A SAMPLE POINTS; K = 6A TONES 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 1ST BAUD xx 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

351.68213 
5379. 5A 

22.9909 = 13.6156 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 2ND BAUD XX 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

353.07666 

56A5.37 

22.0823 = 13.‘4‘40^ DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 3RD BAUD xx 



TOTAL NUMBER OF POIfHS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

3^6.26807 

A283.9A 

27.9886 = 1A.A698 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE ‘^TH BAUD xx 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIAriCE 
BAUD SNROUT 



123 

352.39697 

5A59.51 

22.6220 = 13.5A53 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 5TH BAUD XX 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

335.90625 
720A .71 

15.6610 = 11 .9A82 DB 



XXX TOTAL OVER ALL 5 BAUDS WITH EPSILOfi = 0.50 ,xxx 



TOTAL NUMBER OF BAUDS 
MEAN SNROUT 
STANDARD DEVIATION OF St^ROUT 



5 

13.A039 DB 
0.9112 DB 
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SNROUT VS. DOPPLER TEST — EPSILON = 0.75 

INPUT SNRNB = 31.623 = 15.000 DB 

BAUD TYPE 3: KX = 102A SAMPLE POINTS; K = 6A TONES 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 1ST BAUD XX 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

3A3.56AA5 

7197.11 

16 .A005 = 12.1A86 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 2ND BAUD xx 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

3^5.20166 

7A79.9^ 

15.9312 = 12.0225 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 3RD BAUD XX 



TOTAL NUMBER OF POIfJTS, 2K 
BAUD MEAN 
BAUD VARIAfJCE 
BAUD SNROUT 



128 

336 .72A85 
5392.5^ 

21.0260 = 13.2276 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE ATH BAUD XX 



TOTAL NUMBER OF POIfJTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

3A2. 88135 
7A09 .56 

15.8670 = 12.00A9 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 5TH BAUD xx 



TOTAL NUMBER OF POItJTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

325.28076 

9073.28 

11 .661A = 10.6675 DB 



XXX TOTAL OVER ALL' 5 BAUDS WITH EPSILOfi = 0.75 xxx 



TOTAL NUMBER OF BAUDS 
MEAN SfJROUT 
STANDARD DEVIATION OF SNROUT 



5 

12.01A2 DB 
0.9088 DB 
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SNROUT VS. DOPPLER TEST -- EPSILON = 1.00 

INPUT SNRNB = 31.623 = 15.000 DB 

BAUD TYPE 3; KX = 102A SAMPLE POINTS; K = 6A TONES 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 1ST BAUD xx 



TOTAL NUMBER OF POINTS; 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

332. A69A8 
9581 .lA 

11.5368 = 10.6209 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 2ND BAUD XX 



TOTAL NUMBER OF POINTS; 2K 
BAUD MEAN 
BAUD VARIAfiCE 
BAUD SNROUT 



128 

33A . A6362 
10027 . 16 

11.1563 = 10.A752 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 3RD BAUD xx 



TOTAL NUMBER OF POINTS; 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

323.73706 

7001.29 

1A.9695 = 11.7521 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE ATH BAUD XX 



TOTAL NUMBER OF POItiTS; 2K 
BAUD MEAN 
BAUD VARIAflCE 
BAUD Sf^ROUT 



128 

329.95630 
10013. A6 

10.8725 = 10.3633 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 5TH BAUD xx 



TOTAL NUMBER OF POINTS; 2K 
BAUD MEAf] 
BAUD VARIAflCE 
BAUD Sr^ROUT 



128 

310.86A50 
11A50 . 02 

8.A399 = 9.263A DB 



XXX TOTAL OVER ALL 5 BAUDS MITH EPSILON = 1.00 xxx 



TOTAL NUMBER OF BAUDS 
MEAti SfiROUT 
STANDARD DEVIATIOri OF SrmOUT 



5 

10.A9A9 DB 
0.88A7 DB 
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SNROUT VS. DOPPLER TEST -- EPSILON = 1.25 

INPUT SNRNB = 31.625 = 15.000 DB 

BAUD TYPE 3: KX = 102A SAMPLE POINTS; K = 6A TONES 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 1ST BAUD XX 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

518.5879A 
125A8 .15 

8.0786 = 9.0755 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 2ND BAUD xx 



TOTAL NUMBER OF POItJTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

520.88208 

15305.20 

7.7587 = 8.8867 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 3RD BAUD xx 



TOTAL NUMBER OF POIfJTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

307 .26807 
9126.79 

10.59^7 = 10.1A72 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE ATH BAUD XX 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

515.63057 

15515.15 

7 .3879 = 8.6899 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 5TH BAUD XX 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

292.66162 

19553.81 

5.9671 = 7.7576 DB 



XXX TOTAL OVER ALL 5 BAUDS HITH EPSILON = 1.25 xxx 



TOTAL NUMBER OF BAUDS 
MEAN SNROUT 
STANDARD DEVIATION OF SNROUT 



5 

8.9099 DB 
0.8568 DB 
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SNROUT VS. DOPPLER TEST -- EPSILON = 1.50 

INPUT SNRNB = 31.623 = 15.000 DB 

BAUD TYPE 3: KX = 102^ SAMPLE POINTS; K = 6A TONES 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 1ST BAUD xx 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD StJROUT 



128 

301 .6A111 
16039.55 

5.6727 = 7.5379 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 2ND BAUD xx 



TOTAL NUMBER OF POItJTS, 2K 
BAUD MEAN 
BAUD VARIAt^CE 
BAUD SfiROUT 



128 

30A. 77563 
172A2.^9 

5.3872 = 7 .3136 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 3RD BAUD xx 



TOTAL fJUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SflROUT 



128 

287 .72A12 
11751.02 

7 . 0AA9 = 8 . A788 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE ATH BAUD XX 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAtl 
BAUD VARIAfICE 
BAUD SfiROUT 



128 

29A. 33276 
17252.69 

5.0213 = 7.0082 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 5TH BAUD xx 



TOTAL NUMBER OF POItiTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

271.16650 
17716 .52 

A.150A = 6.1809 DB 



XXX TOTAL OVER ALL 5 BAUDS NITH EPSILON = 1.50 xxx 



TOTAL fJUMBER OF BAUDS 
MEAN SNROUT 
STANDARD DEVIATIOfi OF SNROUT 



5 

7 .3039 DB 
0.S3A2 DB 
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SNROUT VS. DOPPLER TEST -- EPSILON = 1.75 

INPUT SNRNB = 31.623 = 15.000 DB 

BAUD TYPE 3: KX = 102A SAMPLE POINTS; K = 6A TONES 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 1ST BAUD xx 



TOTAL NUMBER OF POINTS. 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

282.71875 

19969.26 

A. 0026 = 6 .0235 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 2ND BAUD XX 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

236 .6A087 
21735.70 

3.7801 = 5.7750 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 3RD BAUD xx 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

265.68335 
1A806 . A8 

A .7673 = 6 .7828 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE ATH BAUD XX 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIAfICE 
BAUD SNROUT 



128 

272.69653 

21720.69 

3.A236 = 5.3AA9 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 5TH BAUD XX 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

2A7 . 10A03 
21A51 .55 

2.8A6A = A.5A30 DB 



XXX TOTAL OVER ALL 5 BAUDS WITH EPSILON = 1.75 



TOTAL NUMBER OF BAUDS 
MEAN SNROUT 
STANDARD DEVIATION OF SNROUT 



5 

5.6938 DB 
0.8287 DB 
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SNROUT VS. DOPPLER TEST — EPSILON = 2.00 

INPUT SNRNB = 31.625 = 15.000 DB 

BAUD TYPE 5: KX = 102^ SAMPLE POINTS; K = 6A TONES 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 1ST BAUD XX 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

261 . 53735 
2A362.90 

2.8076 = A.A83A DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 2ND BAUD XX 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

266 . A6338 
2681A .79 

2.6A79 = 9.2290 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 3RD BAUD xx 



TOTAL NUMBER OF POiriTS, 2K 
BAUD MEAli 
BAUD VARIANCE 
BAUD SNROUT 



128 

292.98599 

18393.12 

3.2188 = 5.0769 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 9TH BAUD XX 



TOTAL NUMBER OF POINTS, 2K 
BAUD 

BAUD VARIANCE 
BAUD StiROUT 



128 

298.73959 

26755.95 

2.3129 = 3.6906 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 5TH BAUD xx 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIAtiCE 
BAUD SNROUT 



128 

239.01096 

25535.32 

2.2328 = 3.9889 DB 



XXX TOTAL OVER ALL 5 BAUDS NITH EPSILON = 2.00 xxx 



TOTAL NUMBER OF BAUDS 
MEAN StIROUT 
STANDARD DEVIATION OF SNROUT 



5 

9.1837 DB 
0.6958 DB 
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SNROUT VS. DOPPLER TEST -- EPSILON = 2.25 

INPUT SNRNB = 31.623 = 15.000 DB 

BAUD TYPE 3: KX = 102A SAMPLE POINTS; K = 6A TONES 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 1ST BAUD XX 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

2A3.6AA71 

29136.31 

2.037A = 3.0908 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 2ND BAUD xx 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SfiROUT 



128 

2A8 .93130 
32359.91 

1 .91A9 = 2.8215 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 3RD BAUD XX 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

2A0, 99722 
22283.86 

2.606A = A. 1603 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE ATH BAUD XX 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD StJROUT 



128 

230,72388 

32258.76 

1.6512 = 2,1781 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 5TH BAUD xx 



TOTAL fiUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

2AA . 09976 
300A1 .91 

1 .983A = 2.97A1 DB 



XXX TOTAL OVER ALL 5 BAUDS NITH EPSILON = 2.25 XXX 



TOTAL NUMBER OF BAUDS 
MEAtJ SNROUT 
STANDARD DEVIATIOti OF SliROUT 



5 

3.0A50 DB 
0.7163 DB 
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5NR0UT VS. DOPPLER TEST — EPSILON = 2.50 

INPUT SNRNB = 31,. 623 = 15.000 DB 

BAUD TYPE 3: KX = 102A SAMPLE POINTS; K = 6A TOtJES 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 1ST BAUD xx 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

2A0 . 0A703 
3A182.69 

1.6857 = 2.2679 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 2ND BAUD XX 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAf^ 
BAUD VARIAt^CE 
BAUD SNROUT 



128 

23A .67^13 
382A9 .63 

1.A398 = 1.5830 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 3RD BAUD XX 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD StmOUT 



128 

25A .86A0A 
26569.81 

2.9952 = 3.8831 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 9TH BAUD xx 



TOTAL NUMBER OF POIfiTS, 2K 
BAUD MEAt^ 
BAUD VARIAfICE 
BAUD SNROUT 



128 

228.89236 

38099.78 

1.3771 = 1.3897 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 5TH BAUD XX 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SfmOUT 



128 

255.39229 

39706.96 

1.8793 = 2.7900 DB 



XXX TOTAL OVER ALL 5 BAUDS WITH EPSILON = 2.50 xxx 



TOTAL NUMBER OF BAUDS 
MEAf^ SfmOUT 
STAfJDARD DEVIATIOt^ OF SNROUT 



5 

2.3727 DB 
1.0021 DB 
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SNROUT VS* DOPPLER TEST — EPSILON = 2.75 

INPUT SNRNB = 31.623 = 15.000 DB 

BAUD TYPE 3: KX = 102^ SAMPLE POINTS; K = 6^ TONES 



OVERALL (REAL + IMAG) STATISTICS FOR THE 1ST BAUD 



TOTAL NUMBER OF POIfJTS, 2K 
BAUD MEAN 
BAUD VARIAfiCE 
BAUD SfiROUT 



128 

250.07709 
39538 . 09 

1 .5817 = 1.9913 DB 



OVERALL (REAL + IMAG) STATISTICS FOR THE 2ND BAUD XX 



TOTAL [DUMBER OF POItJTS, 2K 
BAUD MEAf^ 
BAUD VARIAIiCE 
BAUD SrmOUT 



128 

225.95108 

^^ 515.95 

1.1^69 = 0.5951 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 3RD BAUD XX 



TOTAL NUMBER OF POIf^TS, 2K 
BAUD MEAN 
BAUD VARIAfiCE 
BAUD Sf^ROUT 



128 

265.656A9 

31226.82 

2.2600 = 3.5^11 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE ^TH BAUD xx 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SfIROUT 



128 

236.691^7 
^^ 4222. 02 

1 .2669 = 1 . 0273 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 5TH BAUD xx 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

262.792^8 
39631 .13 

1.7^26 = 2.A119 DB 



XXX TOTAL OVER ALL 5 BAUDS WITH EPSILOfi = 2.75 xxx 



TOTAL [iUMBER OF BAUDS 
MEAN SNROUT 
STANDARD DEVIATIOfi OF SNROUT 



5 

1 .9133 DB 
1 . 16‘^8 DB 
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B. INPUT SNR = 5 dB 



SNROUT VS. DOPPLER TEST -- EPSILON = 0.00 

INPUT SNRNB = 3.162 = 5.000 DB 

BAUD TYPE 3: KX = 1024 SAMPLE POINTS; K = 64 TONES 



XX OVERALL (REAL + IMAG) STATISTICS* FOR THE 1ST BAUD XX 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

348.85937 

36188.43 

3.3630 = 5.2673 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 2ND BAUD xx 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SflROUT 



128 

352.66479 
43009 .19 

2.8918 = 4.6116 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 3RD BAUD xx 



TOTAL NUMBER OF POIfiTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SfiROUT 



128 

335 . 33081 
36318.57 

3.0961 = 4.9082 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 4TH BAUD xx 



TOTAL NUMBER OF POIfiTS, 2K 
BAUD MEAfJ 
BAUD VARIANCE 
BAUD SfiROUT 



128 

354.62695 

33677.65 

3.2515 = 5.1208 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 5TH BAUD XX 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

305.20776 
52039 .11 

1.7900 = 2.5286 DB 



XXX TOTAL OVER ALL 5 BAUDS HITH EPSILON = 0.00 xxx 



TOTAL NUMBER OF BAUDS 
MEAN SfiROUT 
STAfJDARD DEVIATION OF SfJROUT 



5 

4.4873 DB 
1 . 1224 DB 
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SNROUT VS. DOPPLER TEST — EPSILON = 0.25 

INPUT SNRNB = 3.162 = 5.000 DB 

BAUD TYPE 3: KX = 102A SAMPLE POINTS; K = 6A TONES 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 1ST BAUD XX 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

3A7 .25366 
37339.82 

3.229A = 5.0912 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 2ND BAUD xx 



TOTAL NUMBER OF POIt^TS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SriROUT 



128 

351 ,10767 
A31A3 , 09 

2.857A = A. 5597 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 3RD BAUD XX 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SfmOUT 



128 

333. A1895 
35979 .6A 

3.0898 = A. 8992 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE ATH BAUD xx 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIAf^CE 
BAUD SNROUT 



128 

352.7AA38 
392A6 . lA 

3.1705 = 5.0112 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 5TH BAUD xx 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIAf^CE 
BAUD SNROUT 



128 

303. 096AA 
5355A. 99 

1 .7154 = 2.3436 DB 



XXX TOTAL OVER ALL 5 BAUDS WITH EPSILON = 0.25 px 



TOTAL NUMBER OF BAUDS 
MEAN Sr^ROUT 
STANDARD DEVIATIOf^ OF SNROUT 



5 

4.3810 DB 
1.1568 DB 
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SNROUT VS. DOPPLER TEST — EPSILON = 0.50 

INPUT Sf^RNB = 3.162 = 5.000 DB 

BAUD TYPE 3: KX = 102A SAMPLE POINTS; K = 6A TONES 



OVERALL (REAL + IMAG) STATISTICS FOR THE 1ST BAUD 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIAfJCE 
BAUD SrmOUT 



128 

3 ^ 2 . 3^4521 

39210.31 

2.9890 = A. 7553 DB 



OVERALL (REAL + IMAG) STATISTICS FOR THE 2ND BAUD XK 



TOTAL NUMBER OF POIfiTS, 2K 
BAUD MEAfJ 
BAUD VARIANCE 
BAUD SNROUT 



128 

3A6 .3AAA8 
AA183.36 

2.71A9 = ^4.3376 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 3RD BAUD XX 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

327.63818 

36156.82 

2.9689 = A. 7260 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE ^TH BAUD xx 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

3A6 .979A9 
A0537 . A1 

2.9700 = ^.7275 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 5TH BAUD xx 



TOTAL NUMBER OF POIfJTS, 2K 
BAUD MEAN 
BAUD VARIAf^CE 
BAUD srmouT 



128 

296.65332 
55526 .02 

1.58A9 = 2.0000 DB 



XXX TOTAL OVER ALL 5 BAUDS MITH EPSILON = 0.50 xxx 



TOTAL NUMBER OF BAUDS 
MEAN SNROUT 
STANDARD DEVIATION OF SNROUT 



5 

A. 1093 DB 
1 .1917 DB 
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5NR0UT VS. DOPPLER TEST — • EPSILON = 0.75 

INPUT SNRNB = 3.162 = 5.000 DB 

BAUD TYPE 3: KX = 102A SAMPLE POINTS; K = 6A TONES 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 1ST BAUD X5( 



TOTAL NUMBER OF POINTS. 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

33A. 21777 
A1755.36 

2.6751 = A. 2735 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 2ND BAUD xx 



TOTAL NUMBER OF POINTS. 2K 
BAUD MEAN 
BAUD VARIAr^CE 
BAUD SNROUT 



128 

338 . A6997 
A6133.0A 

2.A833 = 3.9503 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 3RD BAUD XX 



TOTAL NUMBER OF POIf^S. 2K 
BAUD MEAfJ 
BAUD VARIANCE 
BAUD SNROUT 



128 

318 . 0981A 
36878.21 

2.7A38 = A. 3835 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE ATH BAUD XX 



TOTAL NUMBER OF POIfiTS. 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

337 . A6533 
A2530 . A1 

2.6777 = A. 2776 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 5TH BAUD xx 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAt^ 
BAUD VARIANCE 
BAUD SNROUT 



128 

286 . 0A077 
57909.20 

1.A129 = 1.5011 DB 



XXX TOTAL OVER ALL 5 BAUDS NITH EPSILON = 0.75 XXX 



TOTAL NUMBER OF BAUDS 
MEAf^ SNROUT 
STANDARD DEVIATIOf^ OF SUROUl 



5 

Z.677Z DB 
1.2273 DB 
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SNROUT VS. DOPPLER TEST — EPSILON = 1.00 

INPUT SURUh = 3.162 = 5.000 DB 

BAUD TYPE 3: KX = 102^ SAMPLE POINTS; K = 6^ TOfiES 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 1ST BAUD xx 



TOTAL NUMBER OF POINTS. 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

323.12573 
<^^ 902. 68 

2.3253 = 3.66^7 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 2ND BAUD xx 



TOTAL NUMBER OF POINTS. 2K 
BAUD MEAN 
BAUD VARIAtiCE 
BAUD StmOUT 



128 

327 .73608 
^890A .0^ 

2.196^ = 3.A170 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 3RD BAUD XX 



TOTAL NUMBER OF POINTS. 2K 
BAUD MEAfJ 
BAUD VARIAf^CE 
BAUD SfiROUT 



128 

305 . 10A98 
38135.35 

2.^‘^10 = 3.8757 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE ^TH BAUD xx 



TOTAL fJUMBER OF POINTS. 2K 
BAUD MEAri 
BAUD VARIAtiCE 
BAUD StJROUT 



128 

32^.53516 

^5177.82 

2.3313 = 3.6760 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 5TH BAUD xx 



TOTAL fJUMBER OF POItiTS. 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

271 .62^76 
60591 .^7 

1 .2177 = 0.8553 DB 



XXX TOTAL OVER ALL 5 BAUDS WITH EPSILOt^ = 1.00 xxx 



TOTAL NUMBER OF BAUDS 
MEAN SNROUT 
STANDARD DEVIATIOf^ OF SNROUT 



5 

3.0977 DB 
1 .26A1 DB 
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SNROUT VS. DOPPLER TEST — EPSILON = 1.25 

INPUT SNRNB = 3.162 = 5.000 DB 

BAUD TYPE 3: KX = 102A SAMPLE POINTS; K = 6A TONES 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 1ST BAUD XX 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAf^ 
BAUD VARIANCE 
BAUD SNROUT 



128 

309 . 05103 
A866A.3A 

1 .9627 = 2.9285 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 2ND BAUD XX 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

31A. 15283 
5252A.69 

1 .8790 = 2.7392 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 3RD BAUD xx 



TOTAL NUMBER OF POirnS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

288.63892 
39937 .87 

2.0861 = 3.1932 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE ATH BAUD xx 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAfJ 
BAUD VARIANCE 
BAUD SNROUT 



128 

308 . 20581 
A8AS2.35 

1.9593 = 2.9210 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 5TH BAUD XX 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIAtJCE 
BAUD srmouT 



128 

253. A2087 
63620.02 

1.0095 = 0.0A09 DB 



XXX TOTAL OVER ALL 5 BAUDS NITH EPSILON = 1.25 XXX 



TOTAL NUMBER OF BAUDS 
MEArj SNROUT 
STANDARD DEVIATION OF SNROUT 



5 

2.36A6 DB 
1.3090 DB 
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SNROUT VS. DOPPLER TEST • EPSILON = 1.50 

INPUT SNRNB = 5.162 = 5.000 DB 

BAUD TYPE 3: KX = 102A SAMPLE POINTS; K = 6A TONES 



OVERALL (REAL + IMAG) STATISTICS FOR THE 1ST BAUD 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

292.29761 

52952.92 

1.6155 = 2.0776 DB 



OVERALL (REAL + IMAG) STATISTICS FOR THE 2ND BAUD 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

298.0^665 

56898.22 

1.5612 = 1.95A7 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 3RD BAUD X)( 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIAfICE 
BAUD SfmOUT 



128 

269 . 08667 
92309.^0 

1.7116 = 2.33A0 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE ATH BAUD X)( 



TOTAL NUMBER OF POIfITS, 2K 
BAUD MEAtJ 
BAUD VARIAfJCE 
BAUD SNROUT 



128 

288.91211 
52A05. 19 

1.5928 = 2.0216 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 5TH BAUD xx 



TOTAL NUMBER OF POirJTS, 2 K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SfJROUT 



128 

237.78751 

66898.56 

0.8^52 = -0.730A DB 



XXX TOTAL OVER ALL 5 BAUDS WITH EPSILON = 1.50 xxx 



TOTAL NUMBER OF BAUDS 
MEAfJ SNROUT 
STANDARD DEVIATION OF StJROUT 



5 

1 .5275 DB 
1.2709 DB 
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SNROUT VS. DOPPLER TEST — EPSILON = 1.75 

INPUT SNRNB = 3.162 = 5.000 DB 

BAUD TYPE 3: KX = 102A SAMPLE POINTS; K = 6A TONES 



OVERALL (REAL + IMAG) STATISTICS FOR THE 1ST BAUD XX 



TOTAL NUMBER OF POINTS^ 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

278.15869 

57652.27 

1.3A20 = 1.2777 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 2ND BAUD XX 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

288.96216 

61909.75 

1.3A87 = 1.2992 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 3RD BAUD XX 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SfJROUT 



128 

2A7 . 0A657 
^5172.92 

1.3511 = 1.3068 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE ATH BAUD XX 



TOTAL NUMBER OF POINTS, 2K 
' BAUD MEAN 
BAUD VARIArjCE 
BAUD SNROUT 
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267 .27612 
56809 . 3A 

1.2575 = 0.9950 DB 



XX OVERALL (REAL + IMAG) STATISTICS FOR THE 5TH BAUD XX 



TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 



128 

2A1 .3063A 
703A5.50 

0.8278 = -0.8210 DB 



XXX TOTAL OVER ALL 5 BAUDS WITH EPSILON = 1.75 ^xxx 



TOTAL NUMBER OF BAUDS 
MEAN SNROUT 
STANDARD DEVIATION OF SNROUT 



5 

0.8115 DB 
0.9218 DB 
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APPENDIX E. STATISTICS OF THE DOPPLER ESTIMATION 
ALGORITHM 



A- INPUT SNR = 10 dB 



ALPHA-HAT VS. ALPHA — INPUT SNR = 10 DB 
BAUD TYPE 3: KX = 102A SAMPLE POINTS; 



K = 6A TONES 



EPSILON =0.0 


DELTA-ALPHA 


= 0.000376506 


LL 


ALPHA 


ALPHA-HAT 


EPSILON-HAT 


1 


0.000 


0.00008162 


0.21678060 


2 


0.000 


-0 .00005518 


-0.1A65710A 


5 


0.000 


0 . 00000018 


0 . 000A6881 




0 .000 


-0 . 0002912A 


-0.77353AA8 


5 


0.000 


-0 .00000678 


-0 . 01801626 





NUMBER OF 


BAUDS = 


5 




MEAN OF EPSILON-HAT = -0.1AA17AA57 


VARIANCE OF EPSILON-HAT = 0.1A0806A37 


XX EPSILON = 0.25 


DELTA-ALPHA 


= 0.000376506 


LL 


ALPHA 


ALPHA-HAT 


EPSILON-HAT 


1 


0 . 00009A128 


0.00017936 


0 . 97636899 


2 


0 . 00009A128 


0 .00003786 


0 . 1 0056025 


3 


0 . 00009A123 


0.00011023 


0.29278039 


A 


0 . 000099128 


-0 .00020189 


-0.53603972 


5 


0.000099128 


0.00009992 


0.25209653 




NUMBER OF 


BAUDS = 


5 




MEAN OF EPSILON-HAT = 0.1171A315A 


VARIANCE OF EPSILON-HAT = 0.151272893 


XX EPSILON = 0.50 


DELTA-ALPHA 


= 0.000376506 


LL 


ALPHA 


ALPHA-HAT 


EPSILON-HAT 


1 


0 .000188255 


-0 .00011252 


-0 . 29886585 


2 


0 .000188255 


0 . 00013292 


0 . 35169991 


3 


0 . 000183255 


0 .00021961 


0.58327237 


A 


0 .000188255 


0.00011711 


0.31109386 


5 


0.000188255 


0 . 00019830 


0.52669251 




NUMBER OF 


BAUDS = 


5 




MEAN OF EPSILON-HAT = 0.299768512 



VARIANCE OF EPSILON-HAT = 0.125232782 
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ALPHA-HAT VS. ALPHA -- lUPVJ SNR = 10 DB 

BAUD TYPE 3: KX = 102A SAMPLE POINTS; K = 6A TONES 



XX EPSILON = 0.75 DELTA-ALPHA = 0.000376506 



XX 



LL 


ALPHA 


ALPHA-HAT 


EPSILON-HAT 


1 


0.000282383 


-0.0002^259 


-0.69918107 


2 


0.000282383 


0.00003668 


0.09791002 


3 


0.000282383 


0.00013307 


0.35399987 




0.000282383 


0.00021555 


0.57250100 


5 


0.000282383 


0.00030271 


0.80399001 




NUMBER OF 


BAUDS = 


5 




MEAN OF EPSIL 


ON-HAT = 0.236632889 


VARIANCE OF EPSILON-HAT = 0.510890436 



EPSILON =1.00 


DELTA-ALPHA 


= 0.000376506 


LL 


ALPHA 


ALPHA-HAT 


EPSILON-HAT 


1 


0 . 000376506 


-0 .00039569 


-0.91803193 


2 


0 .000376506 


0 .00013396 


0.35578859 


3 


0 . 000376506 


0 .00023758 


0.63100237 


9 


0 .000376506 


0 . 00031679 


0.89190372 


5 


0 .000376506 


0 .00019155 


0.50575187 



MEAN 

VARIANCE 



NUMBER OF BAUDS 
OF EPSILON-HAT 
OF EPSILON-HAT 



5 

0.283782721 
0 . A82913A9^ 



EPSILON =1.25 


DELTA-ALPHA 


= 0.000376506 


LL 


ALPHA 


ALPHA-HAT 


EPSILOfHHAT 


1 


0.000970638 


-0.00088195 


-2.39113503 


2 


0 .000970638 


0.00023501 


0.62918330 


3 


0 .000970638 


0.00039109 


0.90593851 


9 


0 .000970638 


0.00023151 


0.61989583 


5 


0 .000970638 


-0.00013302 


-0 . 35329860 



MEAN 

VARIANCE 



NUMBER OF BAUDS 
OF EPSILON-HAT 
OF EPSILOfJ-HAT 



5 

-0 .109883070 
1 .783^^536 
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ALPHA-HAT VS. ALPHA -- INPUT SNR = 10 DB 
BAUD TYPE 3: KX = 102A SAMPLE POINTS; 



K = 6A TONES 



XX EPSILON =1.50 



DELTA-ALPHA = 



0 .000576506 



LL 


ALPHA 


ALPHA-HAT _ 


EPSILON-HAT 


1 


0 . 000564765 


-0 .00100695 


-2.67AA6518 


2 


0.000564765 


-0.00051077 


-1 .35661030 


3 


0 . 000564765 


-0 . 00061805 


-1 .6A155102 


A 


0.000564765 


-0 . 00006197 


-0.16A59151 


5 


0.000564765 


-0 . 00002642 


-0.07016611 




NUMBER OF 


BAUDS = 


5 



MEAN OF EPSILOri-HAT = 
VARIANCE OF EPSILON-HAT = 



-1.181^7659 
1 .185105A2 
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B. INPUT SNR = 15 dB 



ALPHA-HAT VS. ALPHA -- INPUT SNR = 15 DB 
BAUD TYPE 3: KX = 102A SAMPLE POINTS; 



K = 6A TONES 



XX EPSILON = 0.0 



DELTA-ALPHA = 0.000576506 



LL 


ALPHA 


ALPHA-HAT 


EPSILON-HAT 


1 


0 . 000 


0 . 0000A005 


0.10637712 


2 


0 .000 


-0.00002999 


0 .0796A1A6 


3 


0.000 


-0.00000758 


0.02012977 


A 


0 . 000 


-0.00003A88 


0.0926A153 


5 


0.000 


0.00000226 


0.00600086 




rJUMBER OF 


BAUDS = 


5 




MEAN OF EPSILON-HAT = -0.16006 


9A66E-01 


VARIANCE OF EPSILOfi-HAT = 0 . 6 35 03A963E-02 


EPSILON = 0.25 


DELTA-ALPHA = 


0.000376506 


LL 


ALPHA 


ALPHA-HAT 


EPSILON-HAT 


1 


0 . 00009A128 


0 .00013718 


0 . 36A3A567 


2 


0 .000099128 


0.00006A37 


0.17096901 


3 


0.000099128 


0 . 00009A53 


0 .25108A21 


A 


0 .000099128 


0 .00005696 


0 .15127569 


5 


0 . 000099128 


0.00010090 


0.26797980 




NUMBER OF 


BAUDS = 


5 




MEAN OF EPSIL 


ON-HAT = 0.2A1130829 


VARIANCE OF EPSILON-HAT = 0 . 72A961256E-02 


EPSILON = 0.50 


DELTA-ALPHA = 


0 . 000376506 


LL 


ALPHA 


ALPHA-HAT 


EPSILON-HAT 


1 


0 .000188255 


0.00023A08 


0.62171990 


2 


0 .000188255 


0 .00016 005 


0 . «^2509770 


3 


0 .000188255 


0 .000197 05 


0.52335626 


A 


0 . 000188255 


0 .00015115 


0.A01‘^A7A2 


5 


0 . 000188255 


0 . 000200A3 


0.53233200 




rj UMBER OF 


BAUDS = 


5 



MEAti OF EPSILOf^-HAT = 0, 
VARIANCE OF EPSILON-HAT = 0. 



500790179 

793160371E-02 
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ALPHA-HAT VS. ALPHA — INPUT SNR = 15 DB 

BAUD TYPE 3: KX = 102^ SAMPLE POINTS; K = 6^ TONES 



XX EPSILON = 0.75 



DELTA-ALPHA = 0.000376506 



LL 


ALPHA 


ALPHA-HAT 


EPSILON-HAT 


1 


0.000282383 


0.00033001 


0.87650335 


2 


0 .000282383 


0 .00025643 


0.68107283 


3 


0.000282383 


0 .00029901 


0.79417539 




0.000282383 


0.00024693 


0.65585577 


5 


0 .000282383 


0.00029999 


0.79678440 




NUMBER OF 


BAUDS = 


5 




MEAN OF EPSILON-HAT = 0.760878146 


VARIANCE OF EPSILON-HAT = 0 . 829143077 E-02 


EPSI 


LON = 0.80 


DELTA-ALPHA 


= 0.000376506 


LL 


ALPHA 


ALPHA-HAT 


EPSILON-HAT 


1 


0 .000301205 


0.00034829 


0.92506492 


2 


0 .000301205 


0.00027502 


0.73044205 


3 


0 .000301205 


0.00031849 


0.84590751 


4 


0 .000301205 


0 .00026549 


0.70512855 


5 


0 .000301205 


0 .00031912 


0.84758091 




NUMBER OF 


BAUDS = 


5 




KEAH OF EPSILOII-HAT = 0.81082^39^ 


VARI 


AHCE OF EPSILON-HAT = 0 . 831692002E-02 



EPSI 


LON = 0.85 


DELTA-ALPHA 


= 0.000376506 


LL 


ALPHA 


ALPHA-HAT 


EPSILOfi-HAT 


1 


0 .000320030 


0.00036756 


0.97623181 


2 


0 .000320030 


0 .00029466 


0.78261632 


3 


0 .000320030 


0 .00033904 


0.90049165 


4 


0 .000320030 


0 .00028515 


0.75734794 


5 


0 .000320030 


0 .00033927 


0.90109152 



MEAfJ 

VARIANCE 



NUMBER OF BAUDS 
OF EPSILON-HAT 
OF EPSILON-HAT 



5 

0.863555133 
0 .8325092^;9E-02 
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ALPHA-HAT VS. ALPHA — INPUT SNR = 15 DB 

BAUD TYPE 5: KX = 102A SAMPLE POINTS; K = 6A TONES 



XX EPSILON = 0.90 DELTA-ALPHA = 0.000376506 



LL ALPHA 



ALPHA-HAT 



EPSILON-HAT 



1 0.000338855 

2 0.000338855 

3 0.000338855 

A 0.000338855 

5 0.000338855 



0.00038678 
0 . 00031A33 
0 .00035951 
0 . 00030^481 
0.00035939 



1 .02729702 
0.83A85985 
0 .95A85568 
0.80957907 
0 .95A53286 



NUMBER OF BAUDS = 5 

MEAN OF EPSILON-HAT = 0.91622A658 

VARIANCE OF EPSILON-HAT = 0 . 83226 1 1 A5E- 02 



XX EPSILON = 0.95 DELTA-ALPHA = 0.000376506 



LL ALPHA 



ALPHA-HAT EPSILON-HAT 



1 0.000357681 

2 0.000357681 

3 0.000357681 

A 0.000357681 

5 0.000357681 



0.00020686 
0 .00033393 
0 .00037992 
0 . 00032AA9 
0.000379A9 



0 .5A9A3258 
0.88705266 
1 . 009069AA 
0.86183619 
1 .00792885 



NUMBER OF BAUDS 
MEAN OF EPSILON-HAT 
VARIANCE OF EPSILON-HAT 



5 

0 .86306357A 
0 .353112072E-01 



XX EPSILON = 1.00 DELTA-ALPHA = 0.000376506 



LL ALPHA 



ALPHA-HAT EPSILON-HAT 



1 0.000376506 

2 0.000376506 

3 0.000376506 

A 0.000376506 

5 0.000376506 



0 .00003566 
0 .00016103 
0 .00020522 
0.0003A329 
0 .00039861 



0.09A7062A 
0.A2768973 
0 .5A507101 
0.91177A3A 
1 .05371677 



NUMBER OF BAUDS 
MEAri OF EPSILOri-HAT 
VARIANCE OF EPSILON-HAT 



5 

0.607591391 
0 . 1^88A1AA1 
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ALPHA-HAT VS. ALPHA — INPUT SNR = 15 DB 

BAUD TYPE 3: KX = 102A SAMPLE POINTS; K = 6A TONES 



EPSILON = 1.25 


DELTA-ALPHA = 0.000376506 


LL ALPHA 


ALPHA-HAT EPSILON-HAT 


1 0.000A70658 

2 0.000A70658 

3 0.000A70638 

A 0.000A70658 

5 0.000470638 


-0.00049591 -1.51713200 

0.00025919 0.68840146 

0 .00030580 0.81219985 

0.00005479 0.14555374 

0.00028259 0.75056285 


NUMBER OF 


BAUDS = 5 



MEAN OF EPSILON-HAT = 0.215912998 

VARIANCE OF EPSILON-HAT = 0.80A957628 



XX EPSILON = 1.50 


DELTA-ALPHA = 0.000376506 


LL ALPHA 


ALPHA-HAT EPSILON-HAT 


1 0.000564765 

2 0.000564765 

5 0.000564765 

4 0.000564765 

5 0.000564765 


-0.00082413 -2.18889141 
-0.00046605 -1.23785207 
-0.00020983 -0.55731559 
-0.00024753 -0.65690029 
-0 .00004868 -0 .12930280 



NUMBER OF BAUDS = 5 

MEAN OF EPSILON-HAT = -0.95A0A7918 
VARIANCE OF EPSILON-HAT = 0.652816792 
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C. INPUT SNR = 20 dB 



ALPHA-HAT VS. ALPHA — INPUT SNR = 20 DB 

BAUD TYPE 3: KX = 102A SAMPLE POINTS; K = 6A TONES 



EPSILON = 0, 


,0 DELTA-ALPHA 


= 0.000376506 


LL 


ALPHA 


ALPHA-HAT 


EPSILON-HAT 


1 


0.000 


0.00002055* 


0.05A57596 


2 


0.000 


-0.00001621 


-0.0A305860 


3 


0 . 000 


-0 .00000631 


-0 .01675819 


A 


0 .000 


-0.00001915 


-0.05085765 


5 


0.000 


0.00000311 


0.0082601A 




NUMBER OF BAUDS = 


5 




MEAN OF 


EPSILON-HAT = -0 . 956766680E-02 


VARIANCE OF 


EPSILON-HAT = 0 . 182760879E-02 



XX EPSILON = 0.25 
LL ALPHA 



DELTA-ALPHA = 0.000376506 
ALPHA-HAT EPSILON-HAT 



1 


0 . 00009A128 


0 .000116A3 


0.3092A392 


2 


0 . 00009A128 


0 . 00007815 


0 .20756A77 


3 


0 . 00009A128 


0 . 0000920A 


0 .2AAAA658 


A 


0 . 00009A128 


0 . 000073A8 


0 .19517177 


5 


0 . 00009A128 


0 .00009959 


0.26A51528 



MEAN 

VARIANCE 



NUMBER OF BAUDS 
OF EPSILON-HAT 
OF EPSILON-HAT 



5 

0 .2AA188309 
0.20973AA96E-02 



EPSILON = 0.50 


DELTA-ALPHA = 


0 .000376506 


LL 


ALPHA 


ALPHA-HAT 


EPSILON-HAT 


1 


0.000188255 


0 . 00021262 


0 .56A71968 


2 


0 .000188255 


0.00017358 


0 . A6102297 


3 


0.000188255 


0 .00019098 


0 .5072A1A3 


A 


0 .000188255 


0.00016762 


0 . AA520S97 


5 


0.000188255 


0 .00019673 


0 . 52252698 



NUMBER OF BAUDS 
MEAN OF EPSILON-HAT 
VARIANCE OF EPSILON-HAT 



5 

0 .5001A3766 
0 .2317A2300E-02 ' 
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ALPHA-HAT VS. ALPHA 
BAUD TYPE 3: KX 



INPUT SNR = 20 DB 
02A SAMPLE POINTS; K = 6A TONES 



XX EPSILON = 0.75 



DELTA-ALPHA = 0.000376506 



LL ALPHA 



ALPHA-HAT 



EPSILON-HAT 



1 0.000282383 

2 0.000282383 

3 0.000282383 

A 0.000282383 

5 0.000282383 



0.00030833 
0.00026935 
0.00028965 
0 . 00026252 
0.0002937A 



0.81892306 

0.71539086 

0.769311A9 

0.69726127 

0.78016865 



MEAN 

VARIANCE 



NUMBER OF BAUDS 
OF EPSILOf]-HAT 
OF EPSILOfJ-HAT 



5 

0.756210685 
0 .2A5A93A6AE-02 



XX EPSILON =1.00 



DELTA-ALPHA 



0.000376506 



LL ALPHA 



ALPHA-HAT 



EPSILON-HAT 



1 0.000376506 

2 0.000576506 

3 0.000376506 

A 0.000376506 

5 0.000376506 



0.000A0266 
0.00036A61 
0 . 00038703 
0 .00035726 
0 .00038970 



1 . 069A7517 
0 . 968A0972 
1 .0279A170 

0 .9A887573 

1 . 0350A181 



MEAN 

VARIANCE 



NUMBER OF BAUDS 
OF EPSILOf^HAT 
OF EPSILOf^-HAT 



5 

1 .0099A873 
0.2A8S05061E-02 



XX EPSILOfi 



05 DELTA-ALPHA 



000376506 



LL ALPHA 



ALPHA-HAT 



EPSILOf]-HAT 



1 0.000395336 

2 0.000395536 

3 0.000395336 

A 0.000395336 

5 0.000395336 



0 . 000A218A 
0 .00038A08 
0 . 000A068A 
0.00037666 
0 .000A0925 



1 . 120A166A 
1 . 02010632 
1 . 0805616A 
1 . 000A1580 
1 .08697891 



MEAN 

VARIAfJCE 



NUMBER OF BAUDS 
OF EPSILON-HAT 
OF EPSILON-HAT 



5 

1 . 06169510 
0 . 2A8205289E-02 ' 
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ALPHA-HAT VS. 
BAUD TYPE 3: 



ALPHA — INPUT SNR = 20 DB 
KX = 102A SAMPLE POINTS; 



K = 6A TONES 



XX EPSILON 1.10 



DELTA-ALPHA = 0.000376506 



LL 


ALPHA 


ALPHA-HAT 


EPSILON-HAT 


1 


0 . 000A1A161 


0 . 0O0AA105* 


1 . 171A1533 


2 


0 . 000A1A161 


0 . 00021188 


0.5627A819 


3 


0 . 000A1A161 


0 . 000A2667 


1 .13323879 


A 


0 . 000A1A161 


0.00039613 


1 . 05211163 


5 


0.000A1A161 


0 . 000A2883 


1 .13897A19 




NUMBER OF 


BAUDS = 


5 




MEAN OF EPSILON-HAT = 1.01169682 


VARIANCE OF EPSILON-HAT = 0.6A917 


50A8E-01 


XX EPSILON = 1.15 


DELTA-ALPHA = 


0 .000376506 


LL 


ALPHA 


ALPHA-HAT 


EPSILON-HAT 


1 


0 . 000A32987 


0.00026110 


0 .693A7852 


2 


0 . 000A32987 


0 . 0002313A 


0 .61AA5159 


3 


0 . 000A32987 


0 . 00025232 


0 .6701525A 


A 


0 . 000A32987 


0 . 00021900 


0.58167678 


5 


0.000A32987 


0 . OOOAA832 


1 .1907A15A 




NUMBER OF 


BAUDS = 


5 




MEAN OF EPSILON-HAT = 0.750099897 


VARIANCE OF EPSILON-HAT = 0 . 6 26 323223E-01 


XX EPSILOtJ = 1.20 


DELTA-ALPHA = 


0 .000376506 


LL 


ALPHA 


ALPHA-HAT 


EPSILON-HAT 


1 


0 . 000A51812 


0 . 00008A00 


0.22311223 


2 


0 . 000A51812 


0 . 0002A993 


0 .663S03A6 


3 


0 . 000A51812 


0 .00027113 


0 . 720110A8 


A 


0 . 000A51812 


0 . 0000A699 


0 . 12A80A02 


5 


0 . 000A51812 


0 . 000A6693 


1 .2A015808 




NUMBER OF 


BAUDS = 


5 



MEAN OF EPSILON-HAT = 
VARIANCE OF EPSILON-HAT = 



0 . 59A3975A5 
0 . 19899952A 
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ALPHA-HAT VS. ALPHA -- INPUT SNR = 20 DB 

BAUD TYPE 5; KX = 102A SAMPLE POINTS; K = TONES 



XX EPSILON = 1.25 



DELTA-ALPHA = 0.000576506 



LL ALPHA 



ALPHA-HAT EPSILON-HAT 



1 

2 

3 

5 



0. 000A70658 
0 .000^70638 
0 .000A70658 
0 . 000A70638 
0 . 000A70638 



0 . 00010312 
0 . 000269A7 
0 . 000101A8 
0 . 00006651 
0 . 000 ^ 86^6 



0.27389967 
0.71572357 
0.26952112 
0.17663908 
1 . 29209991 



MEAN 

VARIANCE 



NUMBER OF BAUDS 
OF EPSILON-HAT 
OF EPSILON-HAT 



5 

0.595565367 

0.218075991 



XX EPSILON = 1.50 



DELTA-ALPHA 



0.000376506 



LL ALPHA 



ALPHA-HAT 



EPSILON-HAT 



1 

2 

3 

9 

5 



0.000569765 
0 .000569765 
0 .000569765 
0 .000569765 
0 . 000569765 



-0.00062583 
-0 .00029093 
-0 .00022707 
-0.00023890 
-0.00006503 



-1.65688896 

-0.63990915 

-0.60310279 

-0.63318038 

-0.16791729 



NUMBER OF 
MEAN OF EPSIL 
VARIAf^CE OF EPSIL 



BAUDS 

ON-HAT 

ON-HAT 



5 

-0.790099311 

0.302176237 



142 



D. INPUT SNR = 40 dB 



ALPHA-HAT VS. ALPHA — INPUT SHR = AO DB 

BAUD TYPE 3: KX = 102A SAMPLE POINTS; K = 6A TONES 



EPSILON =0.0 


DELTA-ALPHA 


= 0.000376506 


LL 


ALPHA 


ALPHA-HAT 


EPSILON-HAT 


1 


0 . 000 


0 .00000226 * 


0 . 00601395 


2 


0 . 000 


-0.00000109 


-0 . 00289558 


3 


0 .000 


-0 . 00000033 


-0 .00100686 


4 


0 .000 


-0 . 000001A2 


-0 . 00376268 


5 


0 . 000 


0 . 00000092 


0.002A3571 




NUMBER OF 


BAUDS = 


5 




MEAN OF EPSILON-HAT = 0.156 


90870AE-03 


VARIANCE OF EPSILON-HAT = 0.165 


832AA7E-0^ 


EPSILON = 0.25 


DELTA-ALPHA 


= 0.000376506 


LL 


ALPHA 


ALPHA-HAT 


EPSILON-HAT 


1 


0 . 00009^128 


0 .00009623 


0.25558215 


2 


0 . 00009A128 


0 .00009276 


0.2^631^53 


3 


0 . 00009^123 


0 .00009578 


0.2^908829 




0 . 00009^128 


0 .00009219 


0 . 2^A8553A 


5 


0 .00009^128 


0 .00009687 


0 .2519687A 




NUMBER OF 


BAUDS = 


5 




MEAN OF EPSILON-HAT = 0.2^9561667 


VARIANCE OF EPSILON-HAT = 0 . 187 39A526E-0A 



EPSILON = 0.50 


DELTA-ALPHA = 


0 . 000376506 


LL 


ALPHA 


ALPHA-HAT 


EPSILON-HAT 


1 


0 .000188255 


0 . 00019089 


0.50700271 


2 


0 . 000188255 


0 .00018733 


0 .49755794 


3 


0 .000188255 


0 .00018857 


0 .50083911 


4 


0 .000188255 


0 .00018658 


0.49555242 


5 


0 .000188255 


0 . 00018962 


0.50309473 




NUMBER OF 


BAUDS = 


5 




MEAN OF EPSILON-HAT = 0.5008 


09073 



VARIANCE OF EPSILON-HAT = 0 . 20 A A7 1 289 E-0 A ’ 
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ALPHA-HAT VS. ALPHA 
BAUD TYPE 3: KX 



- INPUT SUR = AO DB 
102A SAMPLE POINTS; K 



6A TONES 



XX EPSILON = 0.75 



DELTA-ALPHA = 0.000376506 



LL ALPHA 



ALPHA-HAT EPSILON-HAT 



1 0.000282383 

2 0.000282383 

3 0.000282383 

A 0.000282383 

5 0.000282385 



0 . 000285A7 
0 . 0002819A 
0 .00028319 
0 .00028101 
0 .00028383 



0.75821120 
0.7A88AA33 
0 .752161AA 
0.7A635839 
0.75386528 



MEAN 

VARIANCE 



NUMBER OF BAUDS 
OF EPSILON-HAT 
OF EPSILON-HAT 



5 

0.751887321 
0 .209A99267E-OA 



EPSILON = 1.00 


DELTA-ALPHA = 


0 . 000376506 


LL 


ALPHA 


ALPHA-HAT 


EPSILON-HAT 


1 


0 .000376506 


0 . 00037910 


1 .0068893A 


2 


0 . 000376506 


0 . 00037569 


0.99783307 


5 


0 .000376506 


0 . 00037675 


1 . 00063A19 


A 


0.000376506 


0 . 00037A61 


0 .99A96A96 


5 


0 .000376506 


0 .00037728 


1 .00205612 



NUMBER OF BAUDS 
MEAN OF EPSILOfJ-HAT 
VARIANCE OF EPSILON-HAT 



5 

1 . 000A7A93 
0 .202523515E-0A 



XX EPSILON =1.05 


DELTA-ALPHA 


= 0.000376506 


LL 


ALPHA 


ALPHA-HAT 


EPSILON-HAT 


1 


0 .000395336 


0 .00039818 


1 . 05757332 


2 


0.000395336 


0 . 00039A80 


1 .0A860115 


3 


0 .000395336 


0 .00039580 


1 . 0512A578 


A 


0 .000395336 


0 . 00039370 


1 .0A568005 


5 


0 .000395336 


0 . 0005963A 


1 .05267A29 



MEAti 

VARIANCE 



NUMBER OF BAUDS 
OF EPSILON-HAT 
OF EPSILON-HAT 



5 

1.05115A1A 
0 .200020AA0E-0A 
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ALPHA-HAT VS. ALPHA — INPUT SNR = ^0 DB 

BAUD TYPE 3: KX = 102A SAMPLE POINTS; K = 6^^ TONES 



XX EPSILON = 1.10 DELTA-ALPHA = 0.000376506 



LL 


ALPHA 


ALPHA-HAT 


EPSILON-HAT 


1 


0 . 000<4lA161 


0 .000417 30 


1 .10836124 


2 


0 . 000A1A161 


0 .00041396 


1 .09947 109 


3 


0 . 000A1A161 


0 .00041488 


1 .10192490 


4 


0 . 000A1A161 


0 .00041284 


1 . 09649086 


5 


0 . 00041A161 


0.00041541 


1.10333538 



MEAN 

VARIANCE 



NUMBER OF BAUDS 
OF EPSILOtJ-HAT 
OF EPSILON-HAT 



5 

1.10191631 

0.197^4l3901E“0A 



XX EPSILON =1.15 



DELTA-ALPHA = 0.000376506 



LL ALPHA 



ALPHA-HAT 



EPSILOfHHAT 



1 0.000A32987 

2 0.000432987 

3 0.0004329S7 

4 0.000432987 

5 0.000432987 



0 .00043634 
0 .00043304 
0 .00043389 
0 .00023536 
0 .00043442 



1 . 15891647 
1.15016556 
1 .15240288 
0.62512481 
1 . 15381145 



MEAN 

VARIANCE 



NUMBER OF BAUDS 
OF EPSILON-HAT 
OF EPSILON-HAT 



5 

1 . 04808331 
0 .5591490S6E-01 



XX EPSILON = 1.20 



DELTA-ALPHA 



0 .000376506 



LL ALPHA 



ALPHA-HAT 



EPSILOfi-HAT 



1 0.000451812 

2 0.000451812 

3 0.000451812 

4 0.000451812 

5 0.000451812 



0 . 00045451 
0 .00045124 
0 .00026262 
0 .00006299 
0 .00045256 



1 .20717335 
1 .19849777 
0 .69750923 
0 .16730273 
1 .20199013 



MEAN 

VARIANCE 



NUMBER OF BAUDS 
OF EPSILON-HAT 
OF EPSILON-HAT 



5 

0.894494414 

0.213087618 
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ALPHA-HAT VS. ALPHA -- INPUT SNR = ^^0 DB 

BAUD TYPE 3: KX = 102A SAMPLE POINTS; K = 6A TONES 



XX EPSILON = 1.25 DELTA-ALPHA = 0.000376506 



LL 


ALPHA 


ALPHA-HAT 


EPSILON-HAT 


1 


0 . 000A7063S 


0 . 000A7359 


1.25785255 


2 


0 . 000A70638 


0 . 00008576 


0 .227779A5 


3 


0 . 000A70638 


0 .00028165 


0 .7A806517 


A 


0 . 000A70638 


-0 . 0001169A 


-0.31059778 


5 


0 . 000A70638 


0 . 000A7162 


1.25261021 




NUMBER OF BAUDS = 


5 




MEAN OF EPSILON 


-HAT = 0.6351A1730 


VARIANCE OF EPSILON 


-HAT = 0.A60538337 


EPSILON =1.50 


DELTA-ALPHA 


= 0.000376506 


LL 


ALPHA 


ALPHA-HAT 


EPSILON-HAT 


1 


0 . 00056A765 


-0 . 000AA0A2 


-1 . 16975A03 


2 


0 . 00056A765 


-0 . 00001852 


-0 . 0A919A85 


5 


0 . 00056A765 


-0 . 0002AA38 


-0 .6A908123 


A 


0 . 00056A765 


-0 .00022513 


-0 .59793A7S 


5 


0 . 00056A765 


-0 . 000063A5 


-0 .16852963 



NUMBER OF 
MEAN OF EPSIL 
VARIANCE OF EPSIL 



BAUDS 

ON-HAT 

Ofi-HAT 



5 

-0 .5268987A2 
0 . 197A66671 
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APPENDIX F. THE SIMULATION CODE: RXSIM 



C THIS SIMULATIOri PRODUCES A PREDICTIOfl OF THE RECEIVED SIGfJAL THROUGH 
C AN ACOUSTIC CHANNEL. THE TRAfiSMITTED SIGNAL IS MFQPSK SIGNAL. 

C 

REAL NU(100,1350), PHI ( 1 00 . 1 350 ) , Ul, DELUl 
REAL FREQ, FF, ABSORP, AAUNDB, TL, AA(100,I550) 

REAL YRX, YY, T I ME( A 1 0 0 0 0 ) , YYRX(AIOOOO) 

REAL RXKM(0:A096 ), RXKP ( 0 : A 0 96 ) , RXKPDC 0 : A096 ) , KFREQ ( 0 : A 0 96 ) 

REAL DELT(IOO), DELF(IOO) 

REAL VX(IOO), X(IOO) 

REAL VY(IOO), Y(IOO) 

REAL VZ(IOO), Z(IOO) 

REAL C(IOO), R(IOO), THETAD(IOO) 

REAL TAUL(IOO), ALPHA(IOO), ALPHAM(IOO), DELALF(IOO) 

REAL UHATKIOO), LBAUD(IOO), PTOT(IQO), NOSVAR(IOO) 

REAL TAUKCI350), NOSAVG, NOISE 

INTEGER BDTYPE(IOO), BDTOTL, NBAUDSCSO), IDFT 

INTEGER KMIN(IOO), KMAX(IOO), KX(IOO), M(IOO), KPTS(IOO) 

INTEGER IIPHKIOO, 1350) , IPHI 

COMPLEX RXIN(0:A096), RXOUT ( 0 : A096 ) 

DOUBLE PRECISION DSEED 

CHARACTER^AO PLABEL 

C 

C USER INPUTS FOR THE SIMULATION 
C 

HRITE(6, 1000) 

READ(5,^) XO, YO, ZO 
C 

HRITE(30,1000) 



NRITE(30,^) * XO = * ,X0, ' YO = 
MRITE(30,5() » * 


' , YO, * ZO 


= SZO 


HRITE(6,100I) 

READ(5,X) VXAVG, VXVAR, VYAVG, 


VYVAR, VZAVG, VZVAR 


t-lRITE( 30, 1 001 ) 

WRITE(30,^) * VXAVG = VXAVG,’ 
MRITE(30,X) ‘ VYAVG = VYAVG,’ 
HRITE(30,X) » VZAVG = *,VZAVG,’ 
11RITE(30,X) * * 


' VXVAR = 
' VYVAR = 
' VZVAR = 


‘ , VXVAR 
' , VYVAR 
* , VZVAR 


HRITEC6, 1002) 
READ(5,X) CO, CVAR 






l!RITE( 30, 1 002) 

NRITEC30, ?^) ' CO = * ,C0, * CVAR 


= *,CVAR 





11RITE(30,K) » * 

C 

MRITEC6, 1003) 

READ(5,X) THETAO 
C 

1‘1RITE( 30, 1 003 ) 

tlRITE(30,x) » THETAO = THETAO 
HRITE(30,^) » * 

C 

HRITEC6, lOOA) 

READ(5,X) lAMP 



NRITE(30, lOOA) 

l'IRITE(30, X) * lAMP = MAMP 

HRITEC30, X) * • 

INITIALIZE VARIABLES 

BDTOTL = THE TOTAL NUMBER OF BAUDS IN THE SIGNAL 
RO = CLOSEST POiriT OF APPROACH 
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c 



c 



c 



c 



c 



c 



c 



c 



c 



c 



c 



DDTOTL = 0 
RO = 10MK15.0 
X(l) = XO 
Y(l) = YO 
Z(l) = ZO 



IF 


(XO 


.GT . 


0) 


VXAVG 


IF 


(XO 


.LE. 


0 ) 


VXAVG 


IF 


(YO 


.GT . 


0 ) 


VYAVG 


IF 


(YO 


.LE. 


0) 


VYAVG 


IF 


(ZO 


.GT . 


0) 


VZAVG 


IF 


(ZO 


. LE. 


0) 


VZAVG 



-ABS( VXAVG) 
ABS( VXAVG) 
-ABS( VYAVG) 
ABS( VYAVG) 
-ABS(VZAVG) 
AB5( VZAVG) 



HRITEC6, 1010) 
READ(5,X) NPAKS 



l‘IRITE(30a010) 

WRITE(30,X) * MPAKS = NPAKS 
NRITE(30,X) * » 

10 HRITEC6, 1100) 

WRITE(6ai01 ) 
llRITE(6ai02) 

WRITE(6, 1103) NPAKS 
READ(5,^) METHOD 



WRITE(30,1100) 

NRiTE(3o, non 

WRITEC 30, 1102) 

NRITEC 30,1103) NPAKS 
HRITE(30,^) * METHOD = », METHOD 
NRITE(30,X) * ' 



IF (METHOD .LE. 0) GO TO 10 
IF (METHOD .GE. 3) GO TO 10 
DO 100 I = 1, NPAKS 
50 MRITE(6, 1020) 

l!RITE(6, 1021 ) I 
IIRITE(6, 1 022) 

UR1TE(6, 1023) 

MRITE(6,102^) 

WRITE(6, 1025) 

MR1TE(6, 1026) 

READ(5,^) BDTYPE(I) 

HRITE(30, 1020) 

URITE( 30, 1 021 ) I 
NRITE(30, 1022) 

IIRITE( 50,1 025 ) 

URITE( 30, 1 02<4) 

NRITE( 50, 1025) 

URITE(30, 1026) 

URITE(50,>) ' BDTYPE(\I,n = ’,BDTYFE(I) 

HRITE(30, ^) ‘ * 

IF ( BDTYPE( 1 ) .LE. 0) GO TO 50 

IF ( BDTYPE( 1 ) .GT . 5) GO TO 50 

60 HRITE(6, 1027) I 

READ(5,^) NBAUDS(I) 

11RITE( 30,1 027 ) I 

NRITE(30,X) * NBAUDS(M,') = *,NBAUDS(I) 
HRITE(30,X) » ' 

IF (NBAUDS(I) .LE. 0) GO TO 60 
100 COfJTINUE 

HRITE(6, 1030) 

READ(5,^) IDFT 

MRITE( 30, 1030) 

NRITE(30,^) ’ IDFT = MDFT 
HR1TE(30,^.) ' * 
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C 



C 



C 



C 



C 

C 



IF (IDFT .EQ. I) THEM 
MRITEC6, 10^0) 

READ(5,X) IHNDON 

NRITE(30,10A0) 

HRITE(30,X) ' IHNDON = ',IHflD0H 
HRITE(30,X) • » 

ENDIF 

HRITE(6, 1050) 

READ(5,^) SNRDB 

HRITE(30, 1050) 

NRITEC 30, * ( AlO, FIO . A) * ) * StmOB = 

NRITE(30,X) * • 

SNRIN = 10.0X^(SNRDB/I0.0) 

LL = 0 

DO 300 I = I, NPAKS 

DO 200 J = I, NBAUDS(I) 

LL = LL + I 

IF (BDTYPE(I) .EQ. I) THEM 
DELT(LL) = 1,0 / 2A0.0 
DELF(LL) = 2A0.0 
KniN(LL) =68 
KMAX(LL) =83 
KX(LL) = 256 
ENDIF 

IF (BDTYPE(I) .EQ. 2) THEN 

DELT(LL) = 1.0 / 120.0 
DELF(LL) = 120.0 
KMItKLL) = 135 
KMAX(LL) = 166 
KX(LL) = 512 
EflDIF 

IF (BDTYPE(I) .EQ. 3) THEN 

DELT(LL) =1,0/60.0 
DELF(LL) = 60.0 
KMIN(LL) = 269 
KMAX(LL) = 332 
KX(LL) = I02A 
ENDIF 

IF (BDTYPE(I) .EQ. A) THEN 

DELTC LL ) = 1.0/30.0 
DELF(LL) = 30.0 
KMIN(LL) = 537 
KMAX(LL) = 66A 
KX(LL) = 20A8 
ENDIF 

IF (BDTYPE(I) .EQ. 5) IHEU 

DELT(LL ) =1.0/15.0 
DELF(LL) = 15.0 
KMIN(LL) = 1073 
KMAX(LL) = 1328 
KX(LL) = A096 
EliDIF 

200 COtITIflUE 
300 CONTINUE 



, SNRDB 



INITIALIZE VARIABLES 

PI = A.OXATAN(I.O) 

BDTOTL = LL 

IF (BDTOTL .GT. 100) THEN 

HRITE(6,X)» XXX ERROR: TOTAL fiUMDER OF BAUDS EXCEEDS 100 * 
NRITE(30,X)» XXX ERROR: TOTAL NUMBER OF BAUDS EXCEEDS 100 * 
GO TO 9999 
ENDIF 

ALLOW THE USER TO ENCODE THE PHASES 
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c 



IF (METHOD .EQ. 2) THEN 
LL = 0 

DO ^490 1 = 1, MPAKS 

NRITEC6, 1110) I 
NRITE(6,1111) I, I 
AlO READ(5,X) MTHDPK 

C 

WRITE(30, 1110) I 
HRITE(30,1111) I, I 
HRITE(30,X) • MTHDPK = *, MTHDPK 
NRITE(30,^) • • 

C 

IF (MTHDPK .LE. 0) GO TO AlO 

IF (MTHDPK .GE. 3) GO TO AlO 

IF (MTHDPK .EQ. 1) THEM 

DO A30 J = 1, MBAUDS(I) 

LL = LL + 1 

DO A20 K = KMIM(LL), KMAX^LL) 
DSEED = (K ^ 3.5729) + DSEED 
CALL PHASE( DSEED, RtIDPHi, IPHI ) 
PHI(LL,K) = RMDPHI 
IIPHI(LL,K) = IPHI 
A20 CONTINUE 

A30 CONTINUE 

ENDIF 

IF (MTHDPK .EQ. 2) THEN 

DO A80 J = 1, NBAUDS(I) 

LL = LL + 1 

NRITE(6, 1120) I, J 
MRITE(6,1121 ) I, J, I, J 
AAO READ(5,X) MTHDDD 

C 

URITE(30, 1120) I, J 
WRITE(30, 1121 ) I, J, I, J 
WRITE(30,X) * MTHDBD = *,MTHDBD 
WRITE(30,X) * ‘ 

C 

IF (MTHDBD .LE. 0) GO TO AAO 
IF (flTHDBD .GE. 3) GO TO AAO 
IF (MTHDDD . EQ . 1) THEN 

DO A50 K = KMIN(LL), KMAX(LL) 
DSEED = (K ^ 5.7317) + DSEED 
CALL PHASE(DSEED,RNDPHI,IPHI) 
PH1(LL,K) = RNDPHl 
I1PHI(LL,K) = IPHI 
A50 COriTlNUE 

ENDIF 

IF (MTHDBD . EQ . 2) THEN 
URITE(6, 1100) 

MRITE(6, 1101 ) 

C 

HRITE( 30 ,n 00 ) 

MRITE(30,1101) 



C 



A60 



C 



C 



A70 



DO A70 K = KMIN(LL), KMAX(LL) 
MRITE(6, 1130) I, J, K 
READ(5,>^) IPHI 
IIPHI(LL,K) = IPHI 



IIRITE( 30, 1 1 30) I, J, K 
MR1TE(30.>^) ’ IPHI = *,IPH1 
MR1TE(30,?0 ' 



IF (IPHI . LE. 0) GO TO A60 
IF (IPHI .GE. 5) GO TO A60 
IF (IPHI .EQ. 1) P[1I(LL,K) = 
IF (IPHI .EQ. 2) PH1(LL,K) = 
IF (IPHI .EQ. 3) PH1(LL,K) = 
IF (IPHI .EQ. A) PHI(LL,K) = 
CONTif^UE 
ENDIF 



(A5.0 ^ PI) / 
(135.0 K PI) 
(-135.0 ^ PI) 
(-A5.0 



180.0 
180.0 
/ 180.0 
180.0 
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^80 CONTINUE 

ErJDIF 
COfUINUE 
EI^DIF 

DO 550 LL = 1, BDTOTL 

INITIALIZE VARIABLES 

DSEED = (LL ^ 15.5987) + D5EED 
LBAUD(LL) = LL 

COMPUTE* 

X(LL) = DISTAr^CE THE TRAf^SMITTER TRAVELED l\l THE 
X-DIRECTIOr^ DURING THE LLTH BAUD 

CALL GAUSS ( DSEED, VXAVG,VXVAR,ZRt^D) 

VX(LL) = ZRND 

IF (LL .GT. 1) X(LL) = X(LL-l) + (VX(LL) DELT(LD) 



COMPUTE: 

Y(LL) = DISTAt^CE THE TRAt^SMITTER TRAVELED ir^ THE 
Y-DIRECTIOt^ DURIt^G THE LLTH BAUD 

CALL GAUSSC DSEED, VYAVG,VYVAR,ZRrm) 

VY(LL) = ZRND 

IF (LL .GT. 1) Y(LL) = Y(LL-l) + (VY(LL) DELT(LD) 



COMPUTE: 

Z(LL) = DISTAf^CE THE TRAt^SMITTER TRAVELED IN THE 
Z“DIRECTIOI^ DURING THE LLTH BAUD 

CALL GAUSS( DSEED, VZAVG,VZVAR,ZRf^D) 

VZ(LL) = ZRf^D 

IF (LL .GT. 1) Z(LL) = Z(LL-I) + (VZ(LL) DELT(LD) 



COMPUTE: 

C(LL) = THE SPEED OF SOUt^D DURItiG THE LTH BAUD 

CALL GAUSS(DSEED,CO,CVAR,ZRND) 

C(LL) = ZRtiD 

COMPUTE: 

R(LL) = SLAf^T RAt^GE TO THE RECEIVER 

XXX = ADS(X(l.L)) 

YYY = ABS(Y(LL ) ) 

ZZZ = ABS(Z(LD) 

R(LL) = ( XXX>'/2.0 + YYY^x2.0 + ZZZ?^>2.0 )^xo.5 
IF (R(LL) .LE. RO) RO = R(LL) 

COMPUTE: 

THETAD(LL) = ANGLE BETWEEN THE SLANT RANGE AND ZO 

ARG = ABS(Z(LD) / R(LL) 

THETAD(LL) = ACOS(ARG) (180.0/PI) 

IF (THETAD(LL) .GT. THETAO) THEN 
HRITE(6,X) 

X ’ XX WARNIt^G: THETA IS GREATER THAN THETAO,' DEGREES 

WRITE(6,K) 

X ' xxxxxxxxxxx the receiver CAt^ tlOT RECEIVE THE SIGNAL 
HRITE(6,x)» 

WRITE(30,x) 

X ' XX WARNING: THETA IS GREATER THAt^ ', THETAO,' DEGREES 

NRITE(30,X) 

% • xxxxxxvxxxx THE RECEIVER CAt^ NOT RECEIVE THE SIGf^AL 

WRITE(30,X)’ 

ENDIF 



XX ♦ 
X X » 

XX • 
X X ’ 
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c 

c 



INITIALIZE LIMITS FOR GRAPHS 



C 

C 

C 



C 

C 

C 



C 



C 

C 

C 

C 

C 

C 



C 

C 

C 

C 

C 

C 

C 

C 

c 



c 

c 

c 



IF (LL .EQ. 1) THEM 
XXMIN = INTCXCD) 

XXMAX = INTCXC 1 ) ) + 1.0 
YYMIN = IMTCYC 1 ) ) 

YYMAX = IMTCYC 1) ) + I .0 
ZZMIM = IMTCZC 1 ) ) 

ZZMAX = IMTCZC I ) ) + I . 0 
CCMIM = IMTCCC 1) ) 

CCMAX = IMTCCCD) + 1.0 
RRITIM = IMTCRCI ) ) 

RRMAX = IMTCRCI )) + 1.0 
THTMIM = IMTCTHETADC 1 ) ) 

THTMAX = IMTCTHETADC!)) + 0.5 
EMDIF 



SET LIMITS FOR GRAPHS 



IF 


CXCLL) 


.LT . 


XXMIN) 


XXMIN 


= XCLL) 




IF 


CXCLL) 


.GT. 


XXMAX) 


XXMAX 


= XCLL ) 




IF 


CYCLE) 


.LT. 


YYMIfO 


YYMIN 


= YC LL ) 




IF 


CYCLE) 


.GT. 


YYMAX) 


YYMAX 


= YC LL ) 




IF 


CZCLL) 


.LT. 


ZZMIM) 


ZZMIM 


= ZCLL) 




IF 


CZCLL) 


.GT. 


ZZMAX) 


ZZMAX 


= ZCLL) 




IF 


CCC LL ) 


.LT. 


CCMIfO 


CCMIfi 


= CC LL ) 




IF 


CCCLL ) 


.GT . 


CCMAX) 


ccr.AX 


= CCLL) 




IF 


CRC LL ) 


. LT . 


RRMIN) 


RRMIN 


= RC LL ) 




IF 


CRCLL) 


.GT . 


RRMAX) 


RRMAX 


= RC LL ) 




IF 


CTHETADC LL ) 


.LT. THTfllN) 


THTMIM = 


THETADC LL ) 


IF 


CTHETADCLL) 


.GT. THTMAX) 


THTMAX = 


THEIADCLL) 


ASSIGN 


ALL THE 


PHASES RAIIDOtlLY 


FOR EVERY 


PACKET 


DO 


50 0 K 


= KMINCLL), 


KMAXCLD 





IF CMETHOD . EQ . 1) THEM 

DSEED = CK X 2.767653) + D5EED 
CALL PHASECDSEED,RfiDPHI,IPHI) 

PHICLL.K) = RMDPHI 
IIPHICLL.K) = IPHI 
EMDIF 

500 COfUIMUE 
550 CONTINUE 

DO 660 LL = 1, DDTOTL 

COMPUTE: 

ALPHACLL) = THE DOPPLER COMPRESSlOfJ FACTOR 
DUE TO THE MOVING TRAfJSMITTER 

AVX = VXAVG X XC L L ) 

AVY = VYAVG ^ YC LL ) 

AVZ = VZAVG X ZCLL) 

ALPHACLL) = CAVX + AVY + AVZ) / C R C L L ) >^C C L L ) ) 

COMPUTE: 

DELALFCLL) = THE MAXIMUM CHANGE IN ALPHA FOR THE LTH BAUD 
MCLL) = THE DOPPLER CHANNEL fi NEEDED FOR THE LTH BAUD 
ALPHAMCLL) = THE DOPPLER CHANNEL FACTOR FOR THE LTH BAUD 
NOTE: 

ALPHAMCLL) IS USED TO ESTIMATE THE SAMPLING 
FREQUENCY OF THE RECEIVED SIGfmL 

DELALFCLL) = 1 / C8.0 X KMAXCLD) 

AMIN = CALPHACLL) / DELALFCLL)) - 0.5 
AMAX = CALPHACLL) / DELALFCLL)) + 0.5 
MCLL) = iraCCAMIN + AMAX) / 2.0) 

ALPHAMCLL) = MCLL) X DELALFCLL) 

COMPUTE: 

TAULCLL) = THE TIME DELAY DUE TO THE SIGNAL TRAVELING 
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c 

c 

c 



c 

c 

c 

c 

c 

c 

Cl 

c 

c 

c 

c 



c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



c 

c 

c 



THROUGH THE HEDIUM OR CHANNEL AflD THE RECEIVER 
TOR THE LTH BAUD 

TAUL(LL) = R(LL) / C(LL) 

IF (LL .EQ. 1) JHEU 
ALFfllt^ = INTCALPHAC 1 ) ) 

ALFMAX = INT(ALPHAd)) + 0.1 
ENDIF 

IF (ALPHA(LL) .LT. ALFMIN) ALFMIN = ALPHA(LL) 

IF (ALPHA(LL) .GT. ALFMAX) ALFMAX = ALPHA(LL) 

PTOT(LL) = 0.0 

NOSVARCLL) = 0.0 

DO 6A0 K = KMIN(LL), KMAX(LL) 



COMPUTE: 

Ul = THE RECEIVE SIGNAL SAMPLE START TIME 
DELUl = THE SYNCHRONIZATION ERROR 

NRITE(6, llAO) LL, K 

lAO F0RMAT(/,2X, ‘ENTER THE ESTIMATED RECEIVE SlGflAL START TIME 



' ,/,2X, 



» FOR BAUD NUMBER dIA,/,AX, 

» AND THE TIME DELAY THROUGH THE TRANSMITTER *,/,AX, 

* ELECTRONICS FOR TONE [DUMBER ... ») 

READ(5,)^) UHATKLL), TAUK(K) 

TAUK(K) = 0.005 

Ul = TAUL(LL) + (TAUK(K) 5^ ( 1 + ALPHA(LL))) 

UHATKLL) = Ul 

DELUl = Ul - UHATKLL) 

COMPUTE: 

NU(LL,K) = A RAf^DOM SEQUENCE OF TIMING JITTERS 
NU(LL,K) = (DELUl/DELTCLD) X ( KX ( L L )/ ( 1 + AL PHAMC L L ) ) ) 
COMPUTE: 

AA(LL,K) = THE AMPLITUDE OF THE RECEIVED SIGNAL, NHICH IS 
I FOR ALL LL AND K IF NORMALIZED AMPLITUDES ARE 
DESIRED (I.E., IAMP = 1), OR ATTEfiUATlON FACTOR 
DUE TO THE TRANSMISSION LOSS IF NORMALIZED 
AMPLITUDES ARE NOT DESIRED (I.E., 1AMP=0) 

IF (lAMP .EQ. 0) THEN 

COMPUTE : 

ABSORP = THE AMOUNT OF ABSORPTIOfJ l\l DB/FT AT A DEGREES C 
AT A DEPTH OF APPROXIMATELY 3000 FT 

FREQ = (K DELF(LD) / 1 000.0 
FF = FREQ5^^2.0 

ABSORP = ( 0.0 03 + ( (O.KFF) / (KFF) ) 

+ ( ( AO . 0?^FF) / ( A . 1 00 + FF) ) 

X +( 0 . 000275 ^ FF ))/ 3000 

TXDEP = THE DEPTH THE TRANSMITTER IS BELOH THE SURFACE If^ FEET 

TXDEP = 1000.0 

ZTOTAL = TXDEP + ABS(Z(LD) 



AAUNDB = 10 (ABSORP / 20.0) 

IFAC = INT((ZTOTAL “ 3000.0) / 1000.0) 
IF (IFAC .LT. 0) FACTOR =1.02 
IF (IFAC .GT. 0) FACTOR =0.98 
IIFAC = ABS(IFAC) 

IF (IIFAC .NE. 0) THEN 
DO 620 1=1, IIFAC 

AAUNDB = AAUNDB 5^ FACTOR 
620 CONTINUE 

ENDIF 

ABSORP = 20.0 ^ ALOG10( AAUNDB) 

C 

C COMPUTE; 
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TL = THE TRANSMISSION LOSS IN DB DUE TO SPHERI CAL -SPREADI NG 
At^D ADSORPTIOfJ OF THE ACOUSTIC CHANNEL 

TL = (20.0 ALOGIO(ZTOTAL) ) + (AB30RP ZTOTAL) 

AA(LL,K) = 10 . 0KX(-(TL/20 . 0) ) 

ELSE 

AA(LL,K) =1.0 
EfJDIF 

COMPUTE THE fJOISE VARIANCE FOR THE DESIRED 
HIDE DAt^D SIGNAL-TO-tiOISE RATIO 

PTOT(LL) = THE TOTAL POHER OF THE LLTH BAUD 
NOSVAR(LL) = THE NOISE VARIANCE OF THE LLTH BAUD 

PTOT(LL) = PTOT(LL) + ( ( AA ( L L , K ) XX2 )/2 . 0 ) 

6A0 CONTINUE 

KPTS(LL) = (KMAX(LL) - KMIN(LD) .+ 1 

NOSVAR(LL) = (PTOT(LL) ^ KX(LD) / (2.0 StJRIN X KPTS(LD) 

660 CONTINUE 

II = 0 
RXMAX = 0.0 
RXMIN =0.0 
NOSAVG = 0.0 
NOISE = 0.0 

DO 700 LL = 1, BDTOTL 
DO 690 = I, KX(LL) 

DSEED = (LL X 32.7675) + DSEED 

CALL GAUSS( DSEED, NOSAVG, NOS VAR(LL) , NOISE) 

YRX =0.0 

\\U = N - I 

II = II + I 

TIME(II) = II - I 

DO 660 K = KMIN(LL), KMAX(LL) 



COMPUTE: 

YRX = THE RECEIVED SIGNAL HITH ALL OF THE ABOVE 
PARAMETERS COMBINED FOR THE LTH BAUD 
AT TIME = UU OVER ALL KMIN TO KMAX FREQUENCIES 

YY = AA(LL,K) x COS( ( ( ( 2XPI XK )/KX ( L L ) ) X 
((.1 + ALPHAM(LL)) / ( I + AL PHA ( L L ) ) ) X 
(NN - NU(LL,K))) + PHI(LL,K) ) 

YRX = YRX + YY 
COfJTINUE 

YYRXdl ) = YRX + NOISE 

IF (YYRXdl) .GT. RXMAX) RXMAX = YYRXdl) 

IF (YYRXdl) .LT. RXMIN) RXMIN = YYRXdl) 

CONTINUE 
7 00 COfniNUE 
C 

fiPTS = II 
KLL = 0 

DO 710 LL = I, BDTOTL 
KLL = KLL + KPTS(LL) 

710 COflTINUE 

HRITE(30,2000) KLL 

HRITE(30, 2001 ) ’ LL K PHI(LL,K) IPHI ' 

DO 7A0 LL = 1, BDTOTL 

DO 720 K = KMIN(LL), KMAX(LL) 

HRITE(30,2010) LL, K, PHI(LL,K), IIPHI(LL,K) 

720 CONTINUE 
7A0 COrUINUE 
C 

C COMPUTE THE DFT OF THE OUTPUT (RECEIVED) SIGNAL 
C 

IF (IDFT .EQ. 1) THEN 
RXKMIN =0.0 
RXKMAX =0.0 



680 



690 
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DO 800 LL = I, DDTOTL 
II = 0 

IF CLL .EQ. I) THEN 
JJ = I 
ELSE 

JJ = JJJ + 1 
ENDIF 

KXMI = KX(LL) - I 
JJJ = JJ 4 Kxm 
DO 760 J = JJ, JJJ 

RXimil) = CMPLXCYYRXC J), 0.0) 

II = II 4 I 
760 COMTIfJUE 

CALL DFT(KXMI,RXin,RXOUT) 

C 

C CONVERT OUTPUT DATA TO EXPONENTIAL FORM, I.E., RXKMO XEXP( JmRXKPC ) ) 

C 

DO 770 I = 0, KXMI 

RXKM(I) = SQRTC ( REAL ( RXOUT ( I J ) ) 2 4 ( AIMAG( RXOUT ( I) ) ) ) 

IF (ABS(REAL(RXOUT( I) ) ) .GE. I.E-I5) THEN 

RXKP(I) = A1 AN2(AIMAG(RX0UT(I) ), REAL ( RXOUT ( I) ) ) 

ELSE 

IF (ABSCAIMAGCRXOUTCI))) .LE. I.E-I5) RXKPCI) = 0.0 
IF (AIMAG(RXOUTd)) .GT. I.E-I5) RXKP(I) = PI/2.0 
IF ( AIMAG(RXOUT( I) ) .LT. -1.E-I5) RXKPCI) = -PI/2.0 
ENDIF 

RXKPD(i) = (RXKPCI) X 180.0) / PI 
770 CONTINUE 

C 

C HRITE DFT INPUT AfJD OUTPUT TO A FILE 
C 

l’IRITE( 30,2200 ) LL 
DO 780 I = 0, KXMI 

HRITE(30,22I0) I, REA L C RXI N C I ) ) , AIMAGC RXI N ( I ) ) 

780 COflTirJUE 

IIRITE(30,2220 ) LL 
IF (INNDOl’l .EQ. 0) THEN 
IMIN = 0 

I MAX = KXMI 

MM IN = 0 

MMAX = KXMI 

ELSE 

IMItl = KMIflCLL) 

I MAX = KMAXCLL) 

MMIN = KMINCLL) 

UMAX = KMAXCLL) 

EfJDIF 

KXX = (IMAX - IMIN) 4 1 
DO 7 90 I = IMIfJ, IMAX 

NRI TEC 30, 2230) I, REAL (RXOUT C I ) ) , AIMAGC RXOUT ( I ) ) , 

% RXKMCI), RXKPDCI), RXKPCI) 

IF (RXKMCI) .LT. RXKfUN) RXKMIN = RXKMCI) 

IF (RXKMCI) .GT. RXKMAX) RXKMAX = RXKMCI) 

KFREQCI) = I 
790 COrUINUE 

C 

C PLOT THE DFT OUTPUT 
C 

MRITE(PLABEL,2300) LL 
IF (LL .EQ. 1) CALL COMPRS 
CALL PAGE(11,8.5) 

CALL flODRDR 
CALL AREA2D(8,6) 

CALL XNAMEC ’ FREQ (K) $ MOO) 

CALL YNAMEC* MAGtHTUDE $M00) 

CALL HEADIfIC' DFT OUTPUT OF THE RECEIVED SIGfML $',100, ^,2) 
CALL HEADINC PLABEL , 100, 3, 2) 

CALL GRAFCMMIN, 'SCALE* , MMAX , RXKMI fl , 'SCALE* , RXKMAX) 

CALL GRIDCI,!) 

CALL SETCLRC 'MAGErUA' ) 

IF CKXX LE. 100) THEN 

CALL CURVECKFREQC IMirO , RXKMC IMIfO , KXX, -I ) 
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CALL VLIf^ECiniM, IHAX. KFREQ, RXKM) 

ELSE 

CALL CURVE(KFREQ(IMIN),RXKM(inifO,KXX, 0) 

EUDIF 

CALL EMDGR(O) 

CALL EMDPL(O) 

CALL PAGE(11,8.5) 

CALL NOBRDR 
CALL AREA2D(8,6) 

CALL XflAflEC ' FREQ (K) $M00) 

CALL YNAMEC* PHASE (DEG) $*,100) 

CALL HEADIfK’ DFT OUTPUT OF THE RECEIVED SIGNAL $’,100, A, 2) 
CALL HEADIfKPLABEL, 100,5,2) 

CALL GRAFUIMIN, ’SCALE' , HMAX , - 1 8 0 . 0 , A 5 . 0 , 180.0) 

CALL GRIDCI, 1) 

CALL SETCLRC 'MAGENTA' ) 

IF (KXX .LE. 100) THEN 

CALL CURVECKFREQC IMIN),RXKPD( IMIN),KXX, -1 ) 

CALL VLIfJECIMIN, IMAX, KFREQ, RXKPD ) 

ELSE 

CALL CURVECKFREQCIMIN) ,RXKPD(IMIN),KXX, 0) 

ENDIF 

CALL ENDGR(O) 

CALL ENDPL(O) 

800 CONTINUE 
ENDIF 



FORMATS 

1000 F0RMAT(/,2X, 

X 'PLEASE ENTER THE INITIAL POSITION X0,Y0,Z0 (FT) ’,/,2X, 

X ’ OF THE RECEIVER RELATIVE TO THE TRANSMITTER ...’) 

1001 F0RMAT(/,2X, 

'ENTER THE TRAt^$MITTER"S VELOCITY AVERAGE Af^D VAR I ANC E ’ , / , 2X , 

:: ’ IN THE X, Y, Z-DIRECTIONS (FT/SEC) AND (FT/SEC)^)^2 ...',/,2X, 
X ’ (I.E. VXAVG, VXVAR, VYAVG, VYVAR, VZAVG, VZVAR) ’) 

1002 F0RMAT(/,2X, 

X 'ENTER THE AVERAGE SPEED OF SOUND IN FT/SEC ’,/,2X, 

X ’ AfJD THE VARIANCE IN (FT/SEC)Xk2 ...') 

1003 F0RMAT(/,2X, 

X ’ENTER THE TRANSMI TTER"S DONfJ LINK TRANSMITTED ’,/,2X, 

X ’ HALF BEAM NIDTH ANGLE (DEG) ...’) 
lOOA FORMAK/, 2X, ' DO YOU NAfU NORMALIZED AMPLITUDE ’, 

X 'FOR THE RECEIVED 5 1 GNAL ? ' , / , 2X , 

X ' PLEASE EfHER 1:YE5 OR 0 « NO ...') 

1010 F0RMAT(2X, ' ENTER THE ft OF PACKETS IN THE TRAfiSMlTTED SIGflAL ,..') 

1020 F0RMAT(/,2X, 

X 'ENTER THE BAUD TYPE t CORRESPONDING TO THE FOLLOMItlG ’) 

1021 F0RMAT(2X,’ BAUD LEfIGTH FOR PACKET NUMBER: ’,IA,' ...’) 

1022 FORMAT(/, lOX, ’ 1 : BAUD LEtiGTH (DELT) = 1/2A0 SECONDS ’) 

1 023 FORMATdOX,' 2 : BAUD LENGTH (DELT) = 1/120 SECONDS ’) 

102A FORMATdOX,' 3 : BAUD LENGTH (DELT) = 1/60 SECONDS ') 

1025 FORMATdOX,’ A t BAUD LENGTH (DELT) = 1/30 SECONDS ’) 

1026 FORMATdOX,' 5 : BAUD LENGTH (DELT) = 1/15 SECONDS ’) 

1027 F0RMAT(/,2X, 

X 'ENTER THE NUMBER OF BAUDS IN PACKET NUMBER: ’,IA,' ...') 

1030 FORMAK//, 2X. 'MOULD YOU LIKE THE DISCETE fOURIER TRANSFORM ',/, 

X 2X, ' OF THE OUTPUT SIGNAL ? (ENTER 1:YES OR 0:N0) ... ’) 

1060 FORMAK/, 2X, 

X 'MOULD YOU LIKE THE DFT OUTPUT MINDOMED ? ',/,2X, 

X ' ENTER l-.YES OR 0:NO ... ’ ) 

1050 F0RMAT(/,2X, 

X 'PLEASE ENTER THE DESIRED INPUT MIDE BAND ',/,2X, 

X ' SIGNAL-TO-NOISE RATIO IN DB ... ’) 

1100 F0RMAT(/,2X, 

X 'THIS PROGRAM EfiCODES A QPSK MULT I FREQUENCY 5 1 GNAL . ' , / , 2X , 

X ’ THE PHASES ARE SHOMN BELOM FOR ONE FREQUENCY ...’) 

1101 FORMAT(//, 

X 6X, ' ’ ,/, 

X6X,' >' ^ ’/// 
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X 6X, * 

X ex , ' 
X ex , ' 
X ex, * 
X ex , ' 
X ex, * 
X ex , ' 

X 6X, » 

X ex, » 
X ex , ' 

X 6X, * 



1 

, 


t 1 


1 




1 


1 1 






1 


t 1 






1 


f 1 

/ 






1 

, 


t 1 






1 


f » 

, 






I 

t 


1 I 

f 






t 


1 f 

, 






1 

> 


t t 

, 


X 


* ,/, 


1 


1 f 

r 


A 




f 


t 1 

, 




’,///) 



1102 FORMATC//; 2X, 

X ’SELECT OUE OF THE FOLLOWING METHODS FOR ENCODING ’,/,AX, 

X ’THE PHASES FOR EACH FREQUEfiCY FOR EVERY BAUD ’,/,AX, 

X ’WITHIN EACH PACKET . . . ’ ) 

1103 F0RMAT(/,5X, ’ENTER IX, 

X '1: THE PROGRAM RANDOMLY SELECTS ALL THE PHASES ',/,lX, 

X ’ FOR ALL ’,IA,* PACKETS ^,/,lX, 

X ’2: YOU INDIVIDUALLY SELECT THE PHASES ’) 

1110 FORMATC//, 2X, 

X ’SELECT ONE OF THE FOLLOWING METHODS FOR ENCODING ’,/,AX, 

X ’THE PHASES FOR PACKET NUMBER; ’,IA,’ ...’) 

1111 F0RMAT(/,5X, ’EfiTER IX, 

X '1: THE PROGRAM RANDOMLY SELECTS ALL THE PHASES ',/,lX, 

X ’ FOR PACKET NUMBER; ',IA,/,7X, 

X ’2: YOU INDIVIDUALLY SELECT THE PHASES ^,/,lX, 

X ’ FOR PACKET [DUMBER: ’,IA) 

1120 FORMATC//, 2X, 

X ’SELECT ONE OF THE FOLLOWING METHODS FOR ENCODING ’,/,AX, 

X ’THE PHASES FOR PACKET NUMBER; *,IA,’ BAUD [DUMBER; ’,IA,’ ...» 

1121 F0RMAT(/,5X, ’ENTER IX, 

X '1; THE PROGRAM RANDOMLY SELECTS ALL THE PHASES ' ,/,lX, 

X ’ FOR PACKET NUMBER; ',IA,* BAUD NUMBER; ’,IA,/,7X. 

X ’2; YOU INDIVIDUALLY SELECT THE PHASES ',/,lX, 

X ’ FOR PACKET riUMBER; ’;IA,’ BAUD NUMBER; ’,IA) 

1130 FORMATC//, 2X, 

X ’PLEASE ENTER THE DESIRED QUADRANT FOR THE PHASE S/,AX, 

X ’OF PACKET NUMBER: ’,1A,’ BAUD NUMBER: ’,IA,/,AX, 

X ’FREQUEf^CY NUMBER (K): ’,IA,’ ...*) 



C 



2000 FORMATC18) 

2001 F0RMATCA35) 

2010 F0RMATC2X, I6,2X, 16, 2X, FIO .3,2X, lA) 

2020 F0RMAT(//,2X, 18) 

2100 F0RMAT(2X,I6,AX,F1A.A) 

2200 FORMATC//, 2X, ’ DFT INPUT DATA FOR BAUD fl ’,1A,/,5X, 

X ’N’,5X,’REAL P ART ’ , 5X , ’ I MAG PART’) 

2210 F0RMATC2X, lA ,2X, E12.6,2X,E12.6) 

2220 FORMATC//, 2X, ' DFT OUTPUT DATA FOR BAUD P ’,IA,/,5X, 

X'K’ ,7X, ’REAL PART’,7X, ’I MAG PART’,7X, ' MAGNI TUDE ’ , 16 X . ’PHASE',/, 
X 6AX, ’(DEG)’,11X, ’(RAD)’ ) 

22 30 FORMAT (2X, I A , 2X, El A . 8 , 2X, El A . 8 , 2X, El A . 8 , AX, ElA . 8 , 2X, El A . 8 ) 

2300 FORMATC' FOR BAUD NUMBER ’,IA,’ $’) 



PLOTTING 



C 



xxMir^ = 


INT(XXMIN) 


- 


0.2 


XXMAX = 


INT(XXMAX) 


+ 


0 . 2 


YYMIN = 


INT(YYMIfl) 


- 


0.2 


YYMAX = 


IfiTCYYMAX) 


+ 


0.2 


ZZMIN = 


INTCZZMIN) 


- 


0.2 


ZZMAX = 


It^T(ZZMAX) 


+ 


0.2 


CCMIfI = 


INT(CCMIll) 


- 


1 . 0 


CCMAX = 


nJT(CCMAX) 


+ 


1 . 0 


RRMIf^ = 


INT(RRMIN) 


- 


1 .0 


RRMAX = 


irU(RRMAX) 


+ 


1 . 0 


ALFMIN 


= INTCALFMlf]) 


- 0.2 


ALFMAX 


= IMT(ALFMAX) 


+ 0.2 


THTMHi 


= INT(THTMIN) 


- O.A 


THTMAX 


= INT(THTMAX) 


+ O.A 



IF (IDFT .EQ. 0) CALL COMPRS 
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CALL PAGE(11,8.5) 

CALL NODRDR 
CALL AREA2D(8,6) 

CALL XMAME( ' TIME (N) O' , 100) 

CALL YfJAME( ’ MAGNITUDE O'. 100) 

CALL HEADINC RECEIVED SIGNAL $’,100.5.1) 

CALL GRAF(0. ’SCALE’ , r^PTS . RXMI N , ’SCALE’ .RXMAX) 

CALL GRID(l.l) 

CALL SETCLRC ’GREEN’ ) 

CALL CURVECTlME.YYRX.tJPTS, 0) 

CALL ENDGR(O) 

CALL EflDPL(O) 

IF (BDTOTL .GT. 1) THE^J 
CALL PAGE(11.8.5) 

CALL NODRDR 
CALL AREA2D(8.6) 

CALL XfJAMEC* BAUD (LL) $M00) 

XXM = XXMAX - XXMIN 
IF (XXM .GE. 100.0) THEN 
XXMIfJ = XXMlf^ / 1 000 . 0 
XXMAX = XXMAX / 1000.0 
CALL YNAME(’ X(LL) (KFT) $’,100) 

DO 900 1 = 1. BDTOTL 

X(I) = X(I) / 1000.0 
900 CONTIfiUE 

ELSE 

CALL YNAMEC* X(LL) (FT) $’,100) 

ENDIF 

CALL HEADIfK* X-POSITlOfI $’.100,3,1) 

CALL GRAF( 1 . ’SCALE' , BDTOTL , XXMl fJ , ’SCALE’ .XXMAX) 
CALL GRID(1,1) 

CALL SETCLR( ’CYAN’ ) 

CALL CURVE(LBAUD.X. BDTOTL, 0) 

CALL ENDGR(O) 

CALL ENDPL(O) 

CALL PAGE(11,3.5) 

CALL NODRDR 
CALL AREA2D(8,6) 

CALL XfMME(’ BAUD (LL) $’,100) 

YYM = YYMAX - YYMIfJ 
IF (YYM .GE. 100.0) THEN 
YYMIfJ = YYMIfJ / 1 00 0.0 
YYMAX = YYMAX / 1000.0 
CALL YfJAME(’ Y(LL) (KFT) $’,100) 

DO 905 1 = 1, BDTOTL 

Y( I ) = Y( 1 ) / 1 000 . 0 
905 COfJTIfJUE 

ELSE 

CALL YfJAMEC Y(LL) (FT) $’,100) 

EfJDIF 

CALL HEADINC Y-POSITIOfJ $’.100,5.1) 

CALL GRAFd, ’SCALE’ , DDT OTL .YYMIfJ, ’ SCAL E ’. YYMAX ) 
CALL GRID( 1 , 1 ) 

CALL SETCLR( 'MAGEfJTA ' ) 

CALL CURVE(LBAUD,Y, BDTOTL. 0) 

CALL EfJDGR(O) 

CALL EfJDPL(O) 

CALL PAGE(11.8.5) 

CALL NODRDR 
CALL AREA2D(8,6) 

CALL XfJAME( ’ BAUD (LL ) $ ' , 100) 

ZZM = ZZMAX - ZZMIfJ 
IF (ZZM .GE. 100.0) THEfJ 
ZZMIfJ = ZZMIfI / 1 000.0 
ZZMAX = ZZMAX / 1000.0 
CALL YfJAME(’ Z(LL) (KFT) $’,100) 

DO 910 I = 1. BDTOTL 
Z(I) = Z(I) / 1000.0 
910 COfJTIfJUE 

ELSE 

CALL YtJAME( ’ Z(LL ) (FT) $ ’ , 100) 

EfJDIF 
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CALL HEADIfiC’ Z-POSITIOf^ $S100,3,1) 

CALL GRAFCl, ' SCA L E * , B DT OT L , ZZhI ’ SCAL E ’ , ZZMAX ) 

CALL GRIDCl, 1) 

CALL SETCLRC ' BLUE' ) 

CALL CURVE(LBAUD,Z.DDTOTL,0) 

CALL ENDGRCO) 

CALL ENDPL(O) 

CALL PAGE(11,8.5) 

CALL MODRDR 
CALL AREA2D(8,6) 

CALL XNAMEC * BAUD (LL) $MOO) 
ccn = ccMAX - ccMin 
IF (CCMAX .GE. 100.0) THEf^ 

CCMIfi = CCMIN / 1 000.0 
CCMAX = CCMAX / 1000.0 

CALL YflAMEC' C(LL) (KFT/SEC) $M00) 

DO 915 I = 1, BDTOTL 
C( I ) = C( I ) / 1 000 . 0 
915 COtiTlfiUE 

ELSE 

CALL YNAMEC ' C(LL) (FT/SEC) $M00) 

ENDIF 

CALL HEADIIU * 'SPEED OF SOUND $M00,3,1) 

CALL GRAFCl, * SCALE' , BDTOT L , CCMI N , 'SCALE' , CCMAX) 

CALL GRID(1,1) 

CALL SETCLRC 'RED' ) 

CALL CURVECLBAUD,C, BDTOTL, 0) 

CALL EtJDGRCO) 

CALL ENDPLCO) 

CALL PAGEC11,8.5) 

CALL NOBRDR 
CALL AREA2DC8,6) 

CALL XNAMEC' BAUD CLL) $*,100) 

RRM = RRMAX - RRMIt^ 

IF (RRM .GE. 100.0) THEN 
RRMIfJ = RRMIN / 1 000.0 
RRMAX = RRMAX / 1000.0 
CALL YfiAMEC RCLL) (KFT) $*,100) 

DO 920 I = 1, BDTOTL 
RCI) = RCI) / 1000.0 
920 CONTItiUE 

ELSE 

CALL YNAMEC* RCLL) (FT) $',100) 

ENDIF 

CALL HEADINC* SLANT RANGE TO RECEIVER $*,100,3,1) 

CALL GRAFC 1 , ’ SCALE' , BDTOTL , RRMIN, 'SCALE' , RRMAX) 

CALL GRIDCl,!) 

CALL SEICLRC 'GREEN* ) 

CALL CURVE(LBAUD,R, BDTOTL, 0) 

CALL ENDGRCO) 

CALL EfJDPLCO) 

CALL PAGEC11,5.5) 

CALL NOBRDR 
CALL AREA2D(8,6) 

CALL XtlAMEC * BAUD CLL) $* ,100) 

CALL YNAMEC' ALPHACLL) $',100) 

CALL HEADINC 'COMPRESSION FACTOR DUE TO THE MOVIflG TX $*,100,3,1) 
CALL GRAFCl, 'SCALE* , B DT OT L , A L FMI fJ , 'SCALE* ,ALFMAX) 

CALL GRIDCl, 1) 

CALL SETCLRC *RED* ) 

CALL CURVECLBAUD, ALPHA, BDTOTL, 0) 

CALL ENDGRCO) 

CALL EliDPLCO) 

CALL PAGEC11,8.5) 

CALL liOBRDR 
CALL AREA2DC8,6) 

CALL XNAMEC BAUD CLL) $*,100) 

CALL YNAMEC* THETACLL) (DEG) $*,100) 

CALL HEADINC *ANGLE BETHEEN RCLL) AND ZO $*,100,3,1) 

CALL GRAFCl, * SCAL E * , BDT OT L , THTMI * SCAL E * , T HTMAX!) 

CALL GRIDCl,!) 

CALL SETCLRC ' BLUE* ) 
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CALL CURVEC L BAUDOT HET AD, BDT OIL, 0) 
CALL ENDGR(O) 

CALL ENDPL(O) 

EN DIF 

CALL DONEPL 



9999 STOP 
END 



C 



C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 



SUBROUTINE GAUSSC DSEED, AVG, VAR, ZRND) 

THIS SUBROUTINE GENERATES A GAUSSIAN RANDOM NUMBER HITH 
MEAN = AVG AND VARIANCE = VAR 

INPUTS: DSEED = A DOUBLE PRECISIOf^ SEED THAT MUST BE A VARIABLE 

AVG = THE MEAN OF THE GAUSSIAN RANDOM VARIABLE 
VAR = THE VARlAfJCE OF THE GAUSSIAN RAfJDOM VARIABLE 

OUTPUTS: DSEED = THE SEED IS CHANCED DURING EXECUTION NHICH IS 

REQUIRED FOR THE NEXT SUBROUTINE CALL 
ZRflD = A GAUSSIAri RAfJDOM NUMBER HITH MEAfJ = AVG AND 
VARIANCE = VAR 



REAL AVG, VAR, ZRND 

DOUBLE PRECISION ZZ, DSEED, URND 

ZZ = 0.0 

DO 5000 I = 1, 12 

DSEED = DSEED + 135795^5.1597537 
URND = GGUBFS(DSEED) 

ZZ = (URtJD - 0.5) + ZZ 
5000 CONTINUE 

ZRND = (ZZ (VARX)(0.5)) + AVG 

RETURN 

END 



C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 



SUBROUT I rJE PHASEC DSEED, RNDPHI , IPHI ) 

THIS SUBROUTINE SELECTS A PHASE RANDOMLY FROM QUADRAfUS 1 TO A 

irrPUTS: DSEED = A DOUBLE PREClSIOfi SEED THAT MUST BE A VARIABLE 

OUTPUTS: DSEED = THE SEED IS CHANGED DURING EXECUTION MHICH IS 
REQUIRED FOR THE NEXT SUBROUTINE CALL 
RfJDPHI = A RANDOM PHASE FROM ONE OF THE FOUR QUADRANTS 
IN RADIANS 

IPHI = THE QUADRANT NUMBER OF THE RAfiDOM PHASE 



C 



C 



REAL PHIRND, RfJDPHI 
DOUBLE PRECISIOfl DSEED 
irJTEGER IPHI 



PI = A.O^ATANCl .0) 
PHIRflD = 0.0 
IPHI = 0 

PHIRND = GGUBFSC DSEED) 
DSEED = DSEED + 127 . A55 
PHIRND = (PHIRND > A.O) 
IPHI = INT(PHIRND) + 1 



IF 


(IPHI 


.EQ. 


1 ) 


RNDPHI 


IF 


(IPHI 


.EQ. 


2) 


RNDPHI 


IF 


(IPHI 


.EQ. 


5) 


RNDPHI 


IF 


(IPHI 


• EQ. 


A) 


RNDPHI 


RET 


URN 








END 









(A5.0 ^ PI) / 180.0 
(155.0 ^ PI) / 180.0 
(-155.0 PI) / 180.0 
( -A5 . 0 X PI ) / 180.0 



SUBROUTINE D FT ( NMl , XI N , XOUT ) 

C 

C THIS SUBROUTINE COMPUTES THE DISCRETE FOURIER TRANSFORM OF 
C A COMPLEX DATA SET OF N (= NMl+1) POIfHS STORED IN THE ARRAY XIN 
C THE RESULT IS STORED IN THE COMPLEX ARRAY XOUT 
C 

COMPLEX XIN(0:NM1), XOUT(0:NM1), H, MM 
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o o n 



N = flMl + 1 

PI = ^<.0>^ATAN(1 .0) 

EtI = U 

NM2 = mil - 1 
IF (MMl .EQ. 0) THEN 
XOUT(O) = XIN(O) 

ELSE 

BETA = 2.0 ^ PI / Er^ 

H = CMPLX(COS( BETA), “SIN(BETA)) 

DO 3110 K = 0, NMl 

m - 

XOUT(K) = XIN(NMl) 

DO 3100 L = NM2, 0, -1 

XOUT(K) = XOUT(K)xhM + XltUL) 
5100 CONTIflOE 

5110 CONTINUE 
C 

ENDIF 

RETURN 

END 

C 

SUBROUT I flE VLINEC Jl, JN,XARRAY, B) 

REAL XARRAY(0rA096), B(0:A096) 

DRAH VERTICAL LINES 

DO 5200 JLINE = Jl, JN 
XFROM = XARRAYC JLINE) 

XTO = XFROM 
YFROM = 0.0 
YTO = BCJLItlE) 

CALL RLVEC(XFROM,YFROM, XTO, YTO, 0000) 
3200 CONTINUE 
RETURN 
END 
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